
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

72180
10.1038/s41598-024-72180-x
Article
Differential effects of environmental exposures on clinically relevant endophenotypes between sexes
González Zarzar Tomás tgonza@uc.cl

12
Palmiero Nicole E. 3
Kim Dokyoon 4
Shen Li 4
Hall Molly A. 3
1 https://ror.org/04teye511 grid.7870.8 0000 0001 2157 0406 School of Medicine, Pontificia Universidad Católica de Chile, Santiago, Chile
2 https://ror.org/04teye511 grid.7870.8 0000 0001 2157 0406 Institute for Biological and Medical Engineering, Pontificia Universidad Católica de Chile, Santiago, Chile
3 grid.25879.31 0000 0004 1936 8972 Institute for Biomedical Informatics, Department of Genetics, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA 19104 USA
4 grid.25879.31 0000 0004 1936 8972 Department of Biostatistics, Epidemiology and Informatics, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA 19104 USA
13 9 2024
13 9 2024
2024
14 2145323 1 2024
4 9 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
Sex and gender differences play a crucial role in health and disease outcomes. This study used data from the National Health and Nutrition Examination Survey to explore how environmental exposures affect health-related traits differently in males and females. We utilized a sex-stratified phenomic environment-wide association study (PheEWAS), which allowed the identification of associations across a wide range of phenotypes and environmental exposures. We examined associations between 272 environmental exposures, including smoking-related exposures such as cotinine levels and smoking habits, and 58 clinically relevant blood phenotypes, such as serum albumin and homocysteine levels. Our analysis identified 119 sex-specific associations. For example, smoking-related exposures had a stronger impact on increasing homocysteine, hemoglobin, and hematocrit levels in females while reducing serum albumin and bilirubin levels and increasing c-reactive protein levels more significantly in males. These findings suggest mechanisms by which smoking exposure may pose higher cardiovascular risks and greater induced hypoxia for women, and greater inflammatory and immune responses in men. The results highlight the importance of considering sex differences in biomedical research. Understanding these differences can help develop more personalized and effective health interventions and improve clinical outcomes for both men and women.

Keywords

PheEWAS
Environmental exposures
Smoking exposure
Albumin
Homocysteine
Red blood cell
Inflammation
CVD
Subject terms

Biological techniques
Biomarkers
Medical research
Risk factors
issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Sex and gender differences are crucial in health and disease outcomes1. These differences in disease prevalence, risk factors, progression, and treatment response underscore the importance of considering sex and gender differences in biomedical research. For instance, females are at greater risk of developing immune diseases, while males have a higher risk of non-reproductive cancers2. Additionally, cardiovascular diseases have different pathologic mechanisms across sexes3, and females are more susceptible to chronic obstructive pulmonary disease (COPD)4.

While some sex and gender differences are well-documented across various health outcomes, significant gaps remain in our understanding of how environmental factors contribute to these disparities. This gap is particularly evident in research on environmental influences on health, which has yet to explore how these exposures differentially affect males and females extensively5. For example, sex differences in response to environmental exposures, such as those observed in occupational asthma and smoking-related lung function decline, highlight the need for more comprehensive studies in this area6–10.

Some sex differences, such as the difference in occupational asthma, might be explained by different degrees of environmental exposure between sexes, which might be attributed to gender differences6. In contrast, other differences, such as mortality and lung function differences between male and female smokers, are due to differential effects of the same environmental exposure10. Differential biological pathways, interactions with sex hormones, or genetic and epigenetic mechanisms might mediate the latter type of differences11.

Sex refers to the biological attributes that distinguish male and female organisms, including chromosomes, hormone levels, and reproductive anatomy. These biological differences are typically categorized as male or female but include intersex. Gender, on the other hand, refers to the roles, behaviors, activities, expectations, and norms that cultures and societies consider appropriate for men, women, and gender-diverse people. It is crucial to understand the distinction and the potential for interaction between sex and gender in biomedical research, as both can influence health outcomes and provide different interpretative scopes1,12.

Adopting research methodologies that capture the complex interplay between sex, environmental exposures, and health outcomes is essential to addressing these gaps. Genome-wide association studies (GWAS) and phenome-wide association studies (PheWAS) have been pivotal in identifying genetic and phenotypic associations, respectively, across diverse populations13–20. Environment-wide association studies (EWAS) have expanded on these methodologies and explored the associations between numerous environmental exposures and health outcomes21–25. While evaluating associations between multiple exposures and a single phenotype in EWAS has proven valuable, considering phenotypes in isolation can hinder the ability to detect exposures with cross-phenotype associations or associations between a single environmental exposure and multiple phenotypes.

Understanding the sex-specific effects of environmental exposures is crucial for assessing environmental risks, guiding therapeutic recommendations, and generating new biological hypotheses. To address this knowledge gap, we conducted a sex-stratified phenomic environment-wide association study (PheEWAS) to investigate the associations between environmental exposures and clinically relevant phenotypes. PheEWAS is a hypothesis-free approach that tests for significant associations between multiple environmental exposures and multiple phenotypes. By examining these associations across a broad range of phenotypes, this study aims to uncover sex-specific mechanisms that could inform precision health strategies and improve clinical outcomes. Utilizing data from the National Health and Nutrition Examination Survey (NHANES), we evaluated 18,937 participants across 58 phenotypes and 272 environmental exposures. By classifying sex-different exposure-phenotype associations into pure, quantitative, and qualitative categories26, we aim to elucidate the mechanisms underlying these associations and provide insights that could inform precision health approaches. Our results suggest promising avenues to explore exposures that pose differential risks to the sexes, elucidate biological pathways with significant sex-based differences, and open the road to more personalized therapeutics.

Materials and methods

NHANES data sets

This study utilized data from NHANES to identify sex-specific exposure-phenotype associations. NHANES is a program of studies that uses data from physical examinations and interviews to assess the health and nutritional status of adults and children in the U.S. The physical examination includes medical, dental, and physiological measurements and laboratory tests. The interview includes demographic, socioeconomic, dietary, and health-related questions.

We analyzed 1191 variables and 41,474 participants from four survey periods (1999–2000, 2001–2002, 2003–2004, and 2005–2006),27,28, using a discovery (1999–2000 and 2001–2002) and replication (2003–2004 and 2005–2006) design29,30, split by sex into four cohorts: discovery females (df), discovery males (dm), replication females (rf), and replication males (rm).

To select relevant phenotypes and exposures, we focused on 84 routine clinical measurements, including blood and urine biomarkers and hormones, and 739 exposures related to health outcomes, such as tobacco and alcohol use, heavy metals, and volatile compounds. PheEWAS, as a hypothesis-generating approach, relies on a wide selection of exposures and phenotypes, only restricted to retaining high-quality variables. Data cleaning procedures described below ensured high-quality phenotypes and exposures.

We applied survey weights according to NHANES guidelines, using 4-year weights for the discovery cohorts and adjusted 2-year weights for the replication cohorts. We excluded three phenotypes and 116 exposures lacking survey weight data. Data cleaning, QC, and analyses were done using CLARITE software30. All analytical steps and code are available in Zenodo (DOI: 10.5281/zenodo.6483503), along with replication instructions and package details.

Data quality control

The QC process was conducted independently in each of the four cohorts, but if a variable was removed in one cohort, it was removed from all. Participants aged 18 and older, with no upper age limit, and complete covariate data were included (N=18,937,Ndf=4724,Ndm=4339,Nrf=5123,Nrm=4751). Low-density lipoprotein values (NHANES code name: LBDLDL) were adjusted by dividing by 0.7 for participants on statins.

Variables were categorized using CLARITE as constant, binary, categorical, and continuous. Constants (single-value variables that were not covariates or phenotypes), variables with only missing data, and continuous variables with less than ten unique values were removed. Filters ensured proper sample sizes and high-quality variables in each cohort. Variables with fewer than 200 non-missing values, binary or categorical variables with fewer than 200 observations per category, and continuous variables with more than 90% of values at zero were excluded. This left 59 phenotypes and 281 exposures.

Continuous variables were log2 transformed and z-score normalized to address skewed distributions while retaining biologically meaningful values. Briefly, transforming raw values into their base 2 logarithms provides a way to produce symmetric distributions and keep easily interpretable values because changes in a one-unit increase in a log2-transformed variable correspond to a doubling of the original value31,32. On the other hand, z-score normalization makes variables comparable even when measured on different scales or units, allowing comparisons across exposure-phenotype results. Outliers, defined as values three standard deviations from the mean, were treated as missing and re-filtered, removing nine additional exposures. A Pearson correlation was used to detect identical transformed variables, leading to the removal of one phenotype. The final dataset for PheEWAS included eight covariates, 58 phenotypes, and 272 exposures. Summary statistics for all covariates are available in Table 1.Table 1 Descriptive statistics.

Cohort	Sex	N	Age	BMI	SES	Black	Mexican	Other Hispanic	Other ethnicity	
Discovery	Females	4724	44.69 (19.72)	28.50 (6.87)	0 = 1926 1 = 1557 2 = 1241	911	1194	253	170	
Males	4339	46.21 (19.82)	27.45 (5.47)	0 = 1565 1 = 1515 2 = 1259	840	1091	200	144	
Replication	Females	5123	45.09 (20.34)	28.70 (7.13)	0 = 2349 1 = 1497 2 = 1277	1175	1065	160	227	
Males	4751	46.58 (20.29)	27.95 (5.86)	0 = 1935 1 = 1525 2 = 1291	1089	989	139	179	

Statistical methods

PheEWAS

Following the preparation and cleaning of the dataset, statistical analyses were conducted to explore the associations between exposures and phenotypes. A sex-stratified PheEWAS was conducted to estimate associations between 58 phenotypes and 272 exposures, adjusting for eight covariates, including age, socioeconomic status (SES), body mass index (BMI), survey year, and four race/ethnicity categories. Analyses were performed separately for each cohort using generalized linear models adjusted for survey weights. A minimum of 200 complete observations were required for a successful regression. Finally, observations were dropped in each regression if they had missing values in the phenotype or exposure or missing survey weight information. An inverse variance meta-analysis was used to generate sex-specific results in the female and male cohorts26,33.

Sex difference analysis

Sex differences in exposure-phenotype associations were assessed by comparing β values between males and females. Differences were categorized as pure (significant in one sex but not the other), quantitative (significant in both sexes with differing effect sizes), or qualitative (significant in both sexes with opposing effect directions). Two approaches were used to detect sex differences:Full Evaluation Approach: all possible exposure–phenotype associations were tested for sex differences using the following formula 1 Zdiff=βf-βmsef2+sem2

where βf and βm are the effect sizes in the female and male cohorts, respectively, and sef and sem are their standard errors. Significant sex differences were identified using a Bonferroni-corrected αdiff based on M total number of tests, such that αdiff=0.05M.

Filtering-First Approach: associations significant across both sexes were filtered first using the following formula 2 Zoverall=βfsef2+βmsem21sef2+1sem2

Associations with poverall<10-5 were retained for sex difference testing using formula 1. A Bonferroni-corrected αdiff was applied based on N number of filtered tests, such that αdiff=0.05N.

Both approaches were employed to capture different types of sex differences, providing a comprehensive evaluation. Specifically, the full evaluation approach is better suited to detect qualitative sex differences, whereas the filtering-first approach is better suited to detect quantitative and pure sex differences26. Differences significant by either approach were classified into pure, quantitative, or qualitative categories.

Age stratification

To assess the robustness of our findings and account for potential confounding effects of hormonal changes across the lifespan, we replicated the analysis by excluding participants 46 years or older. We retained 10,149 participants, with sample sizes across cohorts of Ndf=2650;Ndm=2248;Nrf=2822;Nrm=2429. This exclusion criterion removes the influence of age-related hormonal transitions on the observed associations between environmental exposures and health outcomes.

Results

We conducted 15,999 association tests in each of the four NHANES cohorts: discovery females, discovery males, replication females, and replication males. For further details on the sample sizes of each cohort, please refer to the “Materials and methods” section. The number of converged tests varied slightly across cohorts because some tests had fewer than 200 complete observations (between 15,470 in the replication males cohort and 15,994 in the discovery females and discovery males cohorts). From the successful tests, 15,466 exposure-phenotype associations converged across all cohorts and were used to estimate sex differences. We identified 119 significant sex-specific associations, of which 49 were classified as pure, 45 as quantitative, and 25 as qualitative.

Among the 58 phenotypes evaluated, 35 showed sex-specific associations with at least one exposure. Serum albumin levels had the highest number of sex-specific associations, with 30 different exposures, including smoking-related exposures (e.g., blood cotinine levels, number of cigarettes smoked daily, and whether anyone smokes in the household), nutrients (e.g., folate, vitamins A, D, and E), supplements (e.g., folic acid, calcium, vitamin C), and heavy metals and volatile compounds. Other phenotypes with notable associations included red cell count and C-reactive protein levels, which had sex-specific associations with heavy metals like cadmium and lead, and volatile compounds such as ethylbenzene and toluene (Figs. 1, 2).Figure 1 Miami Plot from PheEWAS results stratified by sex. The Y axis shows the -log of the p-value, and the X axis is sorted by exposure categories. Significant sex differences are highlighted in color, and categorized as pure, quantitative, and qualitative sex differences. Relevant exposure-phenotype associations discussed in the manuscript are annotated.

Of the 272 exposures evaluated, 53 had sex-specific associations with at least one phenotype. Exposures to heavy metals (e.g., lead and cadmium), nutrients (e.g., vitamins A and E, lutein and zeaxanthin), and smoking behaviors were most frequently associated with phenotypes such as serum albumin, C-reactive protein, and red cell count. Lead exposure, for example, was associated with 11 different phenotypes, underscoring its widespread impact on various clinically relevant measurements (Figs. 1, 2). Supplementary Table 1 shows the complete set of results in more detail.Figure 2 Circos plot of significant exposure-phenotype associations divided by sex and type of sex difference. In each circos plot, exposures are in the upper section and phenotypes in the lower section of the circle. The presence of a line connecting an exposure and phenotype indicates a significant association. Red lines indicate a negative association and blue lines indicate a positive association. The width of the line indicates the effect size, with greater width indicating a greater absolute effect size.

Pure sex differences

We identified 53 pure sex-specific associations, with a slight bias toward males. In females, exposures in pure associations were primarily heavy metals, food component recall, nutrients, and volatile compounds, associated with phenotypes like red cell count, lymphocyte percent, albumin, HDL cholesterol, and homocysteine. For example, blood levels of the volatile compounds ethylbenzene (βf=0.096,βm=-0.020,Pdiff=1.53E-05) and toluene (βf=0.115,βm=-0.015,Pdiff=3.91E-06) were both associated with red cell counts in females. In contrast, lead was associated with albumin (βf=0.084,βm=-0.028,Pdiff=4.31E-09), creatinine (βf=0.076,βm=-0.005,Pdiff=1.96E-06), and sodium levels (βf=0.07,βm=-0.024,Pdiff=4.49E-06). In males, exposures in pure associations were often categorized as nutrients, smoking behavior, and supplement use, with albumin showing associations with 17 different exposures (Figs. 2, 3; Table 2).Table 2 Pure sex differences results.

Exposure name	Outcome name	βf	Pf	βm	Pm	Pdiff	
Cotinine	Albumin	0.005	4.41E−01	− 0.089	6.80E−16	2.36E−13	
Cotinine	Osmolality	− 0.015	1.23E−01	− 0.086	9.96E−12	8.16E−06	
Cotinine	Red cell count	0.074	4.25E−09	− 0.006	5.96E−01	3.90E−06	
Alcohol	Cholesterol, total	0.0	9.79E−01	0.074	3.69E−07	1.67E−05	
Alcohol	Total cholesterol	0.001	9.48E−01	0.075	4.29E−07	1.77E−05	
Calcium	Albumin	− 0.061	1.85E−12	0.009	4.40E−01	2.36E−06	
Carbohydrate	Lymphocyte percent	− 0.063	3.96E−08	0.026	5.76E−02	7.06E−07	
PFA 22:6 (Docosahexaenoic)	Chloride	− 0.062	1.57E−07	0.018	1.88E−01	8.34E−06	
Cadmium	Albumin	− 0.013	1.65E−01	− 0.098	2.18E−17	6.12E−09	
Cobalt, urine	Hematocrit	− 0.141	2.29E−14	0.013	5.47E−01	7.75E−08	
Cobalt, urine	Hemoglobin	− 0.155	2.32E−18	0.005	8.24E−01	2.40E−08	
Lead	1/Creatinine	− 0.076	3.09E−08	0.005	6.00E−01	1.96E−06	
Lead	Albumin	0.083	3.63E−14	−0.022	1.24E−01	4.31E−09	
Lead	Direct HDL-Cholesterol	− 0.019	2.33E−01	0.095	7.80E−07	4.95E−06	
Lead	Sodium	0.07	4.30E−07	− 0.024	1.12E−01	4.49E−06	
Alpha-carotene	Osmolality	− 0.015	1.85E−01	0.056	1.13E−06	1.01E−05	
Folate, RBC	Bilirubin, total	0.011	3.69E−01	0.097	1.96E−16	4.46E−07	
Lutein and zeaxanthin	Albumin	0.013	3.25E−01	0.123	1.22E−20	1.75E−09	
Lutein and zeaxanthin	Mean cell volume	0.024	5.76E−02	− 0.063	4.42E−10	7.05E−08	
Vitamin A	C-reactive protein	0.014	2.80E−01	− 0.179	6.20E−30	7.25E−21	
Vitamin A	Lymphocyte number	0.099	8.29E−10	− 0.013	3.53E−01	2.03E−07	
Vitamin A	Lymphocyte percent	0.066	1.95E−06	− 0.022	1.23E−01	9.64E−06	
Vitamin D	Albumin	− 0.005	5.98E−01	0.083	6.37E−11	5.63E−08	
Vitamin D	Triglycerides	0.032	7.73E−02	− 0.066	3.56E−07	1.16E−05	
Vitamin E	Albumin	− 0.018	2.81E−01	0.16	1.17E−32	4.02E−17	
Vitamin E	Direct HDL-Cholesterol	0.083	1.77E−07	− 0.031	6.29E−02	7.05E−07	
Vitamin E	Protein, total	− 0.014	3.71E−01	0.088	1.16E−10	7.56E−07	
Vitamin E	Red cell count SI	− 0.028	7.60E−02	0.066	1.05E−08	1.65E−06	
b-Cryptoxanthin	Mean cell hemoglobin	0.003	8.35E−01	− 0.071	4.61E−09	1.89E−05	
# cigarettes smoked per day now	Albumin	0.019	1.98E−01	− 0.121	1.14E−09	1.30E−08	
# cigarettes smoked per day now	Homocysteine	0.152	7.87E−14	0.027	1.69E−01	7.91E−06	
# days smoked cigs during past 30 days	Albumin	0.022	8.40E−02	− 0.123	9.99E−12	5.51E−11	
# days smoked cigs during past 30 days	Homocysteine	0.147	7.18E−15	0.023	2.08E−01	2.33E−06	
Avg # cigarettes/day during past 30 days	Albumin	0.018	1.69E−01	− 0.109	8.91E−09	3.98E−08	
Current Cigarette Smoker?	C-reactive protein	0.061	7.15E−02	0.248	1.05E−21	1.21E−05	
Current or Past Cigarette Smoker?	Albumin	0.012	4.82E−01	− 0.187	4.29E−11	2.25E−09	
FTC Nicotine Content	Albumin	0.018	1.99E−01	− 0.086	1.06E−06	3.72E−06	
FTC Tar Content	Albumin	0.017	2.30E−01	− 0.094	4.60E−08	5.74E−07	
How many years smoked this amount	Albumin	0.005	6.76E−01	− 0.096	1.54E−07	3.36E−06	
Does anyone smoke in home?	Albumin	− 0.02	3.11E−01	− 0.179	2.25E−12	1.06E−06	
Does anyone smoke in home?	Direct HDL-Cholesterol	− 0.28	5.27E−20	− 0.051	1.89E−01	3.38E−06	
Does anyone smoke in home?	Red cell count SI	0.186	2.49E−09	0.002	9.42E−01	4.52E−06	
Calcium	Albumin	− 0.012	2.08E−01	0.06	1.84E−06	4.67E−06	
Copper	Albumin	− 0.011	2.48E−01	0.064	1.24E−07	1.35E−06	
Vitamin B12	Albumin	− 0.01	2.99E−01	0.063	5.69E−07	3.51E−06	
Vitamin C	Albumin	− 0.003	7.83E−01	0.07	4.09E−09	4.43E−06	
Blood Benzene	1/Creatinine	0.015	4.13E−01	0.111	1.60E−16	1.51E−05	
Blood Ethylbenzene	Red cell count SI	0.096	5.62E−07	− 0.02	2.91E−01	1.54E−05	
Blood Toluene	Red cell count SI	0.115	1.90E−09	− 0.015	4.80E−01	3.91E−06	
Weighted generalized linear model results, indicating the environmental exposure (predictor), the phenotype (outcome), and the sex-stratified β and P values, and the P values of the sex difference test. For more detailed information on the results the Supplementary Table 1 shows the complete set of results in each cohort, including the standard errors, and units of measurement

Nineteen pure associations involved smoking-related exposures. For example, cotinine levels were associated with albumin levels (βf=0.005,βm=-0.089,Pdiff=2.36E-13) and osmolality (βf=-0.015,βm=-0.086,Pdiff=8.16E-06) in males and with red cell count (βf=0.074,βm=-0.006,Pdiff=3.9E-06) in females. Similarly, “Does anyone smoke in home?” was associated with albumin levels (βf=-0.020,βm=-0.179;Pdiff=0.1E-06) in males and with HDL cholesterol (βf=-0.280,βm=-0.051,Pdiff=3.38E-06) and red cell count (βf=0.186,βm=0.002,Pdiff=4.52E-06) in females. Two of the ten survey questions related to smoking behavior were positively associated with homocysteine levels in females. In contrast, seven were associated with albumin and one with C-reactive protein in males (Fig. 3).Figure 3 Forest plot showing the effect sizes and 95% confidence intervals in males and females. Exposures are defined on the left side and phenotypes on the right. Plots are divided into the three sex difference categories of pure, quantitative, and qualitative.

Quantitative sex differences

Forty-five quantitative sex-specific associations were identified, with greater effect sizes observed more frequently in males. Notably, heavy metal exposures like cadmium had greater effect on males on phenotypes like monocyte number (βf=0.092,βm=0.166,Pdiff=7.48E-06), segmented neutrophils number (βf=0.111,βm=0.206,Pdiff=2.2E-08), C-reactive protein level (βf=0.045,βm=0.144,Pdiff=6.01E-09), mean cell hemoglobin level (βf=0.047,βm=0.135,Pdiff=1.01E-09), mean cell volume (βf=0.057,βm=0.145,Pdiff=1.27E-10), and white blood cell count (βf=0.151,βm=0.252,Pdiff=9.7E-10); while in females it had greater effect on red cell count (βf=0.127,βm=0.027,Pdiff=6.97E-11). On the other hand, lead exposure had greater effect on females on bone alkaline phosphatase level (βf=0.233,βm=0.049,Pdiff=5.87E-14), homocysteine level (βf=0.210,βm=0.130,Pdiff=5.17E-06), red cell count (βf=0.190,βm=0.071,Pdiff=1.1E-09), and alkaline phosphatase level (βf=0.193,βm=0.064,Pdiff=7.9E-12) (Figs. 2, 3; Table 3).

Smoking-related exposures, including cotinine levels, “Does anyone smoke in home?”, and “Current or Past Cigarette Smoker?” had greater positive effects on hematocrit in females. Similarly, associations between current or past cigarette smoking and hemoglobin level (βf=0.429,βm=0.236,Pdiff=1.25E-06) and between the average number of cigarettes smoked per day during the past 30 days and homocysteine level (βf=0.154,βm=0.038,Pdiff=2.71E-06) had greater effects in females. Conversely, all other quantitative associations between smoking-related exposures and bilirubin level, creatinine level, segmented neutrophils number, mean cell hemoglobin level, and white blood cell count had greater effects in males (Fig. 3).

Quantitative associations involving serum albumin, triglyceride, and GGT levels tended to have larger effects in males. For example, the associations between g-tocopherol and GGT level (βf=0.066,βm=0.129,Pdiff=2.08E-05), g-tocopherol and triglyceride level (βf=0.20,βm=0.319,Pdiff=1.35E-10), retinyl palmitate and triglyceride level (βf=0.273,βm=0.394,Pdiff=7.65E-11), vitamin A and albumin level (βf=0.145,βm=0.258,Pdiff=8.65E-11), and vitamin E and GGT level (βf=0.032,βm=0.145,Pdiff=4.3E-07) had greater effects sizes in males. Conversely, associations with homocysteine and urine albumin tended to have greater effects in females than in males. For example, the associations between iron supplement use and homocysteine level (βf=-0.127,βm=-0.064,Pdiff=1.78E-07) and between cesium level in urine and urine albumin level (βf=0.504,βm=0.362,Pdiff=2.74E-09) had greater effect sizes in females.Table 3 Quantitative sex differences results.

Exposure name	Outcome name	βf	Pf	βm	Pm	Pdiff	
Cotinine	Hematocrit	0.163	8.81E−41	0.09	3.81E−14	1.65E−05	
Alcohol	Direct HDL-Cholesterol	0.114	5.20E−16	0.214	2.75E−32	1.18E−05	
Alcohol	GGT	0.063	7.31E−09	0.142	3.33E−39	3.46E−07	
Alcohol	Iron	0.067	6.19E−09	0.151	1.47E−26	4.65E−06	
Cadmium	C-reactive protein	0.045	1.77E−05	0.144	1.75E−26	6.01E−09	
Cadmium	Mean cell hemoglobin	0.047	2.61E−07	0.135	3.23E−33	1.01E−09	
Cadmium	Mean cell volume	0.057	4.26E−09	0.145	1.12E−51	1.26E−10	
Cadmium	Monocyte number	0.092	4.60E−17	0.166	7.74E−41	7.48E−06	
Cadmium	Red cell count SI	0.127	2.53E−30	0.027	1.22E−02	6.97E−11	
Cadmium	Segmented neutrophils number	0.111	9.17E−20	0.205	1.16E−69	2.20E−08	
Cadmium	White blood cell count (SI)	0.151	9.32E−36	0.252	1.82E−112	9.70E−10	
Cesium, urine	Albumin, urineSI	0.511	1.02E−192	0.38	2.96E−59	7.22E−06	
Cesium, urine	Albumin, urine	0.504	5.45E−245	0.362	8.58E−85	2.74E−09	
Lead	Alkaline phosphotase	0.193	3.84E−53	0.064	4.85E−06	7.90E−12	
Lead	Bone alkaline phosphotase	0.233	1.44E−40	0.049	3.97E−03	5.87E−14	
Lead	Homocysteine	0.21	4.66E−60	0.13	2.82E−27	5.17E−06	
Lead	Red cell count SI	0.19	8.14E−36	0.07	1.52E−08	1.10E−09	
Thallium, urine	Albumin, urine SI	0.457	4.26E−121	0.325	6.08E−45	1.28E−05	
Thallium, urine	Albumin, urine	0.447	4.64E−144	0.333	7.42E−64	1.57E−05	
Folate, serum	Albumin	0.058	5.57E−08	0.142	1.07E−35	8.61E−08	
Lutein and zeaxanthin	Direct HDL-Cholesterol (mg/dL)	0.277	5.63E−61	0.16	1.93E−16	5.32E−06	
Retinyl Palmitate	Triglyceride	0.183	2.57E−21	0.338	1.97E−62	3.10E−08	
Retinyl Palmitate	Triglycerides	0.273	2.11E−75	0.394	6.51E−262	7.65E−11	
Retinyl Stearate	Triglycerides	0.287	2.13E−69	0.407	1.16E−176	3.17E−08	
Vitamin A	Albumin	0.145	1.50E−37	0.258	1.98E−82	8.65E−11	
Vitamin A	Blood urea nitrogen	0.261	1.15E−106	0.174	5.16E−27	1.45E−05	
Vitamin E	GGT	0.032	4.80E−02	0.145	4.63E−20	4.30E−07	
b-Cryptoxanthin	Albumin	0.037	5.03E−03	0.136	1.66E−17	1.52E−06	
g-Tocopherol	GGT	0.066	4.09E−12	0.129	7.50E−30	2.08E−05	
g-Tocopherol	Triglyceride	0.164	1.70E−17	0.309	4.81E−62	7.29E−08	
g-Tocopherol	Triglycerides	0.2	2.72E−51	0.319	5.69E−133	1.35E−10	
Mono-ethyl phthalate	Creatinine, urine	0.47	9.52E−136	0.354	1.68E−80	1.31E−05	
# cigarettes smoked per day now	Bilirubin, total	− 0.077	3.76E−04	− 0.218	4.96E−26	2.42E−06	
Avg # cigarettes/day during past 30 days	Bilirubin, total	− 0.072	4.99E−04	− 0.208	1.32E−23	4.13E−06	
Avg # cigarettes/day during past 30 days	Homocysteine	0.154	4.58E−16	0.038	1.57E−02	2.71E−06	
Current Cigarette Smoker?	Mean cell hemoglobin	0.062	8.80E−03	0.245	5.67E−15	2.93E−06	
Current Cigarette Smoker?	Segmented neutrophils number	0.175	1.37E−07	0.462	1.17E−51	2.36E−10	
Current Cigarette Smoker?	White blood cell count (SI)	0.255	9.13E−13	0.528	3.90E−59	1.68E−08	
Current or Past Cigarette Smoker?	1/Creatinine	0.083	6.93E−04	0.242	3.26E−25	2.72E−06	
Current or Past Cigarette Smoker?	Hematocrit	0.411	1.68E−52	0.232	7.51E−20	1.36E−06	
Current or Past Cigarette Smoker?	Hemoglobin	0.429	3.86E−54	0.236	1.42E−16	1.25E−06	
FTC Carbon Monoxide Content	1/Creatinine	0.031	4.14E−02	0.127	3.37E−18	6.00E−06	
How many years smoked this amount	Bilirubin, total	− 0.046	2.08E−02	− 0.185	6.37E−19	1.18E−06	
Does anyone smoke in home?	Hematocrit	0.351	4.12E−36	0.184	2.62E−13	8.24E−06	
Iron	Homocysteine	− 0.127	8.23E−54	− 0.064	3.17E−13	1.78E−07	
Weighted generalized linear model results, indicating the environmental exposure (predictor), the phenotype (outcome), and the sex-stratified β and P values, and the P values of the sex difference test. For more detailed information on the results the Supplementary Table 1 shows the complete set of results in each cohort, including the standard errors, and units of measurement

Qualitative sex differences

Twenty-five qualitative sex differences were identified, mostly involving nutrient exposure, food component recall, and heavy metals. Exposures to lead, lutein and zeaxanthin, retinyl palmitate, and vitamin A were associated with more than one phenotype, whereas the phenotypes albumin, C-reactive protein, and lymphocyte percent were associated with more than one exposure. Notably, various folate and iron measurements were negatively associated with albumin levels in females and positively associated in males (Fig. 3). Moreover, lead exposure was differentially associated with osmolality (βf=0.044,βm=-0.061,Pdiff=5.91E-08), bilirubin levels (βf=0.068,βm=-0.029,Pdiff=6.11E-07), and C-reactive protein levels (βf=-0.075,βm=0.049,Pdiff=8.48E-10) (Figs. 2, 3; Table 4).Table 4 Qualitative sex differences results.

Exposure name	Outcome name	βf	Pf	βm	Pm	Pdiff	
1,2,3,6,7,8-hxcdd	Chloride	− 0.133	4.55E−05	0.074	2.18E−02	6.46E−06	
Caffeine	C-reactive protein	− 0.058	9.03E−06	0.034	7.33E−03	4.45E−07	
Folate, DFE	Albumin	− 0.06	1.94E−07	0.029	4.70E−02	1.62E−06	
Food folate	Albumin	− 0.035	4.56E−04	0.044	6.20E−04	1.20E−06	
Potassium	Lymphocyte percent	− 0.044	6.40E−04	0.051	2.12E−03	6.13E−06	
Retinol	C-reactive protein	0.034	2.42E−02	− 0.066	9.09E−05	1.01E−05	
Thiamin (Vitamin B1)	Lymphocyte percent	− 0.064	3.94E−05	0.036	3.74E−02	1.82E−05	
Tuna eaten during past 30 days	Cholesterol, total	0.101	8.21E−04	− 0.166	4.31E−04	1.86E−06	
Lead	Bilirubin, total	0.067	2.12E−07	− 0.029	4.37E−02	6.11E−07	
Lead	C-reactive protein	− 0.075	1.03E−07	0.049	7.03E−04	8.48E−10	
Lead	Osmolality	0.044	1.16E−03	− 0.061	1.08E−05	5.90E−08	
Folate, RBC	Albumin	− 0.062	3.23E−07	0.05	3.84E−04	1.73E−09	
Lutein and zeaxanthin	Hematocrit	− 0.05	2.30E−05	0.031	2.20E−02	6.28E−06	
Lutein and zeaxanthin	Mean cell hemoglobin	0.028	2.71E−02	− 0.059	1.18E−07	2.21E−07	
Lutein and zeaxanthin	Red cell count SI	− 0.064	2.79E−05	0.058	4.04E−05	4.53E−09	
Retinyl Palmitate	Albumin	− 0.053	2.07E−05	0.05	1.22E−05	1.08E−09	
Retinyl Palmitate	Iron	0.037	3.43E−02	− 0.065	9.22E−07	3.24E−06	
Vitamin A	Glycohemoglobin	0.021	9.55E−03	− 0.037	3.06E−04	9.10E−06	
Vitamin A	Segmented neutrophils percent	− 0.041	2.19E−03	0.043	2.06E−03	1.40E−05	
PCB101	Bone alkaline phosphotase	− 0.111	8.63E−03	0.148	6.76E−05	4.16E−06	
Mono-n-methyl phthalate	Methylmalonic acid	− 0.032	3.23E−02	0.095	1.67E−05	1.93E−06	
Folic Acid	Albumin	− 0.026	1.28E−02	0.063	1.14E−07	1.75E−08	
Iron	Albumin	− 0.034	3.24E−04	0.043	1.13E−05	1.57E−08	
Blood m-/p-Xylene	Albumin	0.076	5.52E−06	− 0.05	1.41E−02	1.68E−06	
Blood o-Xylene	Albumin	0.071	8.43E−05	− 0.045	2.58E−02	1.81E−05	
Weighted generalized linear model results, indicating the environmental exposure (predictor), the phenotype (outcome), and the sex-stratified β and P values, and the P values of the sex difference test. For more detailed information on the results the Supplementary Table 1 shows the complete set of results in each cohort, including the standard errors, and units of measurement

Age stratification

When analyzing participants younger than 46, we identified 107 significant sex-specific associations: 43 were classified as qualitative, 27 as quantitative, and 37 as pure. Supplementary Table 2 contains the complete set of results for the young cohort. We identified 54 exposure-phenotype pairs replicated from the entire cohort, 34 of which showed the same type of sex difference classification, and 20 were classified differently. Notable examples of replicated findings with the same sex difference classification were the effects of cadmium on segmented neutrophils number (βf=0.11,βm=0.20,Pdiff=1.18E-05), C-reactive protein (βf=0.03,βm=0.12,Pdiff=4.10E-06), mean cell volume (βf=0.05,βm=0.14,Pdiff=5.45E-06), and albumin (βf=-0.002,βm=-0.08,Pdiff=2.47E-06), and the effect of smoking behavior (number of cigarettes smoked per day now and average number of cigarettes per day during past 30 days) and lead on homocysteine (βf=0.21,βm=0.03,Pdiff=3.06E-06;βf=0.20,βm=0.05,Pdiff=2.02E-05;βf=0.23,βm=0.12,Pdiff=7.22E-06). Albumin also showed consistent results, specifically in its association with lead, cadmium, and smoking behavior (current or past cigarette smoker)(βf=0.10,βm=-0.02,Pdiff=1.88E-07;βf=0.002,βm=-0.08,Pdiff=2.47E-06;βf=0.04,βm=-0.17,Pdiff=1.65E-07).

Even though 54 results were replicated in the younger cohort, 20 had different sex difference classifications. Notable examples were the associations of cotinine with red cell count (βf=0.053,βm=-0.04,Pdiff=2.28E-05) and albumin (βf=0.03,βm=-0.08,Pdiff=3.38E-10) which were classified as pure sex differences in the entire cohort, but as qualitative in the young cohort. Similarly, the association between smoking behavior (number of days smoked cigarettes during the past 30 days) and albumin (βf=0.04,βm=-0.10,Pdiff=1.65E-05) was classified as pure in the entire cohort, but qualitative in the young cohort.

On the other hand, we identified 53 novel exposure-phenotype pairs compared to the entire cohort. Of these, 12 associations were classified as quantitative, 22 as qualitative, and 19 as pure. Notable examples were the association between the FTC nicotine and tar content of the participants’ preferred cigarette brand and smoking behavior (how many years smoked this amount) with homocysteine (βf=0.18,βm=0.03,Pdiff=5.44E-06;βf=0.18,βm=0.03,Pdiff=2.54E-06), which showed a significant association only in females. Similarly, FTC carbon monoxide content of participants’ preferred cigarette brand, cotinine, and smoking family (does anyone smoke in home?) were positively associated with homocysteine across sexes but with greater effect sizes in females (βf=0.18,βm=0.04,Pdiff=1.42E-05;βf=0.16,βm=0.07,Pdiff=2.40E-05;βf=0.30,βm=0.08,Pdiff=8.17E-07).

Discussion

We conducted a PheEWAS to identify sex-specific associations between 58 phenotypes and 272 environmental exposures among more than 18,000 NHANES participants. Our study identified environmental exposures with sex-specific effects on several blood phenotypes. Our analysis’s high-throughput and agnostic nature provides avenues for further research into the biological, environmental, and mechanistic aspects of these sex differences. Furthermore, our study demonstrates the utility of using an explicit PheEWAS framework to target potential differences between sexes instead of evaluating sex differences as a complementary analysis.

We identified 119 exposure-phenotype associations that significantly differed between sexes, many of which were related to smoking behaviors, familial smoking, and cotinine levels. Tobacco smoke is a significant source of exposure to heavy metals, such as cadmium and lead34–37, and volatile compounds like benzene, toluene, ethylbenzene, and xylene38. Exposure to these substances can be interpreted as exposure to both active and passive tobacco smoke. Our findings indicate that the effects of heavy metals and volatile compounds on phenotypes, such as serum albumin, closely mirror the impacts of smoking behaviors, as reflected in cotinine levels and self-reported smoking information. For instance, cadmium exposure in males had a similar impact on albumin levels as direct measures of smoking, such as cotinine.

Previous studies have documented sex differences in the prevalence and health impacts of smoking39,40. Our research further supports these differences, showing that smoking-related exposures have distinct associations with various phenotypes in males and females. For example, smoking is associated with elevated homocysteine (Hcy) levels41, and smoking and elevated Hcy levels have been associated with increased risks of cardiovascular disease and COPD42–45. Although males generally have higher Hcy levels than females, smoking exacerbates Hcy levels more in females, potentially due to hormonal mechanisms, including estrogen’s role in methionine metabolism46. Therefore, females might be at greater risk of the harmful effects of elevated Hcy levels than males41,47, including cardiovascular disease3,48 and COPD4. Our study’s findings support sex differences in the effects of tobacco smoke exposure on Hcy. For example, the number of cigarettes currently smoked per day and the frequency of cigarette smoking during the past 30 days were positively associated with plasma Hcy levels only in females. Similarly, the average number of cigarettes smoked per day in the last 30 days had a greater effect on Hcy levels in females.

Smoking exposure has been linked to elevated levels of red blood cell populations, including hematocrit, hemoglobin, and red cell counts49,50. This increase is thought to result from carbon monoxide inhalation, which induces hypoxia and triggers a compensatory rise in hemoglobin, hematocrit, and red cell counts to maintain oxygen transport51–53. Our study found that smoking behavior, familial smoking, cotinine levels, and blood concentrations of cadmium, lead, ethylbenzene, and toluene were positively associated with hematocrit, hemoglobin, and red cell counts in females, consistent with prior research. However, to our knowledge, this study is the first to identify potential sex differences in the impact of tobacco smoke exposure on red blood cell populations.

On the other hand, tobacco smoke exposure significantly impacts inflammation and oxidative stress, primarily due to the presence of reactive oxygen species (ROS) that induce oxidative stress, leading to cellular apoptosis and senescence54. Typically, plasma antioxidants like albumin and bilirubin increase in response to oxidative stress as a protective mechanism55. However, studies have shown that smoking is associated with reduced levels of these antioxidants, even in a dose-dependent manner56–63. The excessive free radical production and extreme oxidative stress from smoking likely lead to the breakdown and depletion of these antioxidants, resulting in their lower levels64,65. Furthermore, smokers have also been found to exhibit elevated levels of inflammatory markers compared to non-smokers, including C-reactive protein (CRP), white blood cell count, and interleukin-662,66–68. Although the exact mechanisms are unclear, albumin levels tend to decrease in response to acute inflammation69. The combined oxidative stress and inflammatory responses to tobacco smoke exposure elevate the risk of chronic diseases such as cancer, COPD, and cardiovascular disease54,70.

In our study, smoking behavior, familial smoking, cotinine levels, and cadmium exposure were negatively associated with serum albumin levels in males but not in females. Additionally, smoking behavior, including the number of cigarettes smoked daily and the average number of cigarettes smoked daily over the past month, was negatively associated with bilirubin levels, with stronger effects observed in males. Moreover, smoking behavior and blood cadmium levels were positively associated with CRP levels and white blood cell counts, again with greater effect sizes in males. These findings align with previous research highlighting the role of tobacco smoke in promoting oxidative stress and inflammation. However, only a few studies have reported sex differences in the associations between tobacco smoke exposure and inflammatory markers or antioxidants56,71. Our findings suggest that males may experience a heightened inflammatory and oxidative stress response to tobacco smoke exposure compared to females. Notably, most observed differences were quantitative, indicating that these effects also occur in females but are less pronounced.

The findings from this study underscore the significant sex differences in the biological mechanisms underlying tobacco smoke exposure. These differences manifest in various processes, including inflammation, oxidative stress, hypoxia, and Hcy metabolism. Future research and prevention strategies must consider these sex-specific variations. For instance, methodological designs should go beyond merely controlling for sex and actively account for these differences. In this study, tobacco smoke exposure was assessed using multiple measures, including direct blood cotinine levels and self-reported smoking behavior. The consistent observation of sex differences across these varied metrics reinforces the robustness of our results. Going forward, researchers can leverage the extensive NHANES dataset to develop more detailed exposure metrics that distinguish between primary and secondary exposure and capture the duration and intensity of exposure.

We found that serum albumin levels seem more susceptible to environmental exposures in males than in females. Of the 22 pure or quantitative sex-specific associations involving serum albumin levels, 20 were exclusive to males or had a greater effect size in males. By contrast, associations involving serum Hcy levels tended to be specific to females. In all five pure or quantitative sex-specific associations involving serum Hcy, the effect was either exclusive to females or stronger in females than males. These results suggest that sex-specific mechanisms might underlie the metabolism of Hcy and serum albumin, implicating sex differences in methionine catabolism and inflammatory responses.

Our study identified three distinct types of sex differences in the associations between environmental exposures and phenotypes: pure, quantitative, and qualitative. Each type of difference has unique implications for understanding sex-specific environmental effects and developing targeted interventions. Pure sex differences occur when a significant association is observed in one sex but not in the other. This type of difference suggests that the biological mechanisms underlying the response to this exposure are either absent or not as pronounced in one sex. Pure sex differences often point to sex-specific biological pathways that may be influenced by factors such as sex hormones or genetic differences. Understanding these mechanisms can help identify potential targets for therapeutic interventions that are tailored to one sex. On the other hand, quantitative sex differences occur when an association is present in both sexes but with varying magnitudes of effect, suggesting that while both sexes are affected, the impact is greater in one sex. Quantitative differences highlight the need for sex-specific thresholds or guidelines in clinical practice. For example, our findings suggest that smoking exposure has a greater effect on Hcy levels in females, which can lead to more tailored interventions. Finally, qualitative sex differences occur when the direction of the association is opposite in males and females. Such opposite effects suggest fundamentally different biological responses to the same exposure in each sex, which have significant implications for tailored interventions and therapeutic guidance.

Of the 119 significant sex-specific associations identified in our PheEWAS, 94 were pure or quantitative, and only 25 were classified as qualitative. This favor toward quantitative and pure sex differences suggests that males and females tend to differ in the degree of associations but not in the direction of effects. Furthermore, quantitative and pure sex differences suggest associations with the same underlying mechanisms in males and females, whereas qualitative differences suggest associations with different mechanisms in males and females. Some differences, which are, in fact, quantitative, might have been misclassified as pure because of a lack of power to detect a significant association in one of the sexes. On the other hand, qualitative differences suggest completely different mechanisms or responses in each sex. The different underlying mechanisms of qualitative sex differences and pure or quantitative sex differences might also explain why qualitative differences were rarer and more heterogeneous than pure or quantitative differences. Pure and quantitative sex differences showed greater similarities than qualitative sex differences in the exposures and phenotypes involved, for example, smoking-related exposures. The observed bias can have implications for future research, particularly in suggesting experimental designs and analyses that privilege detecting differences in effect sizes between sexes. Nevertheless, the bias towards pure and quantitative sex differences might be due to the higher power to detect pure and quantitative sex differences at the expense of qualitative differences.

We replicated 54 exposure-phenotype associations using a younger cohort that retained roughly half of the participants. These replicated findings corroborated the results obtained in the entire cohort regarding the differential effects of tobacco smoke exposure between sexes. The non-replicated exposure-phenotype might be due to the lower sample sizes and, therefore, the reduced power to detect those differences. The unique exposure-phenotype associations that were present only in the young cohort might be due to age-specific differences, including differences in exposure rates in younger versus older participants. Nevertheless, our study did not rigorously test differences between the full and young cohorts, which is beyond the scope of our study. This requires a two-way interaction between sex and age, necessitating substantial sample sizes to estimate and detect significant differences. Therefore, these differences remain speculative, and comparisons between the full and young cohorts are only descriptive.

While this study provides valuable insights into sex-specific associations between environmental exposures and health outcomes, several limitations that could affect the interpretation of our findings must be acknowledged. One of the primary limitations of this study is the exclusion of participants with missing data on covariates, which was necessary to maintain the integrity and consistency of the dataset. While this approach helps ensure that all included participants have complete data for the variables under consideration, it may introduce selection bias, potentially affecting the study’s findings. In short, excluding participants with missing data can lead to a non-random sample that might not fully represent the broader population. This selection bias could manifest in several ways, including underestimating or overestimating associations, reduced generalizability, and losing statistical power. Future studies could employ several strategies to mitigate these potential biases, including multiple imputation techniques and sensitivity analyses.

PheEWAS involves testing thousands of associations simultaneously, which increases the risk of Type I errors (false positives). Although we applied the Bonferroni correction to account for multiple comparisons, a widely used method for controlling the risk of false positives, it is overly conservative, particularly in studies with many comparisons. This conservativeness can increase the risk of Type II errors (false negatives), where true associations may be overlooked because the adjusted p-value threshold is too stringent72. In the context of our study, this could result in missing potentially sex-specific associations between environmental exposures and phenotypes. Given the limitations of the Bonferroni correction, alternative methods such as False Discovery Rate (FDR) can be used in future studies that balance the risk of Type I and Type II errors more effectively but still may lead to some false discoveries73.

Another caveat to our study is that sex differences in the range of variation of the exposures might be observed as sex differences in the association between exposures and phenotypes. For example, if males have a minimal range of variation in a given exposure, there might not be enough variation to detect an association with a given phenotype. In contrast, if females have a greater range of variation in the exposure, it can reveal a positive association between the exposure and the phenotype. In such a case, a phenotype-exposure association would differ between males and females only because of the different degrees of exposure and not necessarily because of different underlying mechanisms. Future research might benefit from detecting significant sex differences in the distributions of the exposures and either removing those exposures or flagging them for potential manual inspection.

Although we identified several exposure-phenotype associations, because of the cross-sectional nature of the data, we can only determine correlations between exposures and phenotypes and not causality74. Furthermore, although we adjusted for several known confounders, residual confounding may still exist due to unmeasured or unknown factors. Finally, reverse causality remains an issue, as in EWAS more generally74. For example, several sex-different associations might be due to differential behaviors that men and women deploy in the face of similar events. In other words, individuals with pre-existing health conditions may alter their behaviors or environments in ways that influence their exposure levels. Such behavioral differences are not trivial in analyses of sex differences because studies have shown gender differences in health-related behaviors such as healthcare seeking and medication intake and adherence75–77. Therefore, our results cannot easily be ascribed simply to sex differences if we understand them solely based on biology.

Finally, the generalizability of our findings may be limited by the demographic characteristics of the NHANES sample, which, while representative of the U.S. population, may not capture the full diversity of global populations. Factors such as ethnicity and environmental exposures vary widely across different populations, and the associations observed in this study may not apply universally22.

Our study unveils novel insights into the sex-specific effects of environmental exposures on health outcomes, highlighting this research’s unique and significant contribution to the field. First, it underscores the need for a deeper understanding of biological mechanisms driving these differences. Identifying sex differences is the first step in discovering potential biological mechanisms that might include hormonal, genetic, and epigenetic influences. With a better understanding of these mechanisms, sex differences become a crude reflection of a more nuanced underlying biology. Second, our study underscores the pressing need for a precision medicine approach that not only relies on genetic background but also on low-throughput demographic information, such as sex, to assess differences in treatments, diagnoses, and prognoses, thereby providing a practical framework for healthcare professionals and policymakers. It also foregrounds the relevance of targeted public health policies. For example, public health policies might address the needs of diverse subgroups differently to mitigate risks produced by environmental exposures, such as tobacco smoke. Finally, future studies have the potential to address two crucial missing pieces in our study: longitudinal variation and multidimensional exposures and phenotypes, thereby paving the way for further significant advancements in this field. We have taken a more straightforward approach in our study by taking a snapshot of the effects of environmental exposures on endophenotypes, but differences in intensities, duration of exposures, and lifespan window can modulate these effects. Furthermore, multiple exposures can affect multiple phenotypes in non-linear and more complex ways.

Supplementary Information

Supplementary Information 1.

Supplementary Information 2.

Supplementary Information 3.

Supplementary Information 4.

Supplementary Information

The online version contains supplementary material available at 10.1038/s41598-024-72180-x.

Acknowledgements

We thank the editor and two anonymous reviewers for their comments and feedback that have improve this manuscript considerably.

Author contributions

T.G.Z., M.A.H., N.E.P., D.K., and L.S. contributed to the conception and design of the study. N.E.P. and T.G.Z. performed the cleaning of the databases and the analysis. T.G.Z. wrote the first draft of the manuscript. All authors contributed to the revision of the manuscript, and read and approved the submitted version.

Data availability

The datasets analysed during the current study are available in the DataDryad repository, under the following DOI: 10.1038/sdata.2016.96. All the code used during the current study are available in the Zenodo repository, under the following DOI: https://zenodo.org/doi/10.5281/zenodo.6483502.

Competing interests

The authors declare no competing interests.

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Mauvais-Jarvis F Sex and gender: Modifiers of health, disease, and medicine Lancet 2020 396 565 582 10.1016/S0140-6736(20)31561-0 32828189
Mauvais-Jarvis, F. et al. Sex and gender: Modifiers of health, disease, and medicine. Lancet 396, 565–582 (2020).32828189 10.1016/S0140-6736(20)31561-0
2. Klein SL Flanagan KL Sex differences in immune responses Nat. Rev. Immunol. 2016 16 626 638 10.1038/nri.2016.90 27546235
Klein, S. L. & Flanagan, K. L. Sex differences in immune responses. Nat. Rev. Immunol. 16, 626–638 (2016).27546235 10.1038/nri.2016.90
3. Regitz-Zagrosek V Kararigas G Mechanistic pathways of sex differences in cardiovascular disease Physiol. Rev. 2017 97 1 37 10.1152/physrev.00021.2015 27807199
Regitz-Zagrosek, V. & Kararigas, G. Mechanistic pathways of sex differences in cardiovascular disease. Physiol. Rev. 97, 1–37 (2017).27807199 10.1152/physrev.00021.2015
4. DeMeo DL Women manifest more severe COPD symptoms across the life course Int. J. Chron. Obstruct. Pulmon. Dis. 2018 13 3021 3029 10.2147/COPD.S160270 30319250
DeMeo, D. L. et al. Women manifest more severe COPD symptoms across the life course. Int. J. Chron. Obstruct. Pulmon. Dis. 13, 3021–3029 (2018).30319250 10.2147/COPD.S160270
5. Mallampalli MP Role of environment and sex differences in the development of autoimmune diseases: A roundtable meeting report J. Womens Health 2013 22 578 586 10.1089/jwh.2013.4462
Mallampalli, M. P. et al. Role of environment and sex differences in the development of autoimmune diseases: A roundtable meeting report. J. Womens Health 22, 578–586 (2013).10.1089/jwh.2013.4462
6. Raghavan D Jain R Increasing awareness of sex differences in airway diseases Respirology 2015 21 449 459 10.1111/resp.12702 26677803
Raghavan, D. & Jain, R. Increasing awareness of sex differences in airway diseases. Respirology 21, 449–459 (2015).26677803 10.1111/resp.12702
7. DeBolt C Harris D The impact of social determinants of health on gender disparities within respiratory medicine Clin. Chest Med. Gender Respir. Dis. 2021 42 407 415 10.1016/j.ccm.2021.04.003
DeBolt, C. & Harris, D. The impact of social determinants of health on gender disparities within respiratory medicine. Clin. Chest Med. Gender Respir. Dis. 42, 407–415 (2021).10.1016/j.ccm.2021.04.003
8. McHugh M.K Symanski E Pompeii L.A Delclos , G.L Prevalence of asthma by industry and occupation in the U.S. working population Am. J. Ind. Med. 2010 53 463 475 10.1002/ajim.20800 20187006
McHugh, M. .K., Symanski, E., Pompeii, L. .A. & Delclos, G. .L. Prevalence of asthma by industry and occupation in the U.S. working population. Am. J. Ind. Med. 53, 463–475 (2010).20187006 10.1002/ajim.20800
9. Centers for Disease Control and Prevention (CDC). Annual smoking-attributable mortality, years of potential life lost, and economic costs—United States, 1995–1999. MMWR Morb. Mortal. Wkly. Rep. 51, 300–303 (2002).
10. Gan WQ Man SP Postma DS Camp P Sin DD Female smokers beyond the perimenopausal period are at increased risk of chronic obstructive pulmonary disease: A systematic review and meta-analysis Respir. Res. 2006 7 52 10.1186/1465-9921-7-52 16571126
Gan, W. Q., Man, S. P., Postma, D. S., Camp, P. & Sin, D. D. Female smokers beyond the perimenopausal period are at increased risk of chronic obstructive pulmonary disease: A systematic review and meta-analysis. Respir. Res. 7, 52 (2006).16571126 10.1186/1465-9921-7-52
11. Pollard KM Gender differences in autoimmunity associated with exposure to environmental factors J. Autoimmun. 2012 38 J177 J186 10.1016/j.jaut.2011.11.007 22137891
Pollard, K. M. Gender differences in autoimmunity associated with exposure to environmental factors. J. Autoimmun. 38, J177–J186 (2012) (Special Issue: Gender, Sex Hormones, Pregnancy and Autoimmunity).22137891 10.1016/j.jaut.2011.11.007
12. Shannon G Gender equality in science, medicine, and global health: Where are we at and why does it matter? Lancet 2019 393 560 569 10.1016/S0140-6736(18)33135-0 30739691
Shannon, G. et al. Gender equality in science, medicine, and global health: Where are we at and why does it matter?. Lancet 393, 560–569 (2019).30739691 10.1016/S0140-6736(18)33135-0
13. Liu LY Schaub MA Sirota M Butte AJ Sex differences in disease risk from reported genome-wide association study findings Hum. Genet. 2012 131 353 364 10.1007/s00439-011-1081-y 21858542
Liu, L. Y., Schaub, M. A., Sirota, M. & Butte, A. J. Sex differences in disease risk from reported genome-wide association study findings. Hum. Genet. 131, 353–364 (2012).21858542 10.1007/s00439-011-1081-y
14. Bush WS Oetjens MT Crawford DC Unravelling the human genome-phenome relationship using phenome-wide association studies Nat. Rev. Genet. 2016 17 129 145 10.1038/nrg.2015.36 26875678
Bush, W. S., Oetjens, M. T. & Crawford, D. C. Unravelling the human genome-phenome relationship using phenome-wide association studies. Nat. Rev. Genet. 17, 129–145 (2016).26875678 10.1038/nrg.2015.36
15. Dumitrescu L Sex differences in the genetic predictors of Alzheimer’s pathology Brain 2019 142 2581 2589 10.1093/brain/awz206 31497858
Dumitrescu, L. et al. Sex differences in the genetic predictors of Alzheimer’s pathology. Brain 142, 2581–2589 (2019).31497858 10.1093/brain/awz206
16. Hu Y Rs34331204 regulates TSPAN13 expression and contributes to Alzheimer’s disease with sex differences Brain 2020 143 e95 10.1093/brain/awaa302 33175954
Hu, Y. et al. Rs34331204 regulates TSPAN13 expression and contributes to Alzheimer’s disease with sex differences. Brain 143, e95 (2020).33175954 10.1093/brain/awaa302
17. Ostrom QT Sex-specific glioma genome-wide association study identifies new risk locus at 3p21.31 in females, and finds sex-differences in risk at 8q24.21 Sci. Rep. 2018 8 7352 10.1038/s41598-018-24580-z 29743610
Ostrom, Q. T. et al. Sex-specific glioma genome-wide association study identifies new risk locus at 3p21.31 in females, and finds sex-differences in risk at 8q24.21. Sci. Rep. 8, 7352 (2018).29743610 10.1038/s41598-018-24580-z
18. Winkler TW The influence of age and sex on genetic associations with adult body size and shape: A large-scale genome-wide interaction study PLoS Genet. 2015 11 e1005378 10.1371/journal.pgen.1005378 26426971
Winkler, T. W. et al. The influence of age and sex on genetic associations with adult body size and shape: A large-scale genome-wide interaction study. PLoS Genet. 11, e1005378 (2015).26426971 10.1371/journal.pgen.1005378
19. Rask-Andersen M Karlsson T Ek WE Johansson Å Genome-wide association study of body fat distribution identifies adiposity loci and sex-specific genetic effects Nat. Commun. 2019 10 339 10.1038/s41467-018-08000-4 30664634
Rask-Andersen, M., Karlsson, T., Ek, W. E. & Johansson, Å. Genome-wide association study of body fat distribution identifies adiposity loci and sex-specific genetic effects. Nat. Commun. 10, 339 (2019).30664634 10.1038/s41467-018-08000-4
20. Choquet H A large multiethnic GWAS meta-analysis of cataract identifies new risk loci and sex-specific effects Nat. Commun. 2021 12 3595 10.1038/s41467-021-23873-8 34127677
Choquet, H. et al. A large multiethnic GWAS meta-analysis of cataract identifies new risk loci and sex-specific effects. Nat. Commun. 12, 3595 (2021).34127677 10.1038/s41467-021-23873-8
21. Zheng Y Design and methodology challenges of environment-wide association studies: A systematic review Environ. Res. 2020 183 109275 10.1016/j.envres.2020.109275 32105887
Zheng, Y. et al. Design and methodology challenges of environment-wide association studies: A systematic review. Environ. Res. 183, 109275 (2020).32105887 10.1016/j.envres.2020.109275
22. Patel CJ Bhattacharya J Butte AJ An environment-wide association study (EWAS) on type 2 diabetes mellitus PLoS ONE 2010 5 e10746 10.1371/journal.pone.0010746 20505766
Patel, C. J., Bhattacharya, J. & Butte, A. J. An environment-wide association study (EWAS) on type 2 diabetes mellitus. PLoS ONE 5, e10746 (2010).20505766 10.1371/journal.pone.0010746
23. Hall MA Environment-wide association study (EWAS) for type 2 diabetes in the Marshfield Personalized Medicine Research Project Biobank Biocomputing 2013 2014 200 211
Hall, M. A. et al. Environment-wide association study (EWAS) for type 2 diabetes in the Marshfield Personalized Medicine Research Project Biobank. Biocomputing 2014, 200–211 (2013).
24. Zhuang X Toward a panoramic perspective of the association between environmental factors and cardiovascular disease: An environment-wide association study from National Health and Nutrition Examination Survey 1999–2014 Environ. Int. 2018 118 146 153 10.1016/j.envint.2018.05.046 29879615
Zhuang, X. et al. Toward a panoramic perspective of the association between environmental factors and cardiovascular disease: An environment-wide association study from National Health and Nutrition Examination Survey 1999–2014. Environ. Int. 118, 146–153 (2018).29879615 10.1016/j.envint.2018.05.046
25. McGinnis DP Brownstein JS Patel CJ Environment-wide association study of blood pressure in the National Health and Nutrition Examination Survey (1999–2012) Sci. Rep. 2016 6 30373 10.1038/srep30373 27457472
McGinnis, D. P., Brownstein, J. S. & Patel, C. J. Environment-wide association study of blood pressure in the National Health and Nutrition Examination Survey (1999–2012). Sci. Rep. 6, 30373 (2016).27457472 10.1038/srep30373
26. Winkler TW Approaches to detect genetic effects that differ between two strata in genome-wide meta-analyses: Recommendations based on a systematic evaluation PLoS ONE 2017 12 e0181038 10.1371/journal.pone.0181038 28749953
Winkler, T. W. et al. Approaches to detect genetic effects that differ between two strata in genome-wide meta-analyses: Recommendations based on a systematic evaluation. PLoS ONE 12, e0181038 (2017).28749953 10.1371/journal.pone.0181038
27. Patel CJ A database of human exposomes and phenomes from the US National Health and Nutrition Examination Survey Sci. Data 2016 3 160096 10.1038/sdata.2016.96 27779619
Patel, C. J. et al. A database of human exposomes and phenomes from the US National Health and Nutrition Examination Survey. Sci. Data 3, 160096 (2016).27779619 10.1038/sdata.2016.96
28. Patel, C. J. et al. Data from: A database of human exposomes and phenomes from the US National Health and Nutrition Examination Survey 2016 (2020).
29. Hall MA Detection of pleiotropy through a phenome-wide association study (PheWAS) of epidemiologic data as part of the Environmental Architecture for Genes Linked to Environment (EAGLE) study PLoS Genet. 2014 10 e1004678 10.1371/journal.pgen.1004678 25474351
Hall, M. A. et al. Detection of pleiotropy through a phenome-wide association study (PheWAS) of epidemiologic data as part of the Environmental Architecture for Genes Linked to Environment (EAGLE) study. PLoS Genet. 10, e1004678 (2014).25474351 10.1371/journal.pgen.1004678
30. Lucas AM CLARITE facilitates the quality control and analysis process for EWAS of metabolic-related traits Front. Genet. 2019 10 1240 10.3389/fgene.2019.01240 31921293
Lucas, A. M. et al. CLARITE facilitates the quality control and analysis process for EWAS of metabolic-related traits. Front. Genet. 10, 1240 (2019).31921293 10.3389/fgene.2019.01240
31. West RM Best practice in statistics: The use of log transformation Ann. Clin. Biochem. 2022 59 162 165 10.1177/00045632211050531 34666549
West, R. M. Best practice in statistics: The use of log transformation. Ann. Clin. Biochem. 59, 162–165 (2022).34666549 10.1177/00045632211050531
32. Archer KJ Dumur CI Ramakrishnan V Graphical technique for identifying a monotonic variance stabilizing transformation for absolute gene intensity signals BMC Bioinform. 2004 5 60 10.1186/1471-2105-5-60
Archer, K. J., Dumur, C. I. & Ramakrishnan, V. Graphical technique for identifying a monotonic variance stabilizing transformation for absolute gene intensity signals. BMC Bioinform. 5, 60 (2004).10.1186/1471-2105-5-60
33. Willer CJ Li Y Abecasis GR METAL: Fast and efficient meta-analysis of genomewide association scans Bioinformatics 2010 26 2190 2191 10.1093/bioinformatics/btq340 20616382
Willer, C. J., Li, Y. & Abecasis, G. R. METAL: Fast and efficient meta-analysis of genomewide association scans. Bioinformatics 26, 2190–2191 (2010).20616382 10.1093/bioinformatics/btq340
34. Nordberg, G. F., Fowler, B. A. & Nordberg, M. Handbook on the Toxicology of Metals (Academic Press, 2014). ISBN: 978-0-12-397339-9.
35. Pappas RS Toxic elements in tobacco and in cigarette smoke: Inflammation and sensitization Metallomics 2011 3 1181 1198 10.1039/c1mt00066g 21799956
Pappas, R. S. Toxic elements in tobacco and in cigarette smoke: Inflammation and sensitization. Metallomics 3, 1181–1198 (2011).21799956 10.1039/c1mt00066g
36. Apostolou A Secondhand tobacco smoke: A source of lead exposure in US children and adolescents Am. J. Public Health 2012 102 714 722 10.2105/AJPH.2011.300161 21852639
Apostolou, A. et al. Secondhand tobacco smoke: A source of lead exposure in US children and adolescents. Am. J. Public Health 102, 714–722 (2012).21852639 10.2105/AJPH.2011.300161
37. Mannino DM Homa DM Matte T Hernandez-Avila M Active and passive smoking and blood lead levels in U.S. adults: Data from the Third National Health and Nutrition Examination Survey Nicotine Tob. Res. Off. J. Soc. Res. Nicotine Tob. 2005 7 557 564 10.1080/14622200500185264
Mannino, D. M., Homa, D. M., Matte, T. & Hernandez-Avila, M. Active and passive smoking and blood lead levels in U.S. adults: Data from the Third National Health and Nutrition Examination Survey. Nicotine Tob. Res. Off. J. Soc. Res. Nicotine Tob. 7, 557–564 (2005).10.1080/14622200500185264
38. Chambers DM Ocariz JM McGuirk MF Blount BC Impact of cigarette smoking on volatile organic compound (VOC) blood levels in the U.S. population: NHANES 2003–2004 Environ. Int. 2011 37 1321 1328 10.1016/j.envint.2011.05.016 21703688
Chambers, D. M., Ocariz, J. M., McGuirk, M. F. & Blount, B. C. Impact of cigarette smoking on volatile organic compound (VOC) blood levels in the U.S. population: NHANES 2003–2004. Environ. Int. 37, 1321–1328 (2011).21703688 10.1016/j.envint.2011.05.016
39. McHugh RK Votaw VR Sugarman DE Greenfield SF Sex and gender differences in substance use disorders Clin. Psychol. Rev. Gender Mental Health 2018 66 12 23
McHugh, R. K., Votaw, V. R., Sugarman, D. E. & Greenfield, S. F. Sex and gender differences in substance use disorders. Clin. Psychol. Rev. Gender Mental Health 66, 12–23 (2018).
40. Appelman Y van Rijn BB ten Haaf ME Boersma E Peters SAE Sex differences in cardiovascular risk factors and disease prevention Atherosclerosis 2015 241 211 218 10.1016/j.atherosclerosis.2015.01.027 25670232
Appelman, Y., van Rijn, B. B., ten Haaf, M. E., Boersma, E. & Peters, S. A. E. Sex differences in cardiovascular risk factors and disease prevention. Atherosclerosis 241, 211–218 (2015).25670232 10.1016/j.atherosclerosis.2015.01.027
41. Nygård O Refsum H Ueland PM Vollset SE Major lifestyle determinants of plasma total homocysteine distribution: The Hordaland Homocysteine Study Am. J. Clin. Nutr. 1998 67 263 270 10.1093/ajcn/67.2.263 9459374
Nygård, O., Refsum, H., Ueland, P. M. & Vollset, S. E. Major lifestyle determinants of plasma total homocysteine distribution: The Hordaland Homocysteine Study. Am. J. Clin. Nutr. 67, 263–270 (1998).9459374 10.1093/ajcn/67.2.263
42. O’Callaghan, P., Meleady, R., Fitzgerald, T., Graham, I. & European COMAC Group. Smoking and plasma homocysteine. Eur. Heart J. 23, 1580–1586 (2002).
43. Chaudhary D Sharma N Senapati S Serum homocysteine could be used as a predictive marker for chronic obstructive pulmonary disease: A meta-analysis Front. Public Health 2019 7 69 10.3389/fpubh.2019.00069 31019906
Chaudhary, D., Sharma, N. & Senapati, S. Serum homocysteine could be used as a predictive marker for chronic obstructive pulmonary disease: A meta-analysis. Front. Public Health 7, 69 (2019).31019906 10.3389/fpubh.2019.00069
44. Forey BA Thornton AJ Lee PN Systematic review with meta-analysis of the epidemiological evidence relating smoking to COPD, chronic bronchitis and emphysema BMC Pulm. Med. 2011 11 36 10.1186/1471-2466-11-36 21672193
Forey, B. A., Thornton, A. J. & Lee, P. N. Systematic review with meta-analysis of the epidemiological evidence relating smoking to COPD, chronic bronchitis and emphysema. BMC Pulm. Med. 11, 36 (2011).21672193 10.1186/1471-2466-11-36
45. Chrysant SG Chrysant GS The current status of homocysteine as a risk factor for cardiovascular disease: A mini review Expert Rev. Cardiovasc. Ther. 2018 16 559 565 10.1080/14779072.2018.1497974 29979619
Chrysant, S. G. & Chrysant, G. S. The current status of homocysteine as a risk factor for cardiovascular disease: A mini review. Expert Rev. Cardiovasc. Ther. 16, 559–565 (2018).29979619 10.1080/14779072.2018.1497974
46. Finkelstein JD Martin JJ Homocysteine Int. J. Biochem. Cell Biol. 2000 32 385 389 10.1016/S1357-2725(99)00138-7 10762063
Finkelstein, J. D. & Martin, J. J. Homocysteine. Int. J. Biochem. Cell Biol. 32, 385–389 (2000).10762063 10.1016/S1357-2725(99)00138-7
47. Stanisławska-Sachadyn A Evidence for sex differences in the determinants of homocysteine concentrations Mol. Genet. Metab. 2008 93 355 362 10.1016/j.ymgme.2007.11.004 18180190
Stanisławska-Sachadyn, A. et al. Evidence for sex differences in the determinants of homocysteine concentrations. Mol. Genet. Metab. 93, 355–362 (2008).18180190 10.1016/j.ymgme.2007.11.004
48. Huxley RR Woodward M Cigarette smoking as a risk factor for coronary heart disease in women compared with men: A systematic review and meta-analysis of prospective cohort studies Lancet (London, England) 2011 378 1297 1305 10.1016/S0140-6736(11)60781-2 21839503
Huxley, R. R. & Woodward, M. Cigarette smoking as a risk factor for coronary heart disease in women compared with men: A systematic review and meta-analysis of prospective cohort studies. Lancet (London, England) 378, 1297–1305 (2011).21839503 10.1016/S0140-6736(11)60781-2
49. Van Tiel ED Quitting smoking may restore hematological characteristics within five years Ann. Epidemiol. 2002 12 378 388 10.1016/S1047-2797(01)00282-4 12160596
Van Tiel, E. D. et al. Quitting smoking may restore hematological characteristics within five years. Ann. Epidemiol. 12, 378–388 (2002).12160596 10.1016/S1047-2797(01)00282-4
50. Malenica M Effect of cigarette smoking on haematological parameters in healthy population Med. Arch. 2017 71 132 136 10.5455/medarh.2017.71.132-136 28790546
Malenica, M. et al. Effect of cigarette smoking on haematological parameters in healthy population. Med. Arch. 71, 132–136 (2017).28790546 10.5455/medarh.2017.71.132-136
51. Helman N Rubenstein LS The effects of age, sex, and smoking on erythrocytes and leukocytes Am. J. Clin. Pathol. 1975 63 35 44 10.1093/ajcp/63.3.35 1111274
Helman, N. & Rubenstein, L. S. The effects of age, sex, and smoking on erythrocytes and leukocytes. Am. J. Clin. Pathol. 63, 35–44 (1975).1111274 10.1093/ajcp/63.3.35
52. Fricker M Chronic cigarette smoke exposure induces systemic hypoxia that drives intestinal dysfunction JCI Insight 2018 3 1 19 10.1172/jci.insight.94040
Fricker, M. et al. Chronic cigarette smoke exposure induces systemic hypoxia that drives intestinal dysfunction. JCI Insight 3, 1–19 (2018).10.1172/jci.insight.94040
53. Pedersen KM Smoking and increased white and red blood cells Arterioscler. Thromb. Vasc. Biol. 2019 39 965 977 10.1161/ATVBAHA.118.312338 30866659
Pedersen, K. M. et al. Smoking and increased white and red blood cells. Arterioscler. Thromb. Vasc. Biol. 39, 965–977 (2019).30866659 10.1161/ATVBAHA.118.312338
54. Caliri AW Tommasi S Besaratinia A Relationships among smoking, oxidative stress, inflammation, macromolecular damage, and cancer Mutat. Res. Rev. Mutat. Res. 2021 787 108365 10.1016/j.mrrev.2021.108365 34083039
Caliri, A. W., Tommasi, S. & Besaratinia, A. Relationships among smoking, oxidative stress, inflammation, macromolecular damage, and cancer. Mutat. Res. Rev. Mutat. Res. 787, 108365 (2021).34083039 10.1016/j.mrrev.2021.108365
55. Amir Aslani B Ghobadi S Studies on oxidants and antioxidants with a brief glance at their relevance to the immune system Life Sci. 2016 146 163 173 10.1016/j.lfs.2016.01.014 26792059
Amir Aslani, B. & Ghobadi, S. Studies on oxidants and antioxidants with a brief glance at their relevance to the immune system. Life Sci. 146, 163–173 (2016).26792059 10.1016/j.lfs.2016.01.014
56. Hoydonck PGV Temme EH Schouten EG Serum bilirubin concentration in a Belgian population: The association with smoking status and type of cigarettes Int. J. Epidemiol. 2001 30 1465 1472 10.1093/ije/30.6.1465 11821365
Hoydonck, P. G. V., Temme, E. H. & Schouten, E. G. Serum bilirubin concentration in a Belgian population: The association with smoking status and type of cigarettes. Int. J. Epidemiol. 30, 1465–1472 (2001).11821365 10.1093/ije/30.6.1465
57. Schwertner HA Association of smoking and low serum bilirubin antioxidant concentrations Atherosclerosis 1998 136 383 387 10.1016/S0021-9150(97)00232-3 9543110
Schwertner, H. A. Association of smoking and low serum bilirubin antioxidant concentrations. Atherosclerosis 136, 383–387 (1998).9543110 10.1016/S0021-9150(97)00232-3
58. Jo J Kimm H Yun JE Lee KJ Jee SH Cigarette smoking and serum bilirubin subtypes in healthy Korean men: The Korea Medical Institute Study J. Prev. Med. Public Health 2012 45 105 112 10.3961/jpmph.2012.45.2.105 22509451
Jo, J., Kimm, H., Yun, J. E., Lee, K. J. & Jee, S. H. Cigarette smoking and serum bilirubin subtypes in healthy Korean men: The Korea Medical Institute Study. J. Prev. Med. Public Health 45, 105–112 (2012).22509451 10.3961/jpmph.2012.45.2.105
59. Kodal JB Smoking reduces plasma bilirubin: Observational and genetic analyses in the Copenhagen General Population Study Nicotine Tob. Res. 2020 22 104 110 30202916
Kodal, J. B. et al. Smoking reduces plasma bilirubin: Observational and genetic analyses in the Copenhagen General Population Study. Nicotine Tob. Res. 22, 104–110 (2020).30202916
60. Nelson JJ Serum albumin level as a predictor of incident coronary heart disease: The Atherosclerosis Risk in Communities (ARIC) Study Am. J. Epidemiol. 2000 151 468 477 10.1093/oxfordjournals.aje.a010232 10707915
Nelson, J. J. et al. Serum albumin level as a predictor of incident coronary heart disease: The Atherosclerosis Risk in Communities (ARIC) Study. Am. J. Epidemiol. 151, 468–477 (2000).10707915 10.1093/oxfordjournals.aje.a010232
61. Shaper AG Wannamethee SG Whincup PH Serum albumin and risk of stroke, coronary heart disease, and mortality: The role of cigarette smoking J. Clin. Epidemiol. 2004 57 195 202 10.1016/j.jclinepi.2003.07.001 15125630
Shaper, A. G., Wannamethee, S. G. & Whincup, P. H. Serum albumin and risk of stroke, coronary heart disease, and mortality: The role of cigarette smoking. J. Clin. Epidemiol. 57, 195–202 (2004).15125630 10.1016/j.jclinepi.2003.07.001
62. Bakhru A Erlinger TP Smoking cessation and cardiovascular disease risk factors: Results from the Third National Health and Nutrition Examination Survey PLoS Med. 2005 2 e160 10.1371/journal.pmed.0020160 15974805
Bakhru, A. & Erlinger, T. P. Smoking cessation and cardiovascular disease risk factors: Results from the Third National Health and Nutrition Examination Survey. PLoS Med. 2, e160 (2005).15974805 10.1371/journal.pmed.0020160
63. Jang ES Effects of coffee, smoking, and alcohol on liver function tests: A comprehensive cross-sectional study BMC Gastroenterol. 2012 12 1 12 10.1186/1471-230X-12-145 22214224
Jang, E. S. et al. Effects of coffee, smoking, and alcohol on liver function tests: A comprehensive cross-sectional study. BMC Gastroenterol. 12, 1–12 (2012).22214224 10.1186/1471-230X-12-145
64. Zhou JF Effects of cigarette smoking and smoking cessation on plasma constituents and enzyme activities related to oxidative stress Biomed. Environ. Sci. 2000 13 44 55 10853840
Zhou, J. F. et al. Effects of cigarette smoking and smoking cessation on plasma constituents and enzyme activities related to oxidative stress. Biomed. Environ. Sci. 13, 44–55 (2000).10853840
65. Yildiz L Kayaoğlu N Aksoy H The changes of superoxide dismutase, catalase and glutathione peroxidase activities in erythrocytes of active and passive smokers Clin. Chem. Lab. Med. 2002 40 612 615 10.1515/CCLM.2002.106 12211658
Yildiz, L., Kayaoğlu, N. & Aksoy, H. The changes of superoxide dismutase, catalase and glutathione peroxidase activities in erythrocytes of active and passive smokers. Clin. Chem. Lab. Med. 40, 612–615 (2002).12211658 10.1515/CCLM.2002.106
66. Sopori M Effects of cigarette smoke on the immune system Nat. Rev. Immunol. 2002 2 372 377 10.1038/nri803 12033743
Sopori, M. Effects of cigarette smoke on the immune system. Nat. Rev. Immunol. 2, 372–377 (2002).12033743 10.1038/nri803
67. Stämpfli MR Anderson GP How cigarette smoke skews immune responses to promote infection, lung disease and cancer Nat. Rev. Immunol. 2009 9 377 384 10.1038/nri2530 19330016
Stämpfli, M. R. & Anderson, G. P. How cigarette smoke skews immune responses to promote infection, lung disease and cancer. Nat. Rev. Immunol. 9, 377–384 (2009).19330016 10.1038/nri2530
68. Wannamethee SG Associations between cigarette smoking, pipe/cigar smoking, and smoking cessation, and haemostatic and inflammatory markers for cardiovascular disease Eur. Heart J. 2005 26 1765 1773 10.1093/eurheartj/ehi183 15817606
Wannamethee, S. G. et al. Associations between cigarette smoking, pipe/cigar smoking, and smoking cessation, and haemostatic and inflammatory markers for cardiovascular disease. Eur. Heart J. 26, 1765–1773 (2005).15817606 10.1093/eurheartj/ehi183
69. Chien S-C Chen C-Y Lin C-F Yeh H-I Critical appraisal of the role of serum albumin in cardiovascular disease Biomark. Res. 2017 5 31 10.1186/s40364-017-0111-x 29152305
Chien, S.-C., Chen, C.-Y., Lin, C.-F. & Yeh, H.-I. Critical appraisal of the role of serum albumin in cardiovascular disease. Biomark. Res. 5, 31 (2017).29152305 10.1186/s40364-017-0111-x
70. Kotlyarov S The role of smoking in the mechanisms of development of chronic obstructive pulmonary disease and atherosclerosis Int. J. Mol. Sci. 2023 24 8725 10.3390/ijms24108725 37240069
Kotlyarov, S. The role of smoking in the mechanisms of development of chronic obstructive pulmonary disease and atherosclerosis. Int. J. Mol. Sci. 24, 8725 (2023).37240069 10.3390/ijms24108725
71. Fröhlich M Independent association of various smoking characteristics with markers of systemic inflammation in men: Results from a Representative Sample of the General Population (MONICA Augsburg Survey 1994/95) Eur. Heart J. 2003 24 1365 1372 10.1016/S0195-668X(03)00260-4 12871694
Fröhlich, M. et al. Independent association of various smoking characteristics with markers of systemic inflammation in men: Results from a Representative Sample of the General Population (MONICA Augsburg Survey 1994/95). Eur. Heart J. 24, 1365–1372 (2003).12871694 10.1016/S0195-668X(03)00260-4
72. Armstrong RA When to use the Bonferroni correction Ophthalmic Physiol. Opt. 2014 34 502 508 10.1111/opo.12131 24697967
Armstrong, R. A. When to use the Bonferroni correction. Ophthalmic Physiol. Opt. 34, 502–508 (2014).24697967 10.1111/opo.12131
73. Benjamini Y Hochberg Y Controlling the false discovery rate: A practical and powerful approach to multiple testing J. R. Stat. Soc. Ser. B (Methodol.) 1995 57 289 300 10.1111/j.2517-6161.1995.tb02031.x
Benjamini, Y. & Hochberg, Y. Controlling the false discovery rate: A practical and powerful approach to multiple testing. J. R. Stat. Soc. Ser. B (Methodol.) 57, 289–300 (1995).10.1111/j.2517-6161.1995.tb02031.x
74. Patel CJ Ioannidis JPA Studying the elusive environment in large scale JAMA 2014 311 2173 2174 10.1001/jama.2014.4129 24893084
Patel, C. J. & Ioannidis, J. P. A. Studying the elusive environment in large scale. JAMA 311, 2173–2174 (2014).24893084 10.1001/jama.2014.4129
75. Thompson AE The influence of gender and other patient characteristics on health care-seeking behaviour: A QUALICOPC study BMC Fam. Pract. 2016 17 1 7 10.1186/s12875-016-0440-0 26739311
Thompson, A. E. et al. The influence of gender and other patient characteristics on health care-seeking behaviour: A QUALICOPC study. BMC Fam. Pract. 17, 1–7 (2016).26739311 10.1186/s12875-016-0440-0
76. Chen S-L Lee W-L Liang T Liao I-C Factors associated with gender differences in medication adherence: A longitudinal study J. Adv. Nurs. 2014 70 2031 2040 10.1111/jan.12361 24506542
Chen, S.-L., Lee, W.-L., Liang, T. & Liao, I.-C. Factors associated with gender differences in medication adherence: A longitudinal study. J. Adv. Nurs. 70, 2031–2040 (2014).24506542 10.1111/jan.12361
77. Manteuffel M Influence of patient sex and gender on medication use, adherence, and prescribing alignment with guidelines J. Womens Health 2014 23 112 119 10.1089/jwh.2012.3972
Manteuffel, M. et al. Influence of patient sex and gender on medication use, adherence, and prescribing alignment with guidelines. J. Womens Health 23, 112–119 (2014).10.1089/jwh.2012.3972
