Genetic variant rs1205 is associated with COVID-19 outcomes: The Strong Heart Study and Strong Heart Family Study.

BACKGROUND: Although COVID-19 infection has been associated with a number of clinical and environmental risk factors, host genetic variation has also been associated with the incidence and morbidity of infection. The CRP gene codes for a critical component of the innate immune system and CRP variants have been reported associated with infectious disease and vaccination outcomes. We investigated possible associations between COVID-19 outcome and a limited number of candidate gene variants including rs1205. METHODOLOGY/PRINCIPAL FINDINGS: The Strong Heart and Strong Heart Family studies have accumulated detailed genetic, cardiovascular risk and event data in geographically dispersed American Indian communities since 1988. Genotypic data and 91 COVID-19 adjudicated deaths or hospitalizations from 2/1/20 through 3/1/23 were identified among 3,780 participants in two subsets. Among 21 candidate variants including genes in the interferon response pathway, APOE, TMPRSS2, TLR3, the HLA complex and the ABO blood group, only rs1205, a 3' untranslated region variant in the CRP gene, showed nominally significant association in T-dominant model analyses (odds ratio 1.859, 95%CI 1.001-3.453, p = 0.049) after adjustment for age, sex, center, body mass index, and a history of cardiovascular disease. Within the younger subset, association with the rs1205 T-Dom genotype was stronger, both in the same adjusted logistic model and in the SOLAR analysis also adjusting for other genetic relatedness. CONCLUSION: A T-dominant genotype of rs1205 in the CRP gene is associated with COVID-19 death or hospitalization, even after adjustment for relevant clinical factors and potential participant relatedness. Additional study of other populations and genetic variants of this gene are warranted.

Omics feature selection with the extended SIS R package: identification of a body mass index epigenetic multi-marker in the Strong Heart Study.

The statistical analysis of omics data poses a great computational challenge given its ultra-high dimensional nature and frequent between-features correlation. In this work, we extended the Iterative Sure Independence Screening (ISIS) algorithm by pairing ISIS with elastic-net (Enet) and two versions of adaptive Enet (AEnet and MSAEnet) to efficiently improve feature selection and effect estimation in omics research. We subsequently used genome-wide human blood DNA methylation data from American Indians of the Strong Heart Study (N=2,235 participants), measured in 1989-1991, to compare the performance (predictive accuracy, coefficient estimation and computational efficiency) of SIS-paired regularization methods to Bayesian shrinkage and traditional linear regression to identify epigenomic multi-marker of body mass index. ISIS-AEnet outperformed the other methods in prediction. In biological pathway enrichment analysis of genes annotated to BMI-related differentially methylated positions, ISIS-AEnet captured most of the enriched pathways in common for at least two of all the evaluated methods. ISIS-AEnet can favor biological discovery because it identifies the most robust biological pathways while achieving an optimal balance between bias and efficient feature selection. In the extended SIS R package, we also implemented ISIS paired with Cox and logistic regression for time-to-event and binary endpoints, respectively, and bootstrap confidence intervals for the estimated regression coefficients.

Blood DNA methylation and liver cancer in American Indians: evidence from the Strong Heart Study.

PURPOSE: Liver cancer incidence among American Indians/Alaska Natives has risen over the past 20 years. Peripheral blood DNA methylation may be associated with liver cancer and could be used as a biomarker for cancer risk. We evaluated the association of blood DNA methylation with risk of liver cancer. METHODS: We conducted a prospective cohort study in 2324 American Indians, between age 45 and 75 years, from Arizona, Oklahoma, North Dakota and South Dakota who participated in the Strong Heart Study between 1989 and 1991. Liver cancer deaths (n = 21) were ascertained using death certificates obtained through 2017. The mean follow-up duration (SD) for non-cases was 25.1 (5.6) years and for cases, 11.0 (8.8) years. DNA methylation was assessed from blood samples collected at baseline using MethylationEPIC BeadChip 850 K arrays. We used Cox regression models adjusted for age, sex, center, body mass index, low-density lipoprotein cholesterol, smoking, alcohol consumption, and immune cell proportions to examine the associations. RESULTS: We identified 9 CpG sites associated with liver cancer. cg16057201 annotated to MRFAP1) was hypermethylated among cases vs. non-cases (hazard ratio (HR) for one standard deviation increase in methylation was 1.25 (95% CI 1.14, 1.37). The other eight CpGs were hypomethylated and the corresponding HRs (95% CI) ranged from 0.58 (0.44, 0.75) for cg04967787 (annotated to PPRC1) to 0.77 (0.67, 0.88) for cg08550308. We also assessed 7 differentially methylated CpG sites associated with liver cancer in previous studies. The adjusted HR for cg15079934 (annotated to LPS1) was 1.93 (95% CI 1.10, 3.39). CONCLUSIONS: Blood DNA methylation may be associated with liver cancer mortality and may be altered during the development of liver cancer.

Dietary magnesium, C-reactive protein and interleukin-6: The Strong Heart Family Study.

OBJECTIVES: To examine the associations of dietary Mg intake with inflammatory biomarkers (C-reactive protein (CRP) and interleukin 6 (IL-6)), and the interaction of dietary Mg intake with single nucleotide polymorphism (SNP) rs3740393, a SNP related to Mg metabolism and transport, on CRP and IL-6 among American Indians (AIs). METHODS: This cross-sectional study included AI participants (n = 1,924) from the Strong Heart Family Study (SHFS). Mg intake from foods and dietary supplements was ascertained using a 119-item Block food frequency questionnaire, CRP and IL-6 were measured from blood, and SNP rs3740393 was genotyped using MetaboChip. Generalized estimating equations were used to examine associations of Mg intake, and the interaction between rs3740393 and dietary Mg, with CRP and IL-6. RESULTS: Reported Mg intake was not associated with CRP or IL-6, irrespective of genotype. A significant interaction (p-interaction = 0.018) was observed between Mg intake and rs3740393 on IL-6. Among participants with the C/C genotype, for every 1 SD higher in log-Mg, log-IL-6 was 0.04 (95% CI: -0.10 to 0.17) pg/mL higher. Among participants with the C/G genotype, for every 1 SD higher in log-Mg, log-IL-6 was 0.08 (95% CI: -0.21 to 0.05) pg/mL lower, and among participants with the G/G genotype, for every 1 SD higher in log-Mg, log-IL-6 was 0.19 (95% CI: -0.38 to -0.01) pg/mL lower. CONCLUSIONS: Mg intake may be associated with lower IL-6 with increasing dosage of the G allele at rs3740393. Future research is necessary to replicate this finding and examine other Mg-related genes that influence associations of Mg intake with inflammation.

Arsenic Exposure and Epigenetic Aging: The Association with Cardiovascular Disease and All-Cause Mortality in the Strong Heart Study.

BACKGROUND: Inorganic arsenic (As) may increase the risk of cardiovascular disease (CVD) and all-cause mortality through accelerated aging, which can be estimated using epigenetic-based measures. OBJECTIVES: We evaluated three DNA methylation-based aging measures (PhenoAge, GrimAge, DunedinPACE) (epigenetic aging measures) as potential mediators of the previously reported association of As exposure with CVD incidence, CVD mortality, and all-cause mortality in the Strong Heart Study (SHS), an epidemiological cohort of American Indian adults. METHODS: Blood DNA methylation and urinary As levels were measured in 2,323 SHS participants (41.5% men, mean age of 55 years old). PhenoAge and GrimAge values were calculated using a residual-based method. We tested the association of urinary As with epigenetic aging measures using linear regression, the association of epigenetic aging measures with the three health outcomes using additive hazards models, and the mediation of As-related CVD incidence, CVD mortality, and all-cause mortality by epigenetic aging measures using the product of coefficients method. RESULTS: SHS participants with higher vs. lower urinary As levels had similar PhenoAge age, older GrimAge age, and faster DunedinPACE. An interquartile range increase in urinary As was associated with higher of PhenoAge age acceleration [mean difference (95% confidence interval)=0.48 (0.17, 0.80) years], GrimAge age acceleration [0.80 (0.60, 1.00) years], and DunedinPACE [0.011 (0.005, 0.018)], after adjusting for age, sex, center location, genetic components, smoking status, and body mass index. Of the 347 incident CVD events per 100,000 person-years associated with a doubling in As exposure, 21.3% (9.1, 57.1) and 22.6% (9.5, 56.9), were attributable to differences in GrimAge and DunedinPACE, respectively. DISCUSSION: Arsenic exposure was associated with older GrimAge and faster DunedinPACE measures of biological age. Furthermore, accelerated biological aging measured from DNA methylation accounted for a relevant fraction of As-associated risk for CVD, CVD mortality, and all-cause mortality in the SHS, supporting the role of As in accelerated aging. Research of the biological underpinnings can contribute to a better understanding of the role of aging in arsenic-related disease. https://doi.org/10.1289/EHP11981.

A Scan of Pleiotropic Immune Mediated Disease Genes Identifies Novel Determinants of Baseline FVIII Inhibitor Status in Hemophilia-A.

Hemophilia-A (HA) is caused by heterogeneous loss-of-function factor (F)VIII gene (F8)-mutations and deficiencies in plasma-FVIII-activity that impair intrinsic-pathway-mediated coagulation-amplification. The standard-of-care for severe-HA-patients is regular infusions of therapeutic-FVIII-proteins (tFVIIIs) but ~30% develop neutralizing-tFVIII-antibodies called "FVIII-inhibitors (FEIs)" and become refractory. We used the PATH study and ImmunoChip to scan immune-mediated-disease (IMD)-genes for novel and/or replicated genomic-sequence-variations associated with baseline-FEI-status while accounting for non-independence of data due to genetic-relatedness and F8-mutational-heterogeneity. The baseline-FEI-status of 450 North American PATH subjects-206 with black-African-ancestry and 244 with white-European-ancestry-was the dependent variable. The F8-mutation-data and a genetic-relatedness matrix were incorporated into a binary linear-mixed model of genetic association with baseline-FEI-status. We adopted a gene-centric-association-strategy to scan, as candidates, pleiotropic-IMD-genes implicated in the development of either ³2 autoimmune-/autoinflammatory-disorders (AADs) or ³1 AAD and FEIs. Baseline-FEI-status was significantly associated with SNPs assigned to NOS2A (rs117382854; p=3.2E-6) and B3GNT2 (rs10176009; p=5.1E-6), which have functions in anti-microbial-/-tumoral-immunity. Among IMD-genes implicated in FEI-risk previously, we identified strong associations with CTLA4 assigned SNPs (p=2.2E-5). The F8-mutation-effect underlies ~15% of the total heritability for baseline-FEI-status. Additive genetic heritability and SNPs in IMD-genes account for >50% of the patient-specific variability in baseline-FEI-status. Race is a significant determinant independent of F8-mutation-effects and non-F8-genetics.

Smoking, blood DNA methylation sites and lung cancer risk.

Altered DNA methylation (DNAm) might be a biological intermediary in the pathway from smoking to lung cancer. In this study, we investigated the contribution of differential blood DNAm to explain the association between smoking and lung cancer incidence. Blood DNAm was measured in 2321 Strong Heart Study (SHS) participants. Incident lung cancer was assessed as time to event diagnoses. We conducted mediation analysis, including validation with DNAm and paired gene expression data from the Framingham Heart Study (FHS). In the SHS, current versus never smoking and pack-years single-mediator models showed, respectively, 29 and 21 differentially methylated positions (DMPs) for lung cancer with statistically significant mediated effects (14 of 20 available, and five of 14 available, positions, replicated, respectively, in FHS). In FHS, replicated DMPs showed gene expression downregulation largely in trans, and were related to biological pathways in cancer. The multimediator model identified that DMPs annotated to the genes AHRR and IER3 jointly explained a substantial proportion of lung cancer. Thus, the association of smoking with lung cancer was partly explained by differences in baseline blood DNAm at few relevant sites. Experimental studies are needed to confirm the biological role of identified eQTMs and to evaluate potential implications for early detection and control of lung cancer.

Metal mixtures and DNA methylation measures of biological aging in American Indian populations.

INTRODUCTION: Native American communities suffer disproportionately from elevated metal exposures and increased risk for cardiovascular diseases and diabetes. DNA methylation is a sensitive biomarker of aging-related processes and novel epigenetic-based "clocks" can be used to estimate accelerated biological aging that may underlie increased risk. Metals alter DNA methylation, yet little is known about their individual and combined impact on epigenetic age acceleration. Our objective was to investigate the associations of metals on several DNA methylation-based aging measures in the Strong Heart Study (SHS) cohort. METHODS: Blood DNA methylation data from 2,301 SHS participants was used to calculate age acceleration of epigenetic clocks (PhenoAge, GrimAge, DunedinPACE, Hannum, Horvath). Urinary metals [arsenic (As), cadmium (Cd), tungsten (W), zinc (Zn), selenium (Se), molybdenum (Mo)] were creatinine-adjusted and categorized into quartiles. We examined associations of individual metals through linear regression models and used Bayesian Kernel Machine Regression (BKMR) for the impact of the total metal mixture on epigenetic age acceleration. RESULTS: The mixture of nonessential metals (W, As, Cd) was associated with greater GrimAge acceleration and DunedinPACE, while the essential metal mixture (Se, Zn, Mo) was associated with lower epigenetic age acceleration. Cd was associated with increased epigenetic age acceleration across all clocks and BKMR analysis suggested nonlinear associations between Se and DunedinPACE, GrimAge, and PhenoAge acceleration. No interactions between individual metals were observed. The associations between Cd, Zn, and epigenetic age acceleration were greater in never smokers in comparison to current/former smokers. CONCLUSION: Nonessential metals were positively associated with greater epigenetic age acceleration, with strongest associations observed between Cd and DunedinPACE and GrimAge acceleration. In contrast, essential metals were associated with lower epigenetic aging. Examining the influence of metal mixtures on epigenetic age acceleration can provide insight into metals and aging-related diseases.

Maternal DNA methylation signatures of arsenic exposure is associated with adult offspring insulin resistance in the Strong Heart Study.

Exposure to low to moderate arsenic (As) levels has been associated with type 2 diabetes (T2D) and other chronic diseases in American Indian communities. Prenatal exposure to As may also increase the risk for T2D in adulthood, and maternal As has been associated with adult offspring metabolic health measurements. We hypothesized that T2D-related outcomes in adult offspring born to women exposed to low to moderate As can be evaluated utilizing a maternally-derived molecular biosignature of As exposure. Herein, we evaluated the association of maternal DNA methylation with incident T2D and insulin resistance (Homeostatic model assessment of insulin resistance [HOMA2-IR]) in adult offspring. For DNA methylation, we used 20 differentially methylated cytosine-guanine dinucleotides (CpG) previously associated with the sum of inorganic and methylated As species (ΣAs) in urine in the Strong Heart Study (SHS). Of these 20 CpGs, we found six CpGs nominally associated (p < 0.05) with HOMA2-IR in a fully adjusted model that included clinically relevant covariates and offspring adiposity measurements; a similar model that adjusted instead for maternal adiposity measurements found three CpGs nominally associated with HOMA2-IR, two of which overlapped the offspring adiposity model. After adjusting for multiple comparisons, cg03036214 remained associated with HOMA2-IR (q < 0.10) in the offspring adiposity model. The odds ratio of incident T2D increased with an increase in maternal DNA methylation at one HOMA2-IR associated CpG in the model adjusting for offspring adiposity, cg12116137, whereas adjusting for maternal adiposity had a minimal effect on the association. Our data suggests offspring adiposity, rather than maternal adiposity, potentially influences the effects of maternal DNAm signatures on offspring metabolic health parameters. Here, we have presented evidence supporting a role for epigenetic biosignatures of maternal As exposure as a potential biomarker for evaluating risk of T2D-related outcomes in offspring later in life.

Sequencing-based fine-mapping and in silico functional characterization of the 10q24.32 arsenic metabolism efficiency locus across multiple arsenic-exposed populations.

Inorganic arsenic is highly toxic and carcinogenic to humans. Exposed individuals vary in their ability to metabolize arsenic, and variability in arsenic metabolism efficiency (AME) is associated with risks of arsenic-related toxicities. Inherited genetic variation in the 10q24.32 region, near the arsenic methyltransferase (AS3MT) gene, is associated with urine-based measures of AME in multiple arsenic-exposed populations. To identify potential causal variants in this region, we applied fine mapping approaches to targeted sequencing data generated for exposed individuals from Bangladeshi, American Indian, and European American populations (n = 2,357, 557, and 648 respectively). We identified three independent association signals for Bangladeshis, two for American Indians, and one for European Americans. The size of the confidence sets for each signal varied from 4 to 85 variants. There was one signal shared across all three populations, represented by the same SNP in American Indians and European Americans (rs191177668) and in strong linkage disequilibrium (LD) with a lead SNP in Bangladesh (rs145537350). Beyond this shared signal, differences in LD patterns, minor allele frequency (MAF) (e.g., rs12573221 ~13% in Bangladesh ~0.2% among American Indians), and/or heterogeneity in effect sizes across populations likely contributed to the apparent population specificity of the additional identified signals. One of our potential causal variants influences AS3MT expression and nearby DNA methylation in numerous GTEx tissue types (with rs4919690 as a likely causal variant). Several SNPs in our confidence sets overlap transcription factor binding sites and cis-regulatory elements (from ENCODE). Taken together, our analyses reveal multiple potential causal variants in the 10q24.32 region influencing AME, including a variant shared across populations, and elucidate potential biological mechanisms underlying the impact of genetic variation on AME.

DNA methylation analysis is used to identify novel genetic loci associated with circulating fibrinogen levels in blood.

BACKGROUND: Fibrinogen plays an essential role in blood coagulation and inflammation. Circulating fibrinogen levels may be determined based on interindividual differences in DNA methylation at cytosine-phosphate-guanine (CpG) sites and vice versa. OBJECTIVES: To perform an EWAS to examine an association between blood DNA methylation levels and circulating fibrinogen levels to better understand its biological and pathophysiological actions. METHODS: We performed an epigenome-wide association study of circulating fibrinogen levels in 18 037 White, Black, American Indian, and Hispanic participants, representing 14 studies from the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium. Circulating leukocyte DNA methylation was measured using the Illumina 450K array in 12 904 participants and using the EPIC array in 5133 participants. In each study, an epigenome-wide association study of fibrinogen was performed using linear mixed models adjusted for potential confounders. Study-specific results were combined using array-specific meta-analysis, followed by cross-replication of epigenome-wide significant associations. We compared models with and without CRP adjustment to examine the role of inflammation. RESULTS: We identified 208 and 87 significant CpG sites associated with fibrinogen levels from the 450K (p < 1.03 × 10-7) and EPIC arrays (p < 5.78 × 10-8), respectively. There were 78 associations from the 450K array that replicated in the EPIC array and 26 vice versa. After accounting for overlapping sites, there were 83 replicated CpG sites located in 61 loci, of which only 4 have been previously reported for fibrinogen. The examples of genes located near these CpG sites were SOCS3 and AIM2, which are involved in inflammatory pathways. The associations of all 83 replicated CpG sites were attenuated after CRP adjustment, although many remained significant. CONCLUSION: We identified 83 CpG sites associated with circulating fibrinogen levels. These associations are partially driven by inflammatory pathways shared by both fibrinogen and CRP.

Epigenetics of type 2 diabetes and diabetes-related outcomes in the Strong Heart Study.

BACKGROUND: The prevalence of type 2 diabetes has dramatically increased in the past years. Increasing evidence supports that blood DNA methylation, the best studied epigenetic mark, is related to diabetes risk. Few prospective studies, however, are available. We studied the association of blood DNA methylation with diabetes in the Strong Heart Study. We used limma, Iterative Sure Independence Screening and Cox regression to study the association of blood DNA methylation with fasting glucose, HOMA-IR and incident type 2 diabetes among 1312 American Indians from the Strong Heart Study. DNA methylation was measured using Illumina's MethylationEPIC beadchip. We also assessed the biological relevance of our findings using bioinformatics analyses. RESULTS: Among the 358 differentially methylated positions (DMPs) that were cross-sectionally associated either with fasting glucose or HOMA-IR, 49 were prospectively associated with incident type 2 diabetes, although no DMPs remained significant after multiple comparisons correction. Multiple of the top DMPs were annotated to genes with relevant functions for diabetes including SREBF1, associated with obesity, type 2 diabetes and insulin sensitivity; ABCG1, involved in cholesterol and phospholipids transport; and HDAC1, of the HDAC family. (HDAC inhibitors have been proposed as an emerging treatment for diabetes and its complications.) CONCLUSIONS: Our results suggest that differences in peripheral blood DNA methylation are related to cross-sectional markers of glucose metabolism and insulin activity. While some of these DMPs were modestly associated with prospective incident type 2 diabetes, they did not survive multiple testing. Common DMPs with diabetes epigenome-wide association studies from other populations suggest a partially common epigenomic signature of glucose and insulin activity.

Epigenetic Biomarkers of Lead Exposure and Cardiovascular Disease: Prospective Evidence in the Strong Heart Study.

Background Lead is a cardiotoxic metal with a variety of adverse health effects. In the absence of data on bone lead exposure, epigenetic biomarkers can serve as indicators of cumulative lead exposure and body burden. Herein, we leveraged novel epigenetic biomarkers of lead exposure to investigate their association with cardiovascular disease (CVD) incidence and mortality. Methods and Results Blood DNA methylation was measured using the Illumina MethylationEPIC BeadChip among 2231 participants of the Strong Heart Study (SHS) at baseline (1989-1991). Epigenetic biomarkers of lead levels in blood, patella, and tibia were estimated using previously identified cytosine-guanine dinucleotide (CpG) sites. CVD incidence and mortality data were available through 2017. Median concentrations of lead epigenetic biomarkers were 13.8 µg/g, 21.3 µg/g, and 2.9 µg/dL in tibia, patella, and blood, respectively. In adjusted models, the hazard ratio (HR) (95% CI) of CVD mortality per doubling increase in lead epigenetic biomarkers were 1.42 (1.07-1.87) for tibia lead, 1.22 (0.93-1.60) for patella lead, and 1.57 (1.16-2.11) for blood lead. The corresponding HRs for incident CVD were 0.99 (0.83-1.19), 1.07 (0.89-1.29), and 1.06 (0.87-1.30). The association between the tibia lead epigenetic biomarker and CVD mortality was modified by sex (interaction P value: 0.014), with men at increased risk (HR, 1.42 [95% CI, 1.17-1.72]) compared with women (HR, 1.04 [95% CI, 0.89-1.22]). Conclusions Tibia and blood epigenetic biomarkers were associated with increased risk of CVD mortality, potentially reflecting the cardiovascular impact of cumulative and recent lead exposures. These findings support that epigenetic biomarkers of lead exposure may capture some of the disease risk associated with lead exposure.

Association of protein function-altering variants with cardiometabolic traits: the strong heart study.

Clinical and biomarker phenotypic associations for carriers of protein function-altering variants may help to elucidate gene function and health effects in populations. We genotyped 1127 Strong Heart Family Study participants for protein function-altering single nucleotide variants (SNV) and indels selected from a low coverage whole exome sequencing of American Indians. We tested the association of each SNV/indel with 35 cardiometabolic traits. Among 1206 variants (average minor allele count = 20, range of 1 to 1064), ~ 43% were not present in publicly available repositories. We identified seven SNV-trait significant associations including a missense SNV at ABCA10 (rs779392624, p = 8 × 10-9) associated with fasting triglycerides, which gene product is involved in macrophage lipid homeostasis. Among non-diabetic individuals, missense SNVs at four genes were associated with fasting insulin adjusted for BMI (PHIL, chr6:79,650,711, p = 2.1 × 10-6; TRPM3, rs760461668, p = 5 × 10-8; SPTY2D1, rs756851199, p = 1.6 × 10-8; and TSPO, rs566547284, p = 2.4 × 10-6). PHIL encoded protein is involved in pancreatic ß-cell proliferation and survival, and TRPM3 protein mediates calcium signaling in pancreatic ß-cells in response to glucose. A genetic risk score combining increasing insulin risk alleles of these four genes was associated with 53% (95% confidence interval 1.09, 2.15) increased odds of incident diabetes and 83% (95% confidence interval 1.35, 2.48) increased odds of impaired fasting glucose at follow-up. Our study uncovered novel gene-trait associations through the study of protein-coding variants and demonstrates the advantages of association screenings targeting diverse and high-risk populations to study variants absent in publicly available repositories.

Arsenic Exposure, Blood DNA Methylation, and Cardiovascular Disease.

BACKGROUND: Epigenetic dysregulation has been proposed as a key mechanism for arsenic-related cardiovascular disease (CVD). We evaluated differentially methylated positions (DMPs) as potential mediators on the association between arsenic and CVD. METHODS: Blood DNA methylation was measured in 2321 participants (mean age 56.2, 58.6% women) of the Strong Heart Study, a prospective cohort of American Indians. Urinary arsenic species were measured using high-performance liquid chromatography coupled to inductively coupled plasma mass spectrometry. We identified DMPs that are potential mediators between arsenic and CVD. In a cross-species analysis, we compared those DMPs with differential liver DNA methylation following early-life arsenic exposure in the apoE knockout (apoE-/-) mouse model of atherosclerosis. RESULTS: A total of 20 and 13 DMPs were potential mediators for CVD incidence and mortality, respectively, several of them annotated to genes related to diabetes. Eleven of these DMPs were similarly associated with incident CVD in 3 diverse prospective cohorts (Framingham Heart Study, Women's Health Initiative, and Multi-Ethnic Study of Atherosclerosis). In the mouse model, differentially methylated regions in 20 of those genes and DMPs in 10 genes were associated with arsenic. CONCLUSIONS: Differential DNA methylation might be part of the biological link between arsenic and CVD. The gene functions suggest that diabetes might represent a relevant mechanism for arsenic-related cardiovascular risk in populations with a high burden of diabetes.

An epigenome-wide study of DNA methylation profiles and lung function among American Indians in the Strong Heart Study.

BACKGROUND: Epigenetic modifications, including DNA methylation (DNAm), are often related to environmental exposures, and are increasingly recognized as key processes in the pathogenesis of chronic lung disease. American Indian communities have a high burden of lung disease compared to the national average. The objective of this study was to investigate the association of DNAm and lung function in the Strong Heart Study (SHS). We conducted a cross-sectional study of American Indian adults, 45-74 years of age who participated in the SHS. DNAm was measured using the Illumina Infinium Human MethylationEPIC platform at baseline (1989-1991). Lung function was measured via spirometry, including forced expiratory volume in 1 s (FEV1) and forced vital capacity (FVC), at visit 2 (1993-1995). Airflow limitation was defined as FEV1 < 70% predicted and FEV1/FVC < 0.7, restriction was defined as FEV1/FVC > 0.7 and FVC < 80% predicted, and normal spirometry was defined as FEV1/FVC > 0.7, FEV1 > 70% predicted, FVC > 80% predicted. We used elastic-net models to select relevant CpGs for lung function and spirometry-defined lung disease. We also conducted bioinformatic analyses to evaluate the biological plausibility of the findings. RESULTS: Among 1677 participants, 21.2% had spirometry-defined airflow limitation and 13.6% had spirometry-defined restrictive pattern lung function. Elastic-net models selected 1118 Differentially Methylated Positions (DMPs) as predictors of airflow limitation and 1385 for restrictive pattern lung function. A total of 12 DMPs overlapped between airflow limitation and restrictive pattern. EGFR, MAPK1 and PRPF8 genes were the most connected nodes in the protein-protein interaction network. Many of the DMPs targeted genes with biological roles related to lung function such as protein kinases. CONCLUSION: We found multiple differentially methylated CpG sites associated with chronic lung disease. These signals could contribute to better understand molecular mechanisms involved in lung disease, as assessed systemically, as well as to identify patterns that could be useful for diagnostic purposes. Further experimental and longitudinal studies are needed to assess whether DNA methylation has a causal role in lung disease.

Genetic variants in ALDH1L1 and GLDC influence the serine-to-glycine ratio in Hispanic children.

BACKGROUND: Glycine is a proteogenic amino acid that is required for numerous metabolic pathways, including purine, creatine, heme, and glutathione biosynthesis. Glycine formation from serine, catalyzed by serine hydroxy methyltransferase, is the major source of this amino acid in humans. Our previous studies in a mouse model have shown a crucial role for the 10-formyltetrahydrofolate dehydrogenase enzyme in serine-to-glycine conversion. OBJECTIVES: We sought to determine the genomic influence on the serine-glycine ratio in 803 Hispanic children from 319 families of the Viva La Familia cohort. METHODS: We performed a genome-wide association analysis for plasma serine, glycine, and the serine-glycine ratio in Sequential Oligogenic Linkage Analysis Routines while accounting for relationships among family members. RESULTS: All 3 parameters were significantly heritable (h2 = 0.22-0.78; P < 0.004). The strongest associations for the serine-glycine ratio were with single nucleotide polymorphisms (SNPs) in aldehyde dehydrogenase 1 family member L1 (ALDH1L1) and glycine decarboxylase (GLDC) and for glycine with GLDC (P < 3.5 × 10-8; effect sizes, 0.03-0.07). No significant associations were found for serine. We also conducted a targeted genetic analysis with ALDH1L1 exonic SNPs and found significant associations between the serine-glycine ratio and rs2886059 (ß = 0.68; SE, 0.25; P = 0.006) and rs3796191 (ß = 0.25; SE, 0.08; P = 0.003) and between glycine and rs3796191 (ß = -0.08; SE, 0.02; P = 0.0004). These exonic SNPs were further associated with metabolic disease risk factors, mainly adiposity measures (P < 0.006). Significant genetic and phenotypic correlations were found for glycine and the serine-glycine ratio with metabolic disease risk factors, including adiposity, insulin sensitivity, and inflammation-related phenotypes [estimate of genetic correlation = -0.37 to 0.35 (P < 0.03); estimate of phenotypic correlation = -0.19 to 0.13 (P < 0.006)]. The significant genetic correlations indicate shared genetic effects among glycine, the serine-glycine ratio, and adiposity and insulin sensitivity phenotypes. CONCLUSIONS: Our study suggests that ALDH1L1 and GLDC SNPs influence the serine-to-glycine ratio and metabolic disease risk.

Replication of Integrative Data Analysis for Adipose Tissue Dysfunction, Low-Grade Inflammation, Postprandial Responses and OMICs Signatures in Symptom-Free Adults.

We previously reported preliminary characterization of adipose tissue (AT) dysfunction through the adiponectin/leptin ratio (ALR) and fasting/postprandial (F/P) gene expression in subcutaneous (SQ) adipose tissue (AT) biopsies obtained from participants in the GEMM study, a precision medicine research project. Here we present integrative data replication of previous findings from an increased number of GEMM symptom-free (SF) adults (N = 124) to improve characterization of early biomarkers for cardiovascular (CV)/immunometabolic risk in SF adults with AT dysfunction. We achieved this goal by taking advantage of the rich set of GEMM F/P 5 h time course data and three tissue samples collected at the same time and frequency on each adult participant (F/P blood, biopsies of SQAT and skeletal muscle (SKM)). We classified them with the presence/absence of AT dysfunction: low (<1) or high (>1) ALR. We also examined the presence of metabolically healthy (MH)/unhealthy (MUH) individuals through low-grade chronic subclinical inflammation (high sensitivity C-reactive protein (hsCRP)), whole body insulin sensitivity (Matsuda Index) and Metabolic Syndrome criteria in people with/without AT dysfunction. Molecular data directly measured from three tissues in a subset of participants allowed fine-scale multi-OMIC profiling of individual postprandial responses (RNA-seq in SKM and SQAT, miRNA from plasma exosomes and shotgun lipidomics in blood). Dynamic postprandial immunometabolic molecular endophenotypes were obtained to move towards a personalized, patient-defined medicine. This study offers an example of integrative translational research, which applies bench-to-bedside research to clinical medicine. Our F/P study design has the potential to characterize CV/immunometabolic early risk detection in support of precision medicine and discovery in SF individuals.

Blood DNA Methylation and Incident Coronary Heart Disease: Evidence From the Strong Heart Study.

Importance: American Indian communities experience a high burden of coronary heart disease (CHD). Strategies are needed to identify individuals at risk and implement preventive interventions. Objective: To investigate the association of blood DNA methylation (DNAm) with incident CHD using a large number of methylation sites (cytosine-phosphate-guanine [CpG]) in a single model. Design, Setting, and Participants: This prospective study, including a discovery cohort (the Strong Heart Study [SHS]) and 4 additional cohorts (the Women's Health Initiative [WHI], the Framingham Heart Study [FHS], the Atherosclerosis Risk in Communities Study ([ARIC]-Black, and ARIC-White), evaluated 12 American Indian communities in 4 US states; African American women, Hispanic women, and White women throughout the US; White men and White women from Massachusetts; and Black men and women and White men and women from 4 US communities. A total of 2321 men and women (mean [SD] follow-up, 19.1 [9.2] years) were included in the SHS, 1874 women (mean [SD] follow-up, 15.8 [5.9] years) in the WHI, 2128 men and women (mean [SD] follow-up, 7.7 [1.8] years) in the FHS, 2114 men and women (mean [SD] follow-up, 20.9 [7.2] years) in the ARIC-Black, and 931 men and women (mean [SD] follow-up, 20.9 [7.2] years) in the ARIC-White. Data were collected from May 1989 to December 2018 and analyzed from February 2019 to May 2021. Exposure: Blood DNA methylation. Main Outcome and Measure: Using a high-dimensional time-to-event elastic-net model for the association of 407 224 CpG sites with incident CHD in the SHS (749 events), this study selected the differentially methylated CpG positions (DMPs) selected in the SHS and evaluated them in the WHI (531 events), FHS (143 events), ARIC-Black (350 events), and ARIC-White (121 events) cohorts. Results: The median (IQR) age of participants in SHS was 55 (49-62) years, and 1359 participants (58.6%) were women. Elastic-net models selected 505 DMPs associated with incident CHD in the SHS beyond established risk factors, center, blood cell counts, and genetic principal components. Among those DMPs, 33 were commonly selected in 3 or 4 of the other cohorts and the pooled hazard ratios from the standard Cox models were significant at P < .05 for 10 of the DMPs. For example, the hazard ratio (95% CI) for CHD comparing the 90th and 10th percentiles of differentially methylated CpGs was 0.86 (0.78-0.95) for cg16604233 (tagged to COL11A2) and 1.23 (1.08-1.39) for cg09926486 (tagged to FRMD5). Some of the DMPs were consistent in the direction of the association; others showed associations in opposite directions across cohorts. Untargeted independent elastic-net models of CHD showed distinct DMPs, genes, and network of genes in the 5 cohorts. Conclusions and Relevance: In this multi-cohort study, blood-based DNAm findings supported an association between a complex blood epigenomic signature and CHD that was largely different across populations.