A quantile integral linear model to quantify genetic effects on phenotypic variability.

Detecting genetic variants associated with the variance of complex traits, that is, variance quantitative trait loci (vQTLs), can provide crucial insights into the interplay between genes and environments and how they jointly shape human phenotypes in the population. We propose a quantile integral linear model (QUAIL) to estimate genetic effects on trait variability. Through extensive simulations and analyses of real data, we demonstrate that QUAIL provides computationally efficient and statistically powerful vQTL mapping that is robust to non-Gaussian phenotypes and confounding effects on phenotypic variability. Applied to UK Biobank (n = 375,791), QUAIL identified 11 vQTLs for body mass index (BMI) that have not been previously reported. Top vQTL findings showed substantial enrichment for interactions with physical activities and sedentary behavior. Furthermore, variance polygenic scores (vPGSs) based on QUAIL effect estimates showed superior predictive performance on both population-level and within-individual BMI variability compared to existing approaches. Overall, QUAIL is a unified framework to quantify genetic effects on the phenotypic variability at both single-variant and vPGS levels. It addresses critical limitations in existing approaches and may have broad applications in future gene-environment interaction studies.

Feature-specific quantile normalization and feature-specific mean-variance normalization deliver robust bi-directional classification and feature selection performance between microarray and RNAseq data.

BACKGROUND: Cross-platform normalization seeks to minimize technological bias between microarray and RNAseq whole-transcriptome data. Incorporating multiple gene expression platforms permits external validation of experimental findings, and augments training sets for machine learning models. Here, we compare the performance of Feature Specific Quantile Normalization (FSQN) to a previously used but unvalidated and uncharacterized method we label as Feature Specific Mean Variance Normalization (FSMVN). We evaluate the performance of these methods for bidirectional normalization in the context of nested feature selection. RESULTS: FSQN and FSMVN provided clinically equivalent bidirectional model performance with and without feature selection for colon CMS and breast PAM50 classification. Using principal component analysis, we determine that these methods eliminate batch effects related to technological platforms. Without feature selection, no statistical difference was identified between the performance of FSQN and FSMVN of cross-platform data compared to within-platform distributions. Under optimal feature selection conditions, balanced accuracy was FSQN and FSMVN were statistically equivalent to the within-platform distribution performance in multivariable linear regression analysis. FSQN and FSMVN also provided similar performance to within-platform distributions as the number of selected genes used to create models decreases. CONCLUSIONS: In the context of generating supervised machine learning classifiers for molecular subtypes, FSQN and FSMVN are equally effective. Under optimal modeling conditions, FSQN and FSMVN provide equivalent model accuracy performance on cross-platform normalization data compared to within-platform data. Using cross-platform data should still be approached with caution as subtle performance differences may exist depending on the classification problem, training, and testing distributions.

Quantile regressions as a tool to evaluate how an exposure shifts and reshapes the outcome distribution: A primer for epidemiologists.

Quantifying how an exposure affects the entire outcome distribution is often important, e.g., for outcomes such as blood pressure which have non-linear effects on long-term morbidity and mortality. Quantile regressions offer a powerful way of estimating an exposure's relationship with the outcome distribution but remain underused in epidemiology. We introduce quantile regressions with a focus on distinguishing estimators for quantiles of the conditional and unconditional outcome distributions. We also present an empirical example in which we fit mean and quantile regressions to investigate educational attainment's association with later-life systolic blood pressure (SBP). We use data on 8,875 US-born respondents aged 50+ years from the US Health and Retirement Study. More education was negatively associated with mean SBP. Conditional and unconditional quantile regressions both suggested a negative association between education and SBP at all levels of SBP, but the absolute magnitudes of these associations were higher at higher SBP quantiles relative to lower quantiles. In addition to showing that educational attainment shifted the SBP distribution left-wards, quantile regression results revealed that education may have reshaped the SBP distribution through larger protective associations in the right tail, thus benefiting those at highest risk of cardiovascular diseases.

Impact of Vietnam-era G.I. Bill eligibility on later-life blood pressure distribution: evidence from the Vietnam draft lottery natural experiment.

We quantified the impact of Vietnam-era G.I Bill eligibility, which subsidized college education for eligible Veterans, on the later-life blood pressure distribution by exploiting the Vietnam draft lottery natural experiment. We restricted Health and Retirement Study data (2006-2018) to men born between 1947-1953 (N=1,970). We estimated intention-to-treat effects at the mean and 1st-99th quantiles of blood pressure using linear and quantile regressions. Our outcomes were systolic blood pressure (SBP), diastolic blood pressure (DBP), hypertension, and self-reported stroke. We proxied G.I. Bill eligibility using lottery-defined draft eligibility. We also conducted analyses stratified by childhood socioeconomic status (cSES) defined based on a previously validated measure. Draft eligibility reduced mean blood pressure outcomes (e.g., effect on SBP: -1.33 [95% confidence interval (CI) -2.85, 0.19]). Draft eligibility also had larger protective effects at higher quantiles of the SBP and DBP distributions relative to lower quantiles (effects on SBP at the 10th and 90th quantiles: -0.33mmHg [95% CI -2.35,1.68]; -3.00mmHg [95% CI -5.68,-0.32]). Draft eligibility had protective effects on blood pressure among low and medium cSES men but opposite effects among high cSES men. G.I. Bill eligibility reshaped the blood pressure distribution to one of lower morbidity risk, particularly among low and medium cSES men.

Distributional data analysis via quantile functions and its application to modeling digital biomarkers of gait in Alzheimer's Disease.

With the advent of continuous health monitoring with wearable devices, users now generate their unique streams of continuous data such as minute-level step counts or heartbeats. Summarizing these streams via scalar summaries often ignores the distributional nature of wearable data and almost unavoidably leads to the loss of critical information. We propose to capture the distributional nature of wearable data via user-specific quantile functions (QF) and use these QFs as predictors in scalar-on-quantile-function-regression (SOQFR). As an alternative approach, we also propose to represent QFs via user-specific L-moments, robust rank-based analogs of traditional moments, and use L-moments as predictors in SOQFR (SOQFR-L). These two approaches provide two mutually consistent interpretations: in terms of quantile levels by SOQFR and in terms of L-moments by SOQFR-L. We also demonstrate how to deal with multi-modal distributional data via Joint and Individual Variation Explained using L-moments. The proposed methods are illustrated in a study of association of digital gait biomarkers with cognitive function in Alzheimers disease. Our analysis shows that the proposed methods demonstrate higher predictive performance and attain much stronger associations with clinical cognitive scales compared to simple distributional summaries.

Among-species variation in six decades of changing migration timings explained through ecology, life-history and local migratory abundance.

Species exploiting seasonal environments must alter timings of key life-history events in response to large-scale climatic changes in order to maintain trophic synchrony with required resources. Yet, substantial among-species variation in long-term phenological changes has been observed. Advancing from simply describing such variation towards predicting future phenological responses requires studies that rigorously quantify and explain variation in the direction and magnitude of changing timings across diverse species in relation to key ecological and life-history variables. Accordingly, we fitted multi-quantile regressions to 59 years of multi-species data on spring and autumn bird migration timings through northern Scotland. We demonstrate substantial variation in changes in timings among 72 species, and tested whether such variation can be explained by species ecology, life-history and changes in local abundance. Consistent with predictions, species that advanced their migration timing in one or both seasons had more seasonally restricted diet types, fewer suitable breeding habitat types, shorter generation lengths and capability to produce multiple offspring broods per year. In contrast, species with less seasonally restricted diet types and that produce single annual offspring broods, showed no change. Meanwhile, contrary to prediction, long-distance and short-distance migrants advanced migration timings similarly. Changes in migration timing also varied with changes in local migratory abundance, such that species with increasing seasonal abundance apparently altered their migration timing, whilst species with decreasing abundance did not. Such patterns broadly concur with expectation given adaptive changes in migration timing. However, we demonstrate that similar patterns can be generated by numerical sampling given changing local abundances. Any apparent phenology-abundance relationships should, therefore, be carefully validated and interpreted. Overall, our results show that migrant bird species with differing ecologies and life-histories showed systematically differing phenological changes over six decades contextualised by large-scale environmental changes, potentially facilitating future predictions and altering temporal dynamics of seasonal species co-occurrences.

Balancing the benefits and risks of China's national salt iodization policy over 30 years using disability-adjusted life years (DALYs): a systematic review and meta-analysis.

Both insufficient and excessive iodine intake can lead to thyroid-related disorders. Although China has made progress in eliminating iodine deficiency over the past few decades, the incidence of thyroid cancer is increasing. Currently, there is a lack of relevant research on the tradeoff between the benefits and risks of salt iodization in China. In this study, we developed a method that combines the total probability algorithm and disease burden to evaluate the appropriate amount of salt iodization. Following the principle of minimizing the comprehensive disease burden and using the metabolic model of human iodine nutrition. Based on the average national iodine level in water, the optimal iodine content in Chinese salt is determined to be 17 mg/kg. However, iodine content in water is not evenly distributed in China. Approximately 3.23% of administrative villages have water iodine concentrations exceeding 80 ug/L, eliminating the need for iodine fortification in salt. Approximately 83.51% of administrative villages need to continue implementing the salt iodization policy, with the optimal iodine content in salt ranging from 15 to 18 mg/kg. In 13.16% of administrative villages, the iodine content in salt is determined based on the local water iodine concentration, ranging from 0 to 15 mg/kg. Our study cracks open a window of insight suggesting that the optimal iodine content for salt is lower than the existing benchmark dictated by the prevailing policy in China. Hence, there is an urgent need to refine and advance the iodine supplementation strategy in salt to pave the way for precision medicine and health-centric iodine supplementation strategies.

Testing conditional quantile independence with functional covariate.

We propose a new non-parametric conditional independence test for a scalar response and a functional covariate over a continuum of quantile levels. We build a Cramer-von Mises type test statistic based on an empirical process indexed by random projections of the functional covariate, effectively avoiding the "curse of dimensionality" under the projected hypothesis, which is almost surely equivalent to the null hypothesis. The asymptotic null distribution of the proposed test statistic is obtained under some mild assumptions. The asymptotic global and local power properties of our test statistic are then investigated. We specifically demonstrate that the statistic is able to detect a broad class of local alternatives converging to the null at the parametric rate. Additionally, we recommend a simple multiplier bootstrap approach for estimating the critical values. The finite-sample performance of our statistic is examined through several Monte Carlo simulation experiments. Finally, an analysis of an EEG data set is used to show the utility and versatility of our proposed test statistic.

A scalar-on-quantile-function approach for estimating short-term health effects of environmental exposures.

Environmental epidemiologic studies routinely utilize aggregate health outcomes to estimate effects of short-term (eg, daily) exposures that are available at increasingly fine spatial resolutions. However, areal averages are typically used to derive population-level exposure, which cannot capture the spatial variation and individual heterogeneity in exposures that may occur within the spatial and temporal unit of interest (eg, within a day or ZIP code). We propose a general modeling approach to incorporate within-unit exposure heterogeneity in health analyses via exposure quantile functions. Furthermore, by viewing the exposure quantile function as a functional covariate, our approach provides additional flexibility in characterizing associations at different quantile levels. We apply the proposed approach to an analysis of air pollution and emergency department (ED) visits in Atlanta over 4 years. The analysis utilizes daily ZIP code-level distributions of personal exposures to 4 traffic-related ambient air pollutants simulated from the Stochastic Human Exposure and Dose Simulator. Our analyses find that effects of carbon monoxide on respiratory and cardiovascular disease ED visits are more pronounced with changes in lower quantiles of the population's exposure. Software for implement is provided in the R package nbRegQF.

High-dimensional sparse vine copula regression with application to genomic prediction.

High-dimensional data sets are often available in genome-enabled predictions. Such data sets include nonlinear relationships with complex dependence structures. For such situations, vine copula-based (quantile) regression is an important tool. However, the current vine copula-based regression approaches do not scale up to high and ultra-high dimensions. To perform high-dimensional sparse vine copula-based regression, we propose 2 methods. First, we show their superiority regarding computational complexity over the existing methods. Second, we define relevant, irrelevant, and redundant explanatory variables for quantile regression. Then, we show our method's power in selecting relevant variables and prediction accuracy in high-dimensional sparse data sets via simulation studies. Next, we apply the proposed methods to the high-dimensional real data, aiming at the genomic prediction of maize traits. Some data processing and feature extraction steps for the real data are further discussed. Finally, we show the advantage of our methods over linear models and quantile regression forests in simulation studies and real data applications.

Doubly robust estimation and sensitivity analysis for marginal structural quantile models.

The marginal structure quantile model (MSQM) provides a unique lens to understand the causal effect of a time-varying treatment on the full distribution of potential outcomes. Under the semiparametric framework, we derive the efficiency influence function for the MSQM, from which a new doubly robust estimator is proposed for point estimation and inference. We show that the doubly robust estimator is consistent if either of the models associated with treatment assignment or the potential outcome distributions is correctly specified, and is semiparametric efficient if both models are correct. To implement the doubly robust MSQM estimator, we propose to solve a smoothed estimating equation to facilitate efficient computation of the point and variance estimates. In addition, we develop a confounding function approach to investigate the sensitivity of several MSQM estimators when the sequential ignorability assumption is violated. Extensive simulations are conducted to examine the finite-sample performance characteristics of the proposed methods. We apply the proposed methods to the Yale New Haven Health System Electronic Health Record data to study the effect of antihypertensive medications to patients with severe hypertension and assess the robustness of the findings to unmeasured baseline and time-varying confounding.

Penalized weighted smoothed quantile regression for high-dimensional longitudinal data.

Quantile regression, known as a robust alternative to linear regression, has been widely used in statistical modeling and inference. In this paper, we propose a penalized weighted convolution-type smoothed method for variable selection and robust parameter estimation of the quantile regression with high dimensional longitudinal data. The proposed method utilizes a twice-differentiable and smoothed loss function instead of the check function in quantile regression without penalty, and can select the important covariates consistently using the efficient gradient-based iterative algorithms when the dimension of covariates is larger than the sample size. Moreover, the proposed method can circumvent the influence of outliers in the response variable and/or the covariates. To incorporate the correlation within each subject and enhance the accuracy of the parameter estimation, a two-step weighted estimation method is also established. Furthermore, we prove the oracle properties of the proposed method under some regularity conditions. Finally, the performance of the proposed method is demonstrated by simulation studies and two real examples.

The spike-and-slab quantile LASSO for robust variable selection in cancer genomics studies.

Data irregularity in cancer genomics studies has been widely observed in the form of outliers and heavy-tailed distributions in the complex traits. In the past decade, robust variable selection methods have emerged as powerful alternatives to the nonrobust ones to identify important genes associated with heterogeneous disease traits and build superior predictive models. In this study, to keep the remarkable features of the quantile LASSO and fully Bayesian regularized quantile regression while overcoming their disadvantage in the analysis of high-dimensional genomics data, we propose the spike-and-slab quantile LASSO through a fully Bayesian spike-and-slab formulation under the robust likelihood by adopting the asymmetric Laplace distribution (ALD). The proposed robust method has inherited the prominent properties of selective shrinkage and self-adaptivity to the sparsity pattern from the spike-and-slab LASSO (Roc̆ková and George, J Am Stat Associat, 2018, 113(521): 431-444). Furthermore, the spike-and-slab quantile LASSO has a computational advantage to locate the posterior modes via soft-thresholding rule guided Expectation-Maximization (EM) steps in the coordinate descent framework, a phenomenon rarely observed for robust regularization with nondifferentiable loss functions. We have conducted comprehensive simulation studies with a variety of heavy-tailed errors in both homogeneous and heterogeneous model settings to demonstrate the superiority of the spike-and-slab quantile LASSO over its competing methods. The advantage of the proposed method has been further demonstrated in case studies of the lung adenocarcinomas (LUAD) and skin cutaneous melanoma (SKCM) data from The Cancer Genome Atlas (TCGA).

Principal stratification for quantile causal effects under partial compliance.

Within the principal stratification framework in causal inference, the majority of the literature has focused on binary compliance with an intervention and modelling means. Yet in some research areas, compliance is partial, and research questions-and hence analyses-are concerned with causal effects on (possibly high) quantiles rather than on shifts in average outcomes. Modelling partial compliance is challenging because it can suffer from lack of identifiability. We develop an approach to estimate quantile causal effects within a principal stratification framework, where principal strata are defined by the bivariate vector of (partial) compliance to the two levels of a binary intervention. We propose a conditional copula approach to impute the missing potential compliance and estimate the principal quantile treatment effect surface at high quantiles, allowing the copula association parameter to vary with the covariates. A bootstrap procedure is used to estimate the parameter to account for inflation due to imputation of missing compliance. Moreover, we describe precise assumptions on which the proposed approach is based, and investigate the finite sample behavior of our method by a simulation study. The proposed approach is used to study the 90th principal quantile treatment effect of executive stay-at-home orders on mitigating the risk of COVID-19 transmission in the United States.

Multi-Factor Coral Disease Risk: A new product for early warning and management.

Ecological forecasts are becoming increasingly valuable tools for conservation and management. However, there are few examples of near-real-time forecasting systems that account for the wide range of ecological complexities. We developed a new coral disease ecological forecasting system that explores a suite of ecological relationships and their uncertainty and investigates how forecast skill changes with shorter lead times. The Multi-Factor Coral Disease Risk product introduced here uses a combination of ecological and marine environmental conditions to predict the risk of white syndromes and growth anomalies across reefs in the central and western Pacific and along the east coast of Australia and is available through the US National Oceanic and Atmospheric Administration Coral Reef Watch program. This product produces weekly forecasts for a moving window of 6 months at a resolution of ~5 km based on quantile regression forests. The forecasts show superior skill at predicting disease risk on withheld survey data from 2012 to 2020 compared with predecessor forecast systems, with the biggest improvements shown for predicting disease risk at mid- to high-disease levels. Most of the prediction uncertainty arises from model uncertainty, so prediction accuracy and precision do not improve substantially with shorter lead times. This result arises because many predictor variables cannot be accurately forecasted, which is a common challenge across ecosystems. Weekly forecasts and scenarios can be explored through an online decision support tool and data explorer, co-developed with end-user groups to improve use and understanding of ecological forecasts. The models provide near-real-time disease risk assessments and allow users to refine predictions and assess intervention scenarios. This work advances the field of ecological forecasting with real-world complexities and, in doing so, better supports near-term decision making for coral reef ecosystem managers and stakeholders. Secondarily, we identify clear needs and provide recommendations to further enhance our ability to forecast coral disease risk.

Inverse association between maternal serum concentrations of trace elements and risk of spontaneous preterm birth: a nested case-control study in China.

Few studies have evaluated the joint effect of trace elements on spontaneous preterm birth (SPTB). This study aimed to examine the relationships between the individual or mixed maternal serum concentrations of Fe, Cu, Zn, Se, Sr and Mo during pregnancy, and risk of SPTB. Inductively coupled plasma MS was employed to determine maternal serum concentrations of the six trace elements in 192 cases with SPTB and 282 controls with full-term delivery. Multivariate logistic regression, weighted quantile sum regression (WQSR) and Bayesian kernel machine regression (BKMR) were used to evaluate the individual and joint effects of trace elements on SPTB. The median concentrations of Sr and Mo were significantly higher in controls than in SPTB group (P < 0·05). In multivariate logistic regression analysis, compared with the lowest quartile levels of individual trace elements, the third- and fourth-quartile Sr or Mo concentrations were significantly associated with reduced risk of SPTB with adjusted OR (aOR) of 0·432 (95 CI < 0·05). In multivariate logistic regression analysis, compared with the lowest quartile levels of individual trace elements, the third- and fourth-quartile Sr or Mo concentrations were significantly associated with reduced risk of SPTB with adjusted aOR of 0·432 (95 % CI 0·247, 0·756), 0·386 (95 % CI 0·213, 0·701), 0·512 (95 % CI 0·297, 0·883) and 0·559 (95 % CI 0·321, 0·972), respectively. WQSR revealed the inverse combined effect of the trace elements mixture on SPTB (aOR = 0·368, 95 % CI 0·228, 0·593). BKMR analysis confirmed the overall mixture of the trace elements was inversely associated with the risk of SPTB, and the independent effect of Sr and Mo was significant. Our findings suggest that the risk of SPTB decreased with concentrations of the six trace elements, with Sr and Mo being the major contributors.

Smoothed quantile residual life regression analysis with application to the Korea HIV/AIDS cohort study.

BACKGROUND: The residual life of a patient with human immunodeficiency virus (HIV) is of major interest to patients and their physicians. While existing analyses of HIV patient survival focus mostly on data collected at baseline, residual life analysis allows for dynamic analysis based on additional data collected over a period of time. As survival times typically exhibit a right-skewed distribution, the median provides a more useful summary of the underlying distribution than the mean. In this paper, we propose an efficient inference procedure that fits a semiparametric quantile regression model assessing the effect of longitudinal biomarkers on the residual life of HIV patients until the development of dyslipidemia, a disease becoming more prevalent among those with HIV. METHODS: For estimation of model parameters, we propose an induced smoothing method that smooths nonsmooth estimating functions based on check functions. For variance estimation, we propose an efficient resampling-based estimator. The proposed estimators are theoretically justified. Simulation studies are used to evaluate their finite sample performances, including their prediction accuracy. We analyze the Korea HIV/AIDS cohort study data to examine the effects of CD4 (cluster of differentiation 4) cell count on the residual life of HIV patients to the onset of dyslipidemia. RESULTS: The proposed estimator is shown to be consistent and normally distributed asymptotically. Under various simulation settings, our estimates are approximately unbiased. Their variances estimates are close to the empirical variances and their computational efficiency is superior to that of the nonsmooth counterparts. Two measures of prediction performance indicate that our method adequately reflects the dynamic character of longitudinal biomarkers and residual life. The analysis of the Korea HIV/AIDS cohort study data shows that CD4 cell count is positively associated with residual life to the onset of dyslipidemia but the effect is not statistically significant. CONCLUSIONS: Our method enables direct prediction of residual lifetimes with a dynamic feature that accommodates data accumulated at different times. Our estimator significantly improves computational efficiency in variance estimation compared to the existing nonsmooth estimator. Analysis of the HIV/AIDS cohort study data reveals dynamic effects of CD4 cell count on the residual life to the onset of dyslipidemia.

Impact of Gestational Exposure to Individual and Combined Per- and Polyfluoroalkyl Substances on a Placental Structure and Efficiency: Findings from the Ma'anshan Birth Cohort.

Prenatal exposure to perfluoroalkyl and polyfluoroalkyl substances (PFASs) is inevitable among pregnant women. Nevertheless, there is a scarcity of research investigating the connections between prenatal PFAS exposure and the placental structure and efficiency. Based on 712 maternal-fetal dyads in the Ma'anshan Birth Cohort, we analyzed associations between individual and mixed PFAS exposure and placental measures. We repeatedly measured 12 PFAS in the maternal serum during pregnancy. Placental weight, scaling exponent, chorionic disc area, and disc eccentricity were used as the outcome variables. Upon adjusting for confounders and implementing corrections for multiple comparisons, we identified positive associations between branched perfluorohexane sulfonate (br-PFHxS) and 6:2 chlorinated polyfluorinated ether sulfonate (6:2 Cl-PFESA) with placental weight. Additionally, a positive association was observed between br-PFHxS and the scaling exponent, where a higher scaling exponent signified reduced placental efficiency. Based on neonatal sex stratification, female infants were found to be more susceptible to the adverse effects of PFAS exposure. Mixed exposure modeling revealed that mixed PFAS exposure was positively associated with placental weight and scaling exponent, particularly during the second and third trimesters. Furthermore, br-PFHxS and 6:2 Cl-PFESA played major roles in the placental measures. This study provides the first epidemiological evidence of the relationship between prenatal PFAS exposure and placental measures.

Individual and Neighborhood Level Predictors of Children's Exposure to Residential Greenspace.

Inequities in urban greenspace have been identified, though patterns by race and socioeconomic status vary across US settings. We estimated the magnitude of the relationship between a broad mixture of neighborhood-level factors and residential greenspace using weighted quantile sum (WQS) regression, and compared predictive models of greenspace using only neighborhood-level, only individual-level, or multi-level predictors. Greenspace measures included the Normalized Difference Vegetation Index (NDVI), tree canopy, and proximity of the nearest park, for residential locations in Shelby County, Tennessee of children in the CANDLE cohort. Neighborhood measures include socioeconomic and education resources, as well as racial composition and racial residential segregation. In this sample of 1012 mother-child dyads, neighborhood factors were associated with higher NDVI and tree canopy (0.021 unit higher NDVI [95% CI: 0.014, 0.028] per quintile increase in WQS index); homeownership rate, proximity of and enrollment at early childhood education centers, and racial composition, were highly weighted in the WQS index. In models constrained in the opposite direction (0.028 unit lower NDVI [95% CI: - 0.036, - 0.020]), high school graduation rate and teacher experience were highly weighted. In prediction models, adding individual-level predictors to the suite of neighborhood characteristics did not meaningfully improve prediction accuracy for greenspace measures. Our findings highlight disparities in greenspace for families by neighborhood socioeconomic and early education factors, and by race, suggesting several neighborhood indicators for consideration both as potential confounders in studies of greenspace and pediatric health as well as in the development of policies and programs to improve equity in greenspace access.

Association between fat-soluble vitamins and metabolic syndromes in US adults: a cross-section study from NHANES database.

BACKGROUND: Previous studies have shown significant associations between individual fat-soluble vitamins (FSVs) and metabolic syndromes (MetS). However, evidence on the multiple FSVs co-exposure and MetS odds is limited. Given that individuals are typically exposed to different levels of FSVs simultaneously, and FSVs can interact with each other. It's necessary to explore the association between multiple FSVs co-exposure and MetS odds. This study aims to address this gap in general U.S. adults aged ≥ 20 years. METHODS: We conducted a cross-sectional study utilizing data from the National Health and Nutrition Examination Surveys (NHANESs) 2003-2006 and 2017-2018. Three FSV, including vitamin A (VA), vitamin E (VE), and vitamin D (VD), and MetS diagnosed according to the ATP III guidelines were selected as exposure and outcome, respectively. Multivariable-adjusted logistic model was used to explore the associations of individual FSV exposure with MetS odds and MetS components. Restricted cubic splines were performed to explore the dose-response relationships among them. The quantile g-computation method was adopted to explore the associations of multiple FSVs co-exposure with MetS odds and MetS components. RESULTS: The presented study included a total of 13,975 individuals, with 2400 (17.17%) were diagnosed with MetS. After adjusting for various confounders, a positive linear pattern was observed for serum VA and VE and MetS associations. Serum VD was found to be negatively associated with MetS in a linear dose-response way. For each component of MetS, higher serum VA and VE were associated with higher triglyceride and high-density lipoprotein; higher serum VD was negatively associated with triglyceride, blood pressure, and fasting plasma glucose. MetS odds increased by 15% and 13%, respectively, in response to one quartile increase in FSVs co-exposure index (qgcomp) in the conditional model (OR = 1.15, 95%CI: 1.06, 1.24) and the marginal structural model (OR = 1.13, 95%CI: 1.06, 1.20). Besides, co-exposure to VA, VE, and VD was positively associated with triglyceride, high-density lipoprotein, and blood pressure levels. CONCLUSION: Findings in the present study revealed that high serum VA and VE levels were associated with elevated MetS odds, while serum VD was inversely associated with MetS odds. FSVs co-exposure was positively associated with MetS odds.