Novel EDGE encoding method enhances ability to identify genetic interactions.

Assumptions are made about the genetic model of single nucleotide polymorphisms (SNPs) when choosing a traditional genetic encoding: additive, dominant, and recessive. Furthermore, SNPs across the genome are unlikely to demonstrate identical genetic models. However, running SNP-SNP interaction analyses with every combination of encodings raises the multiple testing burden. Here, we present a novel and flexible encoding for genetic interactions, the elastic data-driven genetic encoding (EDGE), in which SNPs are assigned a heterozygous value based on the genetic model they demonstrate in a dataset prior to interaction testing. We assessed the power of EDGE to detect genetic interactions using 29 combinations of simulated genetic models and found it outperformed the traditional encoding methods across 10%, 30%, and 50% minor allele frequencies (MAFs). Further, EDGE maintained a low false-positive rate, while additive and dominant encodings demonstrated inflation. We evaluated EDGE and the traditional encodings with genetic data from the Electronic Medical Records and Genomics (eMERGE) Network for five phenotypes: age-related macular degeneration (AMD), age-related cataract, glaucoma, type 2 diabetes (T2D), and resistant hypertension. A multi-encoding genome-wide association study (GWAS) for each phenotype was performed using the traditional encodings, and the top results of the multi-encoding GWAS were considered for SNP-SNP interaction using the traditional encodings and EDGE. EDGE identified a novel SNP-SNP interaction for age-related cataract that no other method identified: rs7787286 (MAF: 0.041; intergenic region of chromosome 7)-rs4695885 (MAF: 0.34; intergenic region of chromosome 4) with a Bonferroni LRT p of 0.018. A SNP-SNP interaction was found in data from the UK Biobank within 25 kb of these SNPs using the recessive encoding: rs60374751 (MAF: 0.030) and rs6843594 (MAF: 0.34) (Bonferroni LRT p: 0.026). We recommend using EDGE to flexibly detect interactions between SNPs exhibiting diverse action.

Dynamic assessment of the COVID-19 vaccine acceptance leveraging social media data.

Vaccination is the most effective way to provide long-lasting immunity against viral infection; thus, rapid assessment of vaccine acceptance is a pressing challenge for health authorities. Prior studies have applied survey techniques to investigate vaccine acceptance, but these may be slow and expensive. This study investigates 29 million vaccine-related tweets from August 8, 2020 to April 19, 2021 and proposes a social media-based approach that derives a vaccine acceptance index (VAI) to quantify Twitter users' opinions on COVID-19 vaccination. This index is calculated based on opinion classifications identified with the aid of natural language processing techniques and provides a quantitative metric to indicate the level of vaccine acceptance across different geographic scales in the U.S. The VAI is easily calculated from the number of positive and negative Tweets posted by a specific users and groups of users, it can be compiled for regions such a counties or states to provide geospatial information, and it can be tracked over time to assess changes in vaccine acceptance as related to trends in the media and politics. At the national level, it showed that the VAI moved from negative to positive in 2020 and maintained steady after January 2021. Through exploratory analysis of state- and county-level data, reliable assessments of VAI against subsequent vaccination rates could be made for counties with at least 30 users. The paper discusses information characteristics that enable consistent estimation of VAI. The findings support the use of social media to understand opinions and to offer a timely and cost-effective way to assess vaccine acceptance.

Embracing Complex Associations in Common Traits: Critical Considerations for Precision Medicine.

Genome-wide association studies (GWAS) have identified numerous loci associated with human phenotypes. This approach, however, does not consider the richly diverse and complex environment with which humans interact throughout the life course, nor does it allow for interrelationships between genetic loci and across traits. As we move toward making precision medicine a reality, whereby we make predictions about disease risk based on genomic profiles, we need to identify improved predictive models of the relationship between genome and phenome. Methods that embrace pleiotropy (the effect of one locus on more than one trait), and gene-environment (G×E) and gene-gene (G×G) interactions, will further unveil the impact of alterations in biological pathways and identify genes that are only involved with disease in the context of the environment. This valuable information can be used to assess personal risk and choose the most appropriate medical interventions based on the genotype and environment of an individual, the whole premise of precision medicine.

Identifying gene-gene interactions that are highly associated with four quantitative lipid traits across multiple cohorts.

Genetic loci explain only 25-30 % of the heritability observed in plasma lipid traits. Epistasis, or gene-gene interactions may contribute to a portion of this missing heritability. Using the genetic data from five NHLBI cohorts of 24,837 individuals, we combined the use of the quantitative multifactor dimensionality reduction (QMDR) algorithm with two SNP-filtering methods to exhaustively search for SNP-SNP interactions that are associated with HDL cholesterol (HDL-C), LDL cholesterol (LDL-C), total cholesterol (TC) and triglycerides (TG). SNPs were filtered either on the strength of their independent effects (main effect filter) or the prior knowledge supporting a given interaction (Biofilter). After the main effect filter, QMDR identified 20 SNP-SNP models associated with HDL-C, 6 associated with LDL-C, 3 associated with TC, and 10 associated with TG (permutation P value <0.05). With the use of Biofilter, we identified 2 SNP-SNP models associated with HDL-C, 3 associated with LDL-C, 1 associated with TC and 8 associated with TG (permutation P value <0.05). In an independent dataset of 7502 individuals from the eMERGE network, we replicated 14 of the interactions identified after main effect filtering: 11 for HDL-C, 1 for LDL-C and 2 for TG. We also replicated 23 of the interactions found to be associated with TG after applying Biofilter. Prior knowledge supports the possible role of these interactions in the genetic etiology of lipid traits. This study also presents a computationally efficient pipeline for analyzing data from large genotyping arrays and detecting SNP-SNP interactions that are not primarily driven by strong main effects.

Informatics and Data Analytics to Support Exposome-Based Discovery for Public Health.

The complexity of the human exposome-the totality of environmental exposures encountered from birth to death-motivates systematic, high-throughput approaches to discover new environmental determinants of disease. In this review, we describe the state of science in analyzing the human exposome and provide recommendations for the public health community to consider in dealing with analytic challenges of exposome-based biomedical research. We describe extant and novel analytic methods needed to associate the exposome with critical health outcomes and contextualize the data-centered challenges by drawing parallels to other research endeavors such as human genomics research. We discuss efforts for training scientists who can bridge public health, genomics, and biomedicine in informatics and statistics. If an exposome data ecosystem is brought to fruition, it will likely play a role as central as genomic science has had in molding the current and new generations of biomedical researchers, computational scientists, and public health research programs.

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.

We performed a Phenome-wide association study (PheWAS) utilizing diverse genotypic and phenotypic data existing across multiple populations in the National Health and Nutrition Examination Surveys (NHANES), conducted by the Centers for Disease Control and Prevention (CDC), and accessed by the Epidemiological Architecture for Genes Linked to Environment (EAGLE) study. We calculated comprehensive tests of association in Genetic NHANES using 80 SNPs and 1,008 phenotypes (grouped into 184 phenotype classes), stratified by race-ethnicity. Genetic NHANES includes three surveys (NHANES III, 1999-2000, and 2001-2002) and three race-ethnicities: non-Hispanic whites (n = 6,634), non-Hispanic blacks (n = 3,458), and Mexican Americans (n = 3,950). We identified 69 PheWAS associations replicating across surveys for the same SNP, phenotype-class, direction of effect, and race-ethnicity at p<0.01, allele frequency >0.01, and sample size >200. Of these 69 PheWAS associations, 39 replicated previously reported SNP-phenotype associations, 9 were related to previously reported associations, and 21 were novel associations. Fourteen results had the same direction of effect across more than one race-ethnicity: one result was novel, 11 replicated previously reported associations, and two were related to previously reported results. Thirteen SNPs showed evidence of pleiotropy. We further explored results with gene-based biological networks, contrasting the direction of effect for pleiotropic associations across phenotypes. One PheWAS result was ABCG2 missense SNP rs2231142, associated with uric acid levels in both non-Hispanic whites and Mexican Americans, protoporphyrin levels in non-Hispanic whites and Mexican Americans, and blood pressure levels in Mexican Americans. Another example was SNP rs1800588 near LIPC, significantly associated with the novel phenotypes of folate levels (Mexican Americans), vitamin E levels (non-Hispanic whites) and triglyceride levels (non-Hispanic whites), and replication for cholesterol levels. The results of this PheWAS show the utility of this approach for exposing more of the complex genetic architecture underlying multiple traits, through generating novel hypotheses for future research.

Biology-Driven Gene-Gene Interaction Analysis of Age-Related Cataract in the eMERGE Network.

Bioinformatics approaches to examine gene-gene models provide a means to discover interactions between multiple genes that underlie complex disease. Extensive computational demands and adjusting for multiple testing make uncovering genetic interactions a challenge. Here, we address these issues using our knowledge-driven filtering method, Biofilter, to identify putative single nucleotide polymorphism (SNP) interaction models for cataract susceptibility, thereby reducing the number of models for analysis. Models were evaluated in 3,377 European Americans (1,185 controls, 2,192 cases) from the Marshfield Clinic, a study site of the Electronic Medical Records and Genomics (eMERGE) Network, using logistic regression. All statistically significant models from the Marshfield Clinic were then evaluated in an independent dataset of 4,311 individuals (742 controls, 3,569 cases), using independent samples from additional study sites in the eMERGE Network: Mayo Clinic, Group Health/University of Washington, Vanderbilt University Medical Center, and Geisinger Health System. Eighty-three SNP-SNP models replicated in the independent dataset at likelihood ratio test P < 0.05. Among the most significant replicating models was rs12597188 (intron of CDH1)-rs11564445 (intron of CTNNB1). These genes are known to be involved in processes that include: cell-to-cell adhesion signaling, cell-cell junction organization, and cell-cell communication. Further Biofilter analysis of all replicating models revealed a number of common functions among the genes harboring the 83 replicating SNP-SNP models, which included signal transduction and PI3K-Akt signaling pathway. These findings demonstrate the utility of Biofilter as a biology-driven method, applicable for any genome-wide association study dataset.

Phenome-Wide Association Studies: Embracing Complexity for Discovery.

The inherent complexity of biological systems can be leveraged for a greater understanding of the impact of genetic architecture on outcomes, traits, and pharmacological response. The genome-wide association study (GWAS) approach has well-developed methods and relatively straight-forward methodologies; however, the bigger picture of the impact of genetic architecture on phenotypic outcome still remains to be elucidated even with an ever-growing number of GWAS performed. Greater consideration of the complexity of biological processes, using more data from the phenome, exposome, and diverse -omic resources, including considering the interplay of pleiotropy and genetic interactions, may provide additional leverage for making the most of the incredible wealth of information available for study. Here, we describe how incorporating greater complexity into analyses through the use of additional phenotypic data and widespread deployment of phenome-wide association studies may provide new insights into genetic factors influencing diseases, traits, and pharmacological response.

Electronic medical records and genomics (eMERGE) network exploration in cataract: several new potential susceptibility loci.

PURPOSE: Cataract is the leading cause of blindness in the world, and in the United States accounts for approximately 60% of Medicare costs related to vision. The purpose of this study was to identify genetic markers for age-related cataract through a genome-wide association study (GWAS). METHODS: In the electronic medical records and genomics (eMERGE) network, we ran an electronic phenotyping algorithm on individuals in each of five sites with electronic medical records linked to DNA biobanks. We performed a GWAS using 530,101 SNPs from the Illumina 660W-Quad in a total of 7,397 individuals (5,503 cases and 1,894 controls). We also performed an age-at-diagnosis case-only analysis. RESULTS: We identified several statistically significant associations with age-related cataract (45 SNPs) as well as age at diagnosis (44 SNPs). The 45 SNPs associated with cataract at p<1×10(-5) are in several interesting genes, including ALDOB, MAP3K1, and MEF2C. All have potential biologic relationships with cataracts. CONCLUSIONS: This is the first genome-wide association study of age-related cataract, and several regions of interest have been identified. The eMERGE network has pioneered the exploration of genomic associations in biobanks linked to electronic health records, and this study is another example of the utility of such resources. Explorations of age-related cataract including validation and replication of the association results identified herein are needed in future studies.

Sex modifies effects of imaging and CSF biomarkers on cognitive and functional outcomes: a study of Alzheimer's disease.

Alzheimer's disease (AD) is a neurodegenerative disorder characterized by memory and functional impairments. Two of 3 patients with AD are biologically female; therefore, the biological underpinnings of this diagnosis disparity may inform interventions slowing the AD progression. To bridge this gap, we conducted analyses of 1078 male and female participants from the Alzheimer's Disease Neuroimaging Initiative to examine associations between levels of cerebral spinal fluid (CSF)/neuroimaging biomarkers and cognitive/functional outcomes. The Chow test was used to quantify sex differences by determining if biological sex affects relationships between the studied biomarkers and outcomes. Multiple magnetic resonance imaging (whole brain, entorhinal cortex, middle temporal gyrus, fusiform gyrus, hippocampus), position emission tomography (AV45), and CSF (P-TAU, TAU) biomarkers were differentially associated with cognitive and functional outcomes. Post-hoc bootstrapped and association analyses confirmed these differential effects and emphasized the necessity of using separate, sex-stratified models. The studied imaging/CSF biomarkers may account for some of the sex-based variation in AD pathophysiology. The identified sex-varying relationships between CSF/imaging biomarkers and cognitive/functional outcomes warrant future biological investigation in independent cohorts.

Longitudinal method comparison: modeling polygenic risk for post-traumatic stress disorder over time in individuals of African and European ancestry.

Cross-sectional data allow the investigation of how genetics influence health at a single time point, but to understand how the genome impacts phenotype development, one must use repeated measures data. Ignoring the dependency inherent in repeated measures can exacerbate false positives and requires the utilization of methods other than general or generalized linear models. Many methods can accommodate longitudinal data, including the commonly used linear mixed model and generalized estimating equation, as well as the less popular fixed-effects model, cluster-robust standard error adjustment, and aggregate regression. We simulated longitudinal data and applied these five methods alongside naïve linear regression, which ignored the dependency and served as a baseline, to compare their power, false positive rate, estimation accuracy, and precision. The results showed that the naïve linear regression and fixed-effects models incurred high false positive rates when analyzing a predictor that is fixed over time, making them unviable for studying time-invariant genetic effects. The linear mixed models maintained low false positive rates and unbiased estimation. The generalized estimating equation was similar to the former in terms of power and estimation, but it had increased false positives when the sample size was low, as did cluster-robust standard error adjustment. Aggregate regression produced biased estimates when predictor effects varied over time. To show how the method choice affects downstream results, we performed longitudinal analyses in an adolescent cohort of African and European ancestry. We examined how developing post-traumatic stress symptoms were predicted by polygenic risk, traumatic events, exposure to sexual abuse, and income using four approaches-linear mixed models, generalized estimating equations, cluster-robust standard error adjustment, and aggregate regression. While the directions of effect were generally consistent, coefficient magnitudes and statistical significance differed across methods. Our in-depth comparison of longitudinal methods showed that linear mixed models and generalized estimating equations were applicable in most scenarios requiring longitudinal modeling, but no approach produced identical results even if fit to the same data. Since result discrepancies can result from methodological choices, it is crucial that researchers determine their model a priori, refrain from testing multiple approaches to obtain favorable results, and utilize as similar as possible methods when seeking to replicate results.

Neuronally enriched microvesicle RNAs are differentially expressed in the serums of Parkinson's patients.

Background: Circulating small RNAs (smRNAs) originate from diverse tissues and organs. Previous studies investigating smRNAs as potential biomarkers for Parkinson's disease (PD) have yielded inconsistent results. We investigated whether smRNA profiles from neuronally-enriched serum exosomes and microvesicles are altered in PD patients and discriminate PD subjects from controls. Methods: Demographic, clinical, and serum samples were obtained from 60 PD subjects and 40 age- and sex-matched controls. Exosomes and microvesicles were extracted and isolated using a validated neuronal membrane marker (CD171). Sequencing and bioinformatics analyses were used to identify differentially expressed smRNAs in PD and control samples. SmRNAs also were tested for association with clinical metrics. Logistic regression and random forest classification models evaluated the discriminative value of the smRNAs. Results: In serum CD171 enriched exosomes and microvesicles, a panel of 29 smRNAs was expressed differentially between PD and controls (false discovery rate (FDR) < 0.05). Among the smRNAs, 23 were upregulated and 6 were downregulated in PD patients. Pathway analysis revealed links to cellular proliferation regulation and signaling. Least absolute shrinkage and selection operator adjusted for the multicollinearity of these smRNAs and association tests to clinical parameters via linear regression did not yield significant results. Univariate logistic regression models showed that four smRNAs achieved an AUC ≥ 0.74 to discriminate PD subjects from controls. The random forest model had an AUC of 0.942 for the 29 smRNA panel. Conclusion: CD171-enriched exosomes and microvesicles contain the differential expression of smRNAs between PD and controls. Future studies are warranted to follow up on the findings and understand the scientific and clinical relevance.

Loss of selenoprotein W in murine macrophages alters the hierarchy of selenoprotein expression, redox tone, and mitochondrial functions during inflammation.

Macrophages play a pivotal role in mediating inflammation and subsequent resolution of inflammation. The availability of selenium as a micronutrient and the subsequent biosynthesis of selenoproteins, containing the 21st amino acid selenocysteine (Sec), are important for the physiological functions of macrophages. Selenoproteins regulate the redox tone in macrophages during inflammation, the early onset of which involves oxidative burst of reactive oxygen and nitrogen species. SELENOW is a highly expressed selenoprotein in bone marrow-derived macrophages (BMDMs). Beyond its described general role as a thiol and peroxide reductase and as an interacting partner for 14-3-3 proteins, its cellular functions, particularly in macrophages, remain largely unknown. In this study, we utilized Selenow knock-out (KO) murine bone marrow-derived macrophages (BMDMs) to address the role of SELENOW in inflammation following stimulation with bacterial endotoxin lipopolysaccharide (LPS). RNAseq-based temporal analyses of expression of selenoproteins and the Sec incorporation machinery genes suggested no major differences in the selenium utilization pathway in the Selenow KO BMDMs compared to their wild-type counterparts. However, selective enrichment of oxidative stress-related selenoproteins and increased ROS in Selenow-/- BMDMs indicated anomalies in redox homeostasis associated with hierarchical expression of selenoproteins. Selenow-/- BMDMs also exhibited reduced expression of arginase-1, a key enzyme associated with anti-inflammatory (M2) phenotype necessary to resolve inflammation, along with a significant decrease in efferocytosis of neutrophils that triggers pathways of resolution. Parallel targeted metabolomics analysis also confirmed an impairment in arginine metabolism in Selenow-/- BMDMs. Furthermore, Selenow-/- BMDMs lacked the ability to enhance characteristic glycolytic metabolism during inflammation. Instead, these macrophages atypically relied on oxidative phosphorylation for energy production when glucose was used as an energy source. These findings suggest that SELENOW expression in macrophages may have important implications on cellular redox processes and bioenergetics during inflammation and its resolution.

Nitric oxide regulates metabolism in murine stress erythroid progenitors to promote recovery during inflammatory anemia.

Inflammation skews bone marrow hematopoiesis increasing the production of myeloid effector cells at the expense of steady-state erythropoiesis. A compensatory stress erythropoiesis response is induced to maintain homeostasis until inflammation is resolved. In contrast to steady-state erythroid progenitors, stress erythroid progenitors (SEPs) utilize signals induced by inflammatory stimuli. However, the mechanistic basis for this is not clear. Here we reveal a nitric oxide (NO)-dependent regulatory network underlying two stages of stress erythropoiesis, namely proliferation, and the transition to differentiation. In the proliferative stage, immature SEPs and cells in the niche increased expression of inducible nitric oxide synthase ( Nos2 or iNOS ) to generate NO. Increased NO rewires SEP metabolism to increase anabolic pathways, which drive the biosynthesis of nucleotides, amino acids and other intermediates needed for cell division. This NO-dependent metabolism promotes cell proliferation while also inhibiting erythroid differentiation leading to the amplification of a large population of non-committed progenitors. The transition of these progenitors to differentiation is mediated by the activation of nuclear factor erythroid 2-related factor 2 (Nfe2l2 or Nrf2). Nrf2 acts as an anti-inflammatory regulator that decreases NO production, which removes the NO-dependent erythroid inhibition and allows for differentiation. These data provide a paradigm for how alterations in metabolism allow inflammatory signals to amplify immature progenitors prior to differentiation. Key points: Nitric-oxide (NO) dependent signaling favors an anabolic metabolism that promotes proliferation and inhibits differentiation.Activation of Nfe2l2 (Nrf2) decreases NO production allowing erythroid differentiation.

Activation of GPR44 decreases severity of myeloid leukemia via specific targeting of leukemia initiating stem cells.

Relapse of acute myeloid leukemia (AML) remains a significant concern due to persistent leukemia-initiating stem cells (LICs) that are typically not targeted by most existing therapies. Using a murine AML model, human AML cell lines, and patient samples, we show that AML LICs are sensitive to endogenous and exogenous cyclopentenone prostaglandin-J (CyPG), Δ12-PGJ2, and 15d-PGJ2, which are increased upon dietary selenium supplementation via the cyclooxygenase-hematopoietic PGD synthase pathway. CyPGs are endogenous ligands for peroxisome proliferator-activated receptor gamma and GPR44 (CRTH2; PTGDR2). Deletion of GPR44 in a mouse model of AML exacerbated the disease suggesting that GPR44 activation mediates selenium-mediated apoptosis of LICs. Transcriptomic analysis of GPR44-/- LICs indicated that GPR44 activation by CyPGs suppressed KRAS-mediated MAPK and PI3K/AKT/mTOR signaling pathways, to enhance apoptosis. Our studies show the role of GPR44, providing mechanistic underpinnings of the chemopreventive and chemotherapeutic properties of selenium and CyPGs in AML.

In Situ Thermoset Cure Sensing: A Review of Correlation Methods.

Thermoset polymer composites have increased in use across multiple industries, with recent applications consisting of high-complexity and large-scale parts. As applications expand, the emphasis on accurate process-monitoring techniques has increased, with a variety of in situ cure-monitoring sensors being investigated by various research teams. To date, a wide range of data analysis techniques have been used to correlate data collected from thermocouple, dielectric, ultrasonic, and fibre-optic sensors to information on the material cure state. The methods used in existing publications have not been explicitly differentiated between, nor have they been directly compared. This paper provides a critical review of the different data collection and cure state correlation methods for these sensor types. The review includes details of the relevant sensor configurations and governing equations, material combinations, data verification techniques, identified potential research gaps, and areas of improvement. A wide range of both qualitative and quantitative analysis methods are discussed for each sensing technology. Critical analysis is provided on the capability and limitations of these methods to directly identify cure state information for the materials under investigation. This paper aims to provide the reader with sufficient background on available analysis techniques to assist in selecting the most appropriate method for the application.

Sex Differences in the Metabolome of Alzheimer's Disease Progression.

Alzheimer's disease (AD) is the leading cause of dementia; however, men and women face differential AD prevalence, presentation, and progression risks. Characterizing metabolomic profiles during AD progression is fundamental to understand the metabolic disruptions and the biological pathways involved. However, outstanding questions remain of whether peripheral metabolic changes occur equally in men and women with AD. Here, we evaluated differential effects of metabolomic and brain volume associations between sexes. We used three cohorts from the Alzheimer's Disease Neuroimaging Initiative (ADNI), evaluated 1,368 participants, two metabolomic platforms with 380 metabolites in total, and six brain segment volumes. Using dimension reduction techniques, we took advantage of the correlation structure of the brain volume phenotypes and the metabolite concentration values to reduce the number of tests while aggregating relevant biological structures. Using WGCNA, we aggregated modules of highly co-expressed metabolites. On the other hand, we used partial least squares regression-discriminant analysis (PLS-DA) to extract components of brain volumes that maximally co-vary with AD diagnosis as phenotypes. We tested for differences in effect sizes between sexes in the association between single metabolite and metabolite modules with the brain volume components. We found five metabolite modules and 125 single metabolites with significant differences between sexes. These results highlight a differential lipid disruption in AD progression between sexes. Men showed a greater negative association of phosphatidylcholines and sphingomyelins and a positive association of VLDL and large LDL with AD progression. In contrast, women showed a positive association of triglycerides in VLDL and small and medium LDL with AD progression. Explicitly identifying sex differences in metabolomics during AD progression can highlight particular metabolic disruptions in each sex. Our research study and strategy can lead to better-tailored studies and better-suited treatments that take sex differences into account.

Obesity and accelerated epigenetic aging in a high-risk cohort of children.

New insights into mechanisms linking obesity to poor health outcomes suggest a role for cellular aging pathways, casting obesity as a disease of accelerated biological aging. Although obesity has been linked to accelerated epigenetic aging in middle-aged adults, the impact during childhood remains unclear. We tested the association between body mass index (BMI) and accelerated epigenetic aging in a cohort of high-risk children. Participants were children (N = 273, aged 8 to 14 years, 82% investigated for maltreatment) recruited to the Child Health Study, an ongoing prospective study of youth investigated for maltreatment and a comparison youth. BMI was measured as a continuous variable. Accelerated epigenetic aging of blood leukocytes was defined as the age-adjusted residuals of several established epigenetic aging clocks (Horvath, Hannum, GrimAge, PhenoAge) along with a newer algorithm, the DunedinPoAm, developed to quantify the pace-of-aging. Hypotheses were tested with generalized linear models. Higher age-and sex- adjusted z-scored BMI was significantly correlated with household income, blood cell counts, and three of the accelerated epigenetic aging measures: GrimAge (r = 0.31, P < .0001), PhenoAge (r = 0.24, P < .0001), and DunedinPoAm (r = 0.38, P < .0001). In fully adjusted models, GrimAge (ß = 0.07; P = .0009) and DunedinPoAm (ß = 0.0017; P < .0001) remained significantly associated with higher age- and sex-adjusted z-scored BMI. Maltreatment-status was not associated with accelerated epigenetic aging. In a high-risk cohort of children, higher BMI predicted epigenetic aging as assessed by two epigenetic aging clocks. These results suggest the association between obesity and accelerated epigenetic aging begins in early life, with implications for future morbidity and mortality risk.

Semi-automated NMR Pipeline for Environmental Exposures: New Insights on the Metabolomics of Smokers versus Non-smokers.

Environmental exposure pathophysiology related to smoking can yield metabolic changes that are difficult to describe in a biologically informative fashion with manual proprietary software. Nuclear magnetic resonance (NMR) spectroscopy detects compounds found in biofluids yielding a metabolic snapshot. We applied our semi-automated NMR pipeline for a secondary analysis of a smoking study (MTBLS374 from the MetaboLights repository) (n = 112). This involved quality control (in the form of data preprocessing), automated metabolite quantification, and analysis. With our approach we putatively identified 79 metabolites that were previously unreported in the dataset. Quantified metabolites were used for metabolic pathway enrichment analysis that replicated 1 enriched pathway with the original study as well as 3 previously unreported pathways. Our pipeline generated a new random forest (RF) classifier between smoking classes that revealed several combinations of compounds. This study broadens our metabolomic understanding of smoking exposure by 1) notably increasing the number of quantified metabolites with our analytic pipeline, 2) suggesting smoking exposure may lead to heterogenous metabolic responses according to random forest modeling, and 3) modeling how newly quantified individual metabolites can determine smoking status. Our approach can be applied to other NMR studies to characterize environmental risk factors, allowing for the discovery of new biomarkers of disease and exposure status.

What about the environment? Leveraging multi-omic datasets to characterize the environment's role in human health.

The environment plays an important role in mediating human health. In this session we consider research addressing ways to overcome the challenges associated with studying the multifaceted and ever-changing environment. Environmental health research has a need for technological and methodological advances which will further our knowledge of how exposures precipitate complex phenotypes and exacerbate disease.