Publications by Year: 2023

Applying causal inference methods, such as weighting and matching methods, to a survey sampled population requires properly incorporating the survey weights and design to obtain effect estimates that are representative of the target population and correct standard errors (SEs). With a simulation study, we compared various approaches for incorporating the survey weights and design into weighting and matching-based causal inference methods. When the models were correctly specified, most approaches performed well. However, when a variable was treated as an unmeasured confounder and the survey weights were constructed to depend on this variable, only the matching methods that used the survey weights in causal estimation and as a covariate in matching continued to perform well. If unmeasured confounders are potentially associated with the survey sample design, we recommend that investigators include the survey weights as a covariate in matching, in addition to incorporating them in causal effect estimation. Finally, we applied the various approaches to the Hispanic Community Health Study/Study of Latinos (HCHS/SOL) and found that insomnia has a causal association with both mild cognitive impairment (MCI) and incident hypertension 6-7 years later in the US Hispanic/Latino population.

Both short (≤6 h per night) and long sleep duration (≥9 h per night) are associated with increased risk of chronic diseases. Despite evidence linking habitual sleep duration and risk of disease, the genetic determinants of sleep duration in the general population are poorly understood, especially outside of European (EUR) populations. Here, we report that a polygenic score of 78 European ancestry sleep duration single-nucleotide polymorphisms (SNPs) is associated with sleep duration in an African (n = 7288; P = 0.003), an East Asian (n = 13 618; P = 6 × 10-4) and a South Asian (n = 7485; P = 0.025) genetic ancestry cohort, but not in a Hispanic/Latino cohort (n = 8726; P = 0.71). Furthermore, in a pan-ancestry (N = 483 235) meta-analysis of genome-wide association studies (GWAS) for habitual sleep duration, 73 loci are associated with genome-wide statistical significance. Follow-up of five loci (near HACD2, COG5, PRR12, SH3RF1 and KCNQ5) identified expression-quantitative trait loci for PRR12 and COG5 in brain tissues and pleiotropic associations with cardiovascular and neuropsychiatric traits. Overall, our results suggest that the genetic basis of sleep duration is at least partially shared across diverse ancestry groups.

Rationale: Obstructive sleep apnea is characterized by frequent reductions in ventilation, leading to oxygen desaturations and/or arousals. Objectives: In this study, association of hypoxic burden with incident cardiovascular disease (CVD) was examined and compared with that of "ventilatory burden" and "arousal burden." Finally, we assessed the extent to which the ventilatory burden, visceral obesity, and lung function explain variations in hypoxic burden. Methods: Hypoxic, ventilatory, and arousal burdens were measured from baseline polysomnograms in the Multi-Ethnic Study of Atherosclerosis (MESA) and the Osteoporotic Fractures in Men (MrOS) studies. Ventilatory burden was defined as event-specific area under ventilation signal (mean normalized, area under the mean), and arousal burden was defined as the normalized cumulative duration of all arousals. The adjusted hazard ratios for incident CVD and mortality were calculated. Exploratory analyses quantified contributions to hypoxic burden of ventilatory burden, baseline oxygen saturation as measured by pulse oximetry, visceral obesity, and spirometry parameters. Measurements and Main Results: Hypoxic and ventilatory burdens were significantly associated with incident CVD (adjusted hazard ratio [95% confidence interval] per 1 SD increase in hypoxic burden: MESA, 1.45 [1.14, 1.84]; MrOS, 1.13 [1.02, 1.26]; ventilatory burden: MESA, 1.38 [1.11, 1.72]; MrOS, 1.12 [1.01, 1.25]), whereas arousal burden was not. Similar associations with mortality were also observed. Finally, 78% of variation in hypoxic burden was explained by ventilatory burden, whereas other factors explained only <2% of variation. Conclusions: Hypoxic and ventilatory burden predicted CVD morbidity and mortality in two population-based studies. Hypoxic burden is minimally affected by measures of adiposity and captures the risk attributable to ventilatory burden of obstructive sleep apnea rather than a tendency to desaturate.

Sleep-disordered breathing (SDB) is a prevalent disorder characterized by recurrent episodic upper airway obstruction. In a dataset from the Hispanic Community Health Study/Study of Latinos (HCHS/SOL), we applied principal component analysis (PCA) on seven measures characterizing SDB-associated respiratory events. We estimated the association of the top two SDB PCs with serum levels of 617 metabolites, in both single-metabolite analysis, and a joint, penalized regression analysis using the least absolute shrinkage and selection operator (LASSO). Discovery analysis included n = 3,299 HCHS/SOL individuals; associations were validated in a separate dataset of n = 1,522 HCHS/SOL individuals. Seven metabolite associations with SDB PCs were discovered and replicated. Metabolite risk scores (MRSs) developed based on LASSO association results and representing metabolite signatures associated with the two SDB PCs were associated with 6-year incident hypertension and incident diabetes. MRSs have the potential to serve as biomarkers for SDB, guiding risk stratification and treatment decisions.

STUDY OBJECTIVES: Rest-activity rhythms (RAR) may mark development, aging, and physical and mental health. Understanding how they differ between people may inform intervention and health promotion efforts. However, RAR characteristics across the lifespan have not been well-studied. Therefore, we investigated the association between RAR measures with demographic and lifestyle factors in a US nationally representative study.

METHODS: RAR metrics of interdaily stability (IS), intradaily variability (IV), relative amplitude (RA), and mean amplitude and timing of high (M10) and low (L5) activity were derived from 2011 to 2012 and 2013 to 2014 National Health and Nutrition Examination Survey (NHANES) actigraphy data. Population-weighted linear and logistic regression models were fit to examine the associations of age, gender, smoking, alcohol, season, body mass index (BMI), income-to-poverty ratio, and race/ethnicity with RAR. Significance was based on a false-discovery rate-corrected P-value of <0.05.

RESULTS: Among n = 12 526 NHANES participants (3-≥80 years), IS (higher = greater day-to-day regularity) and RA (higher = greater rhythm strength) generally decreased with age and were lower among males, whereas IV (higher = greater rhythm fragmentation) increased with age (p < 0.05). Dynamic changes in RAR trajectories were observed during childhood and adolescence. Income, BMI, smoking, and alcohol use were associated with RAR metrics, as well as season among children and teenagers (p < 0.05). RAR also differed by race/ethnicity (p < 0.05), with trajectories initially diverging in childhood and continuing into adulthood.

CONCLUSIONS: RAR differed by demographic and health-related factors, representing possible windows for public health intervention and sleep health promotion. RAR differences by race/ethnicity begin in childhood, are evident in early adolescence, and persist throughout adulthood.

BACKGROUND: Identifying novel epigenetic signatures associated with serum immunoglobulin E (IgE) may improve our understanding of molecular mechanisms underlying asthma and IgE-mediated diseases.

METHODS: We performed an epigenome-wide association study using whole blood from Framingham Heart Study (FHS; n = 3,471, 46% females) participants and validated results using the Childhood Asthma Management Program (CAMP; n = 674, 39% females) and the Genetic Epidemiology of Asthma in Costa Rica Study (CRA; n = 787, 41% females). Using the closest gene to each IgE-associated CpG, we highlighted biologically plausible pathways underlying IgE regulation and analyzed the transcription patterns linked to IgE-associated CpGs (expression quantitative trait methylation loci; eQTMs). Using prior UK Biobank summary data from genome-wide association studies of asthma and allergy, we performed Mendelian randomization (MR) for causal inference testing using the IgE-associated CpGs from FHS with methylation quantitative trait loci (mQTLs) as instrumental variables.

FINDINGS: We identified 490 statistically significant differentially methylated CpGs associated with IgE in FHS, of which 193 (39.3%) replicated in CAMP and CRA (FDR < 0.05). Gene ontology analysis revealed enrichment in pathways related to transcription factor binding, asthma, and other immunological processes. eQTM analysis identified 124 cis-eQTMs for 106 expressed genes (FDR < 0.05). MR in combination with drug-target analysis revealed CTSB and USP20 as putatively causal regulators of IgE levels (Bonferroni adjusted P < 7.94E-04) that can be explored as potential therapeutic targets.

INTERPRETATION: By integrating eQTM and MR analyses in general and clinical asthma populations, our findings provide a deeper understanding of the multidimensional inter-relations of DNA methylation, gene expression, and IgE levels.

FUNDING: US NIH/NHLBI grants: P01HL132825, K99HL159234. N01-HC-25195 and HHSN268201500001I.

INTRODUCTION: Polygenic Risk Scores (PRSs) are summaries of genetic risk alleles for an outcome.

METHODS: We used summary statistics from five GWASs of AD to construct PRSs in 4,189 diverse Hispanics/Latinos (mean age 63 years) from the Study of Latinos-Investigation of Neurocognitive Aging (SOL-INCA). We assessed the PRS associations with MCI in the combined set of people and in diverse subgroups, and when including and excluding the APOE gene region. We also assessed PRS associations with MCI in an independent dataset from the Mass General Brigham Biobank.

RESULTS: A simple sum of 5 PRSs ("PRSsum"), each constructed based on a different AD GWAS, was associated with MCI (OR = 1.28, 95% CI [1.14, 1.41]) in a model adjusted for counts of the APOE-[Formula: see text] and APOE-[Formula: see text] alleles. Associations of single-GWAS PRSs were weaker. When removing SNPs from the APOE region from the PRSs, the association of PRSsum with MCI was weaker (OR = 1.17, 95% CI [1.04,1.31] with adjustment for APOE alleles). In all association analyses, APOE-[Formula: see text] and APOE-[Formula: see text] alleles were not associated with MCI.

DISCUSSION: A sum of AD PRSs is associated with MCI in Hispanic/Latino older adults. Despite no association of APOE-[Formula: see text] and APOE-[Formula: see text] alleles with MCI, the association of the AD PRS with MCI is stronger when including the APOE region. Thus, APOE variants different than the classic APOE alleles may be important predictors of MCI in Hispanic/Latino adults.

OBJECTIVE: Exposure to light at night (LAN) may influence sleep timing and regularity. Here, we test whether greater light exposure during sleep (LEDS) associates with greater irregularity in sleep onset timing in a large cohort of older adults.

METHODS: Light exposure and activity patterns, measured via wrist-worn actigraphy (ActiWatch Spectrum), were analyzed in 1,933 participants with 6+ valid days of data in the Multi-Ethnic Study of Atherosclerosis (MESA) Exam 5 Sleep Study. Summary measures of LEDS averaged across nights were evaluated in linear and logistic regression analyses to test the association with standard deviation (SD) in sleep onset timing (continuous variable) and irregular sleep onset timing (SD≥1.36 hours, binary). Night-to-night associations between LEDS and absolute differences in nightly sleep onset timing were also evaluated with distributed lag non-linear models and mixed models.

RESULTS: In between-individual linear and logistic models adjusted for demographic, health, and seasonal factors, every 5-lux unit increase in LEDS was associated with an increase of 7.8 minutes in sleep onset SD (β=0.13 hours, 95%CI:0.09-0.17) and 40% greater odds (OR=1.40, 95%CI:1.24-1.60) of irregular sleep onset. In within-individual night-to-night mixed model analyses, every 5-lux unit increase in LEDS the night prior (lag0) was associated with a 2.2-minute greater deviation of sleep onset the next night (β=0.036 hours, p<0.05). Conversely, every 1-hour increase in sleep deviation (lag0) was associated with a 0.35-lux increase in future LEDS (β=0.347 lux, p<0.05).

CONCLUSION: LEDS was associated with greater irregularity in sleep onset in between-individual analyses and subsequent deviation in sleep timing in within-individual analyses, supporting a role for LEDS in exacerbating irregular sleep onset timing. Greater deviation in sleep onset was also associated with greater future LEDS, suggesting a bidirectional relationship. Maintaining a dark sleeping environment and preventing LEDS may promote sleep regularity and following a regular sleep schedule may limit LEDS.