Rare variant contribution to the heritability of coronary artery disease.

Whole genome sequences (WGS) enable discovery of rare variants which may contribute to missing heritability of coronary artery disease (CAD). To measure their contribution, we apply the GREML-LDMS-I approach to WGS of 4949 cases and 17,494 controls of European ancestry from the NHLBI TOPMed program. We estimate CAD heritability at 34.3% assuming a prevalence of 8.2%. Ultra-rare (minor allele frequency ≤ 0.1%) variants with low linkage disequilibrium (LD) score contribute ~50% of the heritability. We also investigate CAD heritability enrichment using a diverse set of functional annotations: i) constraint; ii) predicted protein-altering impact; iii) cis-regulatory elements from a cell-specific chromatin atlas of the human coronary; and iv) annotation principal components representing a wide range of functional processes. We observe marked enrichment of CAD heritability for most functional annotations. These results reveal the predominant role of ultra-rare variants in low LD on the heritability of CAD. Moreover, they highlight several functional processes including cell type-specific regulatory mechanisms as key drivers of CAD genetic risk.

The Genetic Determinants and Genomic Consequences of Non-Leukemogenic Somatic Point Mutations.

Clonal hematopoiesis (CH) is defined by the expansion of a lineage of genetically identical cells in blood. Genetic lesions that confer a fitness advantage, such as point mutations or mosaic chromosomal alterations (mCAs) in genes associated with hematologic malignancy, are frequent mediators of CH. However, recent analyses of both single cell-derived colonies of hematopoietic cells and population sequencing cohorts have revealed CH frequently occurs in the absence of known driver genetic lesions. To characterize CH without known driver genetic lesions, we used 51,399 deeply sequenced whole genomes from the NHLBI TOPMed sequencing initiative to perform simultaneous germline and somatic mutation analyses among individuals without leukemogenic point mutations (LPM), which we term CH-LPMneg. We quantified CH by estimating the total mutation burden. Because estimating somatic mutation burden without a paired-tissue sample is challenging, we developed a novel statistical method, the Genomic and Epigenomic informed Mutation (GEM) rate, that uses external genomic and epigenomic data sources to distinguish artifactual signals from true somatic mutations. We performed a genome-wide association study of GEM to discover the germline determinants of CH-LPMneg. After fine-mapping and variant-to-gene analyses, we identified seven genes associated with CH-LPMneg (TCL1A, TERT, SMC4, NRIP1, PRDM16, MSRA, SCARB1), and one locus associated with a sex-associated mutation pathway (SRGAP2C). We performed a secondary analysis excluding individuals with mCAs, finding that the genetic architecture was largely unaffected by their inclusion. Functional analyses of SMC4 and NRIP1 implicated altered HSC self-renewal and proliferation as the primary mediator of mutation burden in blood. We then performed comprehensive multi-tissue transcriptomic analyses, finding that the expression levels of 404 genes are associated with GEM. Finally, we performed phenotypic association meta-analyses across four cohorts, finding that GEM is associated with increased white blood cell count and increased risk for incident peripheral artery disease, but is not significantly associated with incident stroke or coronary disease events. Overall, we develop GEM for quantifying mutation burden from WGS without a paired-tissue sample and use GEM to discover the genetic, genomic, and phenotypic correlates of CH-LPMneg.

Region-based analysis with functional annotation identifies genes associated with cognitive function in South Asians from India.

The prevalence of dementia among South Asians across India is approximately 7.4% in those 60 years and older, yet little is known about genetic risk factors for dementia in this population. Most known risk loci for Alzheimer's disease (AD) have been identified from studies conducted in European Ancestry (EA) but are unknown in South Asians. Using whole-genome sequence data from 2680 participants from the Diagnostic Assessment of Dementia for the Longitudinal Aging Study of India (LASI-DAD), we performed a gene-based analysis of 84 genes previously associated with AD in EA. We investigated associations with the Hindi Mental State Examination (HMSE) score and factor scores for general cognitive function and five cognitive domains. For each gene, we examined missense/loss-of-function (LoF) variants and brain-specific promoter/enhancer variants, separately, both with and without incorporating additional annotation weights (e.g., deleteriousness, conservation scores) using the variant-Set Test for Association using Annotation infoRmation (STAAR). In the missense/LoF analysis without annotation weights and controlling for age, sex, state/territory, and genetic ancestry, three genes had an association with at least one measure of cognitive function (FDR q<0.1). APOE was associated with four measures of cognitive function, PICALM was associated with HMSE score, and TSPOAP1 was associated with executive function. The most strongly associated variants in each gene were rs429358 (APOE ε4), rs779406084 (PICALM), and rs9913145 (TSPOAP1). rs779406084 is a rare missense mutation that is more prevalent in LASI-DAD than in EA (minor allele frequency=0.075% vs. 0.0015%); the other two are common variants. No genes in the brain-specific promoter/enhancer analysis met criteria for significance. Results with and without annotation weights were similar. Missense/LoF variants in some genes previously associated with AD in EA are associated with measures of cognitive function in South Asians from India. Analyzing genome sequence data allows identification of potential novel causal variants enriched in South Asians.

Identification of novel neuroprotectants against vincristine-induced neurotoxicity in iPSC-derived neurons.

Chemotherapy-induced peripheral neuropathy (CIPN) is a disabling side effect of cancer chemotherapy that can often limit treatment options for cancer patients or have life-long neurodegenerative consequences that reduce the patient's quality of life. CIPN is caused by the detrimental actions of various chemotherapeutic agents on peripheral axons. Currently, there are no approved preventative measures or treatment options for CIPN, highlighting the need for the discovery of novel therapeutics and improving our understanding of disease mechanisms. In this study, we utilized human-induced pluripotent stem cell (hiPSC)-derived motor neurons as a platform to mimic axonal damage after treatment with vincristine, a chemotherapeutic used for the treatment of breast cancers, osteosarcomas, and leukemia. We screened a total of 1902 small molecules for neuroprotective properties in rescuing vincristine-induced axon growth deficits. From our primary screen, we identified 38 hit compounds that were subjected to secondary dose response screens. Six compounds showed favorable pharmacological profiles - AZD7762, A-674563, Blebbistatin, Glesatinib, KW-2449, and Pelitinib, all novel neuroprotectants against vincristine toxicity to neurons. In addition, four of these six compounds also showed efficacy against vincristine-induced growth arrest in human iPSC-derived sensory neurons. In this study, we utilized high-throughput screening of a large library of compounds in a therapeutically relevant assay. We identified several novel compounds that are efficacious in protecting different neuronal subtypes from the toxicity induced by a common chemotherapeutic agent, vincristine which could have therapeutic potential in the clinic.

Identification of novel neuroprotectants against vincristine-induced neurotoxicity in iPSC-derived neurons.

Chemotherapy-induced peripheral neuropathy (CIPN) is a disabling side effect of cancer chemotherapy that can often limit treatment options for cancer patients or have life-long neurodegenerative consequences that reduce the patient's quality of life. CIPN is caused by the detrimental actions of various chemotherapeutic agents on peripheral axons. Currently, there are no approved preventative measures or treatment options for CIPN, highlighting the need for the discovery of novel therapeutics and improving our understanding of disease mechanisms. In this study, we utilized human-induced pluripotent stem cell (hiPSC)-derived motor neurons as a platform to mimic axonal damage after treatment with vincristine, a chemotherapeutic used for the treatment of breast cancers, osteosarcomas, and leukemia. We screened a total of 1902 small molecules for neuroprotective properties in rescuing vincristine-induced axon growth deficits. From our primary screen, we identified 38 hit compounds that were subjected to secondary dose response screens. Six compounds showed favorable pharmacological profiles - AZD7762, A-674563, Blebbistatin, Glesatinib, KW-2449, and Pelitinib, all novel neuroprotectants against vincristine toxicity to neurons. In addition, four of these six compounds also showed efficacy against vincristine-induced growth arrest in human iPSC-derived sensory neurons. In this study, we utilized high-throughput screening of a large library of compounds in a therapeutically relevant assay. We identified several novel compounds that are efficacious in protecting different neuronal subtypes from the toxicity induced by a common chemotherapeutic agent, vincristine which could have therapeutic potential in the clinic.

From hit to vial: Precision discovery and development of an imidazopyrimidine TLR7/8 agonist adjuvant formulation.

Vaccination can help prevent infection and can also be used to treat cancer, allergy, and potentially even drug overdose. Adjuvants enhance vaccine responses, but currently, the path to their advancement and development is incremental. We used a phenotypic small-molecule screen using THP-1 cells to identify nuclear factor-κB (NF-κB)-activating molecules followed by counterscreening lead target libraries with a quantitative tumor necrosis factor immunoassay using primary human peripheral blood mononuclear cells. Screening on primary cells identified an imidazopyrimidine, dubbed PVP-037. Moreover, while PVP-037 did not overtly activate THP-1 cells, it demonstrated broad innate immune activation, including NF-κB and cytokine induction from primary human leukocytes in vitro as well as enhancement of influenza and SARS-CoV-2 antigen-specific humoral responses in mice. Several de novo synthesis structural enhancements iteratively improved PVP-037's in vitro efficacy, potency, species-specific activity, and in vivo adjuvanticity. Overall, we identified imidazopyrimidine Toll-like receptor-7/8 adjuvants that act in synergy with oil-in-water emulsion to enhance immune responses.

Biological Age, Chronological Age, and Survival in Pulmonary Fibrosis: A Causal Mediation Analysis.

Rationale: Accelerated biological aging has been implicated in the development of interstitial lung disease (ILD) and other diseases of aging but remains poorly understood. Objectives: To identify plasma proteins that mediate the relationship between chronological age and survival association in patients with ILD. Methods: Causal mediation analysis was performed to identify plasma proteins that mediated the chronological age-survival relationship in an idiopathic pulmonary fibrosis discovery cohort. Proteins mediating this relationship after adjustment for false discovery were advanced for testing in an independent ILD validation cohort and explored in a chronic obstructive pulmonary disease cohort. A proteomic-based measure of biological age was constructed and survival analysis performed, assessing the impact of biological age and peripheral blood telomere length on the chronological age-survival relationship. Measurements and Main Results: Twenty-two proteins mediated the chronological age-survival relationship after adjustment for false discovery in the idiopathic pulmonary fibrosis discovery cohort (n = 874), with 19 remaining significant mediators of this relationship in the ILD validation cohort (n = 983) and one mediating this relationship in the chronic obstructive pulmonary disease cohort. Latent transforming growth factor-ß binding protein 2 and ectodysplasin A2 receptor showed the strongest mediation across cohorts. A proteomic measure of biological age completely attenuated the chronological age-survival association and better discriminated survival than chronological age. Results were robust to adjustment for peripheral blood telomere length, which did not mediate the chronological age-survival relationship. Conclusions: Molecular measures of aging completely mediate the relationship between chronological age and survival, suggesting that chronological age has no direct effect on ILD survival.

Effect of apolipoprotein E ε4 and its modification by sociodemographic characteristics on cognitive measures in South Asians from LASI-DAD.

BACKGROUND: We investigated the effects of apolipoprotein E (APOE) ε4 and its interactions with sociodemographic characteristics on cognitive measures in South Asians from the Diagnostic Assessment of Dementia for the Longitudinal Aging Study of India (LASI-DAD). METHODS: Linear regression was used to assess the association between APOE ε4 and global- and domain-specific cognitive function in 2563 participants (mean age 69.6 ± 7.3 years; 53% female). Effect modification by age, sex, and education were explored using interaction terms and subgroup analyses. RESULTS: APOE ε4 was inversely associated with most cognitive measures (p < 0.05). This association was stronger with advancing age for the Hindi Mental State Examination (HMSE) score (ßε4×age = -0.44, p = 0.03), orientation (ßε4×age = -0.07, p = 0.01), and language/fluency (ßε4×age = -0.07, p = 0.01), as well as in females for memory (ßε4×male = 0.17, p = 0.02) and language/fluency (ßε4×male = 0.12, p = 0.03). DISCUSSION: APOE Îµ4 is associated with lower cognitive function in South Asians from India, with a more pronounced impact observed in females and older individuals. HIGHLIGHTS: APOE Îµ4 carriers had lower global and domain-specific cognitive performance. Females and older individuals may be more susceptible to ε4 effects. For most cognitive measures, there was no interaction between ε4 and education.

DNA Methylation at C-Reactive Protein-Associated CpG Sites May Mediate the Pathway Between Educational Attainment and Cognition.

Growing evidence has linked inflammatory processes to cognitive decline and dementia. This work examines whether an epigenetic marker of C-reactive protein (CRP), a common clinical inflammatory biomarker, may mediate the relationship between educational attainment and cognition. We first evaluated whether 53 previously reported CRP-associated DNA methylation sites (CpGs) are associated with CRP, both individually and aggregated into a methylation risk score (MRSCRP), in 3 298 participants from the Health and Retirement Study (HRS, mean age = 69.7 years). Forty-nine CpGs (92%) were associated with the natural logarithm of CRP in HRS after adjusting for age, sex, smoking, BMI, genetic ancestry, and white blood cell counts (p < .05), and each standard deviation increase in MRSCRP was associated with a 0.38 unit increase in lnCRP (p = 4.02E-99). In cross-sectional analysis, for each standard deviation increase in MRSCRP, total memory score and total cognitive score decreased, on average, by 0.28 words and 0.43 items, respectively (p < .001). Further, MRSCRP mediated 6.9% of the relationship between high school education and total memory score in a model adjusting for age, sex, and genetic ancestry (p < .05); this was attenuated to 2.4% with additional adjustment for marital status, APOE ε4 status, health behaviors, and comorbidities (p < .05). Thus, CRP-associated methylation may partially mediate the relationship between education and cognition at older ages. Further research is warranted to determine whether DNA methylation at these sites may improve current prediction models for cognitive impairment in older adults.

Aire mediates tolerance to insulin through thymic trimming of high-affinity T cell clones.

Insulin is a central autoantigen in the pathogenesis of T1D, and thymic epithelial cell expression of insulin under the control of the Autoimmune Regulator (Aire) is thought to be a key component of maintaining tolerance to insulin. In spite of this general working model, direct detection of this thymic selection on insulin-specific T cells has been somewhat elusive. Here, we used a combination of highly sensitive T cell receptor transgenic models for detecting thymic selection and sorting and sequencing of Insulin-specific CD4+ T cells from Aire-deficient mice as a strategy to further define their selection. This analysis revealed a number of unique t cell receptor (TCR) clones in Aire-deficient hosts with high affinity for insulin/major histocompatibility complex (MHC) ligands. We then modeled the thymic selection of one of these clones in Aire-deficient versus wild-type hosts and found that this model clone could escape thymic negative selection in the absence of thymic Aire. Together, these results suggest that thymic expression of insulin plays a key role in trimming and removing high-affinity insulin-specific T cells from the repertoire to help promote tolerance.

Association analysis between an epigenetic risk score and blood pressure.

Background: Epigenome-wide association studies have revealed multiple DNA methylation sites (CpGs) associated with alcohol consumption, an important lifestyle risk factor for cardiovascular diseases. Results: We generated an alcohol consumption epigenetic risk score (ERS) based on previously reported 144 alcohol-associated CpGs and examined the association of the ERS with systolic blood pressure (SBP), diastolic blood pressure (DBP), and hypertension (HTN) in 3,898 Framingham Heart Study (FHS) participants. We found an association of alcohol intake with the ERS in the meta-analysis with 0.09 units higher ERS per drink consumed per day (p < 0.0001). Cross-sectional analyses in FHS revealed that a one-unit increment of the ERS was associated with 1.93 mm Hg higher SBP (p = 4.64E-07), 0.68 mm Hg higher DBP (p = 0.006), and an odds ratio of 1.78 for HTN (p < 2E-16). Meta-analysis of the cross-sectional association of the ERS with BP traits in eight independent external cohorts (n = 11,544) showed similar relationships with blood pressure levels, i.e., a one-unit increase in ERS was associated with 0.74 (p = 0.002) and 0.50 (p = 0.0006) mm Hg higher SBP and DBP, but could not confirm the association with hypertension. Longitudinal analyses in FHS (n = 3,260) and five independent external cohorts (n = 4,021) showed that the baseline ERS was not associated with a change in blood pressure over time or with incident HTN. Conclusions: Our findings provide proof-of-concept that utilizing an ERS is a useful approach to capture the recent health consequences of lifestyle behaviors such as alcohol consumption.

Supplementary feeding of birds during the winter influences measures of avian community structure in yards in a subtropical city.

Supplementary feeding, the intentional provision of food to wild birds is a common activity in developed nations during the winter. The energy inputs represented by supplementary feeding are vast, and thus it is likely an important mechanism shaping bird communities in urban areas. However, research in this regard has mainly occurred in temperate and non-urban settings. Moreover, few studies have been informed by supplementary feeding habits of local community members limiting their inference. We evaluated the effects of two commonly provided wild bird foods on the abundance and species diversity of birds in yards over two winters in San Antonio, Texas, United States, a city located in a subtropical region. We used a reversed Before-After-Control-Impact experimental design in which yards were randomly allocated either mixed seed, Nyjer, or no food (control) between November 2019 and March 2020 (Year One). Between November 2020 and March 2021 (Year Two) supplementary food was not provided in any yards. Point counts conducted during both years of the study revealed that overall bird abundance was consistent between years in control yards and yards provided with Nyjer. In contrast, overall bird abundance was statistically significantly higher when supplementary food was present in mixed seed yards, driven by an increase in granivorous and omnivorous species. In contrast, supplementary feeding had no statistically significant effect on the abundance of insectivorous species or on species diversity, although species diversity tended to be higher in the presence of mixed seed. Our study demonstrates that wild bird food commonly provided by community members influences measures of avian community structure during the winter in urban yards in a subtropical city. However, these results depend on the type of bird food provided. Our results provide insight into the processes underlying the effects of urbanization on bird communities, and thus have implications for the management of urban birds more broadly.

Racial and Socioeconomic Disparities Exist in Patients Pursuing Peanut Oral Immunotherapy.

BACKGROUND: Oral immunotherapy (OIT) is a promising treatment for food allergy. Prior studies demonstrate significant differences among food-allergic individuals across race, ethnicity, and socioeconomic groups. Disparities in OIT have not been evaluated. OBJECTIVE: We assessed disparities in the use of OIT in patients with peanut allergy based on race, ethnicity, and socioeconomic status at a single academic medical center. METHODS: We identified 1028 peanut-allergic patients younger than 18 years receiving care in the University of Michigan food allergy clinics. Of these, 148 patients who underwent peanut OIT (treatment group) were compared with the 880 patients who avoided peanut (control group). Pertinent demographic and socioeconomic characteristics were compared. RESULTS: There were no differences in gender or ethnicity between the OIT and control groups. However, Black patients comprised 18% of the control group but only 4.1% of the OIT treatment group (P < .0001). The proportion of patients with private insurance was significantly higher in the treatment group compared with the control group (93.2% vs 82.2%, P = .0004). Finally, the neighborhood affluence index, a census-based measure of the relative socioeconomic prosperity of a neighborhood, was significantly higher in the OIT group than the control group (0.51 ± 0.18 vs 0.47 ± 0.19, P = .015), whereas the neighborhood disadvantage index, a census-based measure of the relative socioeconomic disadvantage of a neighborhood, was significantly lower (0.082 ± 0.062 vs 0.10 ± 0.093, P = .020). CONCLUSIONS: Significant racial and economic disparities exist at our institution between peanut-allergic individuals who receive OIT and those who do not. Efforts to understand the basis for these disparities are important to ensure that patients have equitable access to OIT.

Epigenetic age acceleration is associated with blood lipid levels in a multi-ancestry sample of older U.S. adults.

BACKGROUND: Dyslipidemia, which is characterized by an unfavorable lipid profile, is a key risk factor for cardiovascular disease (CVD). Understanding the relationships between epigenetic aging and lipid levels may help guide early prevention and treatment efforts for dyslipidemia. METHODS: We used weighted linear regression to cross-sectionally investigate the associations between five measures of epigenetic age acceleration estimated from whole blood DNA methylation (HorvathAge Acceleration, HannumAge Acceleration, PhenoAge Acceleration, GrimAge Acceleration, and DunedinPACE) and four blood lipid measures (total cholesterol (TC), LDL-C, HDL-C, and triglycerides (TG)) in 3,813 participants (mean age = 70 years) from the Health and Retirement Study (HRS). As a sensitivity analysis, we examined the same associations in participants who fasted prior to the blood draw (n = 2,531) and in participants who did not take lipid-lowering medication (n = 1,869). Using interaction models, we also examined whether demographic factors including age, sex, and educational attainment modified the relationships between epigenetic age acceleration and blood lipids. RESULTS: After adjusting for age, race/ethnicity, sex, fasting status, and lipid-lowering medication use, greater epigenetic age acceleration was associated with lower TC, HDL-C, and LDL-C, and higher TG (p < 0.05), although the effect sizes were relatively small (e.g., < 7 mg/dL of TC per standard deviation in epigenetic age acceleration). GrimAge acceleration and DunedinPACE associations with all lipids remained significant after further adjustment for body mass index, smoking status, and educational attainment. These associations were stronger in participants who fasted and who did not use lipid-lowering medication, particularly for LDL-C. We observed the largest number of interactions between DunedinPACE and demographic factors, where the associations with lipids were stronger in younger participants, females, and those with higher educational attainment. CONCLUSION: Multiple measures of epigenetic age acceleration are associated with blood lipid levels in older adults. A greater understanding of how these associations differ across demographic groups can help shed light on the relationships between aging and downstream cardiovascular diseases. The inverse associations between epigenetic age and TC and LDL-C could be due to sample limitations or non-linear relationships between age and these lipids, as both TC and LDL-C decrease faster at older ages.

Machine learning models for predicting blood pressure phenotypes by combining multiple polygenic risk scores.

We construct non-linear machine learning (ML) prediction models for systolic and diastolic blood pressure (SBP, DBP) using demographic and clinical variables and polygenic risk scores (PRSs). We developed a two-model ensemble, consisting of a baseline model, where prediction is based on demographic and clinical variables only, and a genetic model, where we also include PRSs. We evaluate the use of a linear versus a non-linear model at both the baseline and the genetic model levels and assess the improvement in performance when incorporating multiple PRSs. We report the ensemble model's performance as percentage variance explained (PVE) on a held-out test dataset. A non-linear baseline model improved the PVEs from 28.1 to 30.1% (SBP) and 14.3% to 17.4% (DBP) compared with a linear baseline model. Including seven PRSs in the genetic model computed based on the largest available GWAS of SBP/DBP improved the genetic model PVE from 4.8 to 5.1% (SBP) and 4.7 to 5% (DBP) compared to using a single PRS. Adding additional 14 PRSs computed based on two independent GWASs further increased the genetic model PVE to 6.3% (SBP) and 5.7% (DBP). PVE differed across self-reported race/ethnicity groups, with primarily all non-White groups benefitting from the inclusion of additional PRSs. In summary, non-linear ML models improves BP prediction in models incorporating diverse populations.

Methylation patterns associated with C-reactive protein in racially and ethnically diverse populations.

Systemic low-grade inflammation is a feature of chronic disease. C-reactive protein (CRP) is a common biomarker of inflammation and used as an indicator of disease risk; however, the role of inflammation in disease is not completely understood. Methylation is an epigenetic modification in the DNA which plays a pivotal role in gene expression. In this study we evaluated differential DNA methylation patterns associated with blood CRP level to elucidate biological pathways and genetic regulatory mechanisms to improve the understanding of chronic inflammation. The racially and ethnically diverse participants in this study were included as 50% White, 41% Black or African American, 7% Hispanic or Latino/a, and 2% Native Hawaiian, Asian American, American Indian, or Alaska Native (total n = 13,433) individuals. We replicated 113 CpG sites from 87 unique loci, of which five were novel (CADM3, NALCN, NLRC5, ZNF792, and cg03282312), across a discovery set of 1,150 CpG sites associated with CRP level (p < 1.2E-7). The downstream pathways affected by DNA methylation included the identification of IFI16 and IRF7 CpG-gene transcript pairs which contributed to the innate immune response gene enrichment pathway along with NLRC5, NOD2, and AIM2. Gene enrichment analysis also identified the nuclear factor-kappaB transcription pathway. Using two-sample Mendelian randomization (MR) we inferred methylation at three CpG sites as causal for CRP levels using both White and Black or African American MR instrument variables. Overall, we identified novel CpG sites and gene transcripts that could be valuable in understanding the specific cellular processes and pathogenic mechanisms involved in inflammation.

Epigenetic age acceleration is associated with blood lipid levels in a multi-ancestry sample of older U.S. adults.

Background: Dyslipidemia, which is characterized by an unfavorable lipid profile, is a key risk factor for cardiovascular disease (CVD). Understanding the relationships between epigenetic aging and lipid levels may help guide early prevention and treatment efforts for dyslipidemia. Methods: We used weighted linear regression to cross-sectionally investigate the associations between five measures of epigenetic age acceleration estimated from whole blood DNA methylation (HorvathAge Acceleration, HannumAge Acceleration, PhenoAge Acceleration, GrimAge Acceleration, and DunedinPACE) and four blood lipid measures (total cholesterol (TC), LDL-C, HDL-C, and triglycerides (TG)) in 3,813 participants (mean age = 70 years) from the Health and Retirement Study (HRS). As a sensitivity analysis, we examined the same associations in participants who fasted prior to the blood draw (n = and f) and in participants who did not take lipid-lowering medication (n = 1,869). Using interaction models, we also examined whether the relationships between epigenetic age acceleration and blood lipids differ by demographic factors including age, sex, and educational attainment. Results: After adjusting for age, race/ethnicity, sex, fasting status, and lipid-lowering medication use, greater epigenetic age acceleration was associated with lower TC, HDL-C, and LDL-C, and higher TG (p < 0.05). GrimAge acceleration and DunedinPACE associations with all lipids remained significant after further adjusting for body mass index, smoking status, and educational attainment. These associations were stronger in participants who fasted and who did not use lipid-lowering medication, particularly for LDL-C. We observed the largest number of interactions between DunedinPACE and demographic factors, where the associations with lipids were stronger in younger participants, females, and those with higher educational attainment. Conclusion: Epigenetic age acceleration, a powerful biomarker of cellular aging, is highly associated with blood lipid levels in older adults. A greater understanding of how these associations differ across demographic groups can help shed light on the relationships between aging and downstream cardiovascular diseases. The inverse associations between epigenetic age and TC and LDL-C could be due to sample limitations or the non-linear relationship between age and these lipids, as both TC and LDL-C decrease faster at older ages. More studies are needed to further understand the temporal relationships between epigenetic age acceleration on blood lipids and other health outcomes.

Epigenome-wide association study of long-term psychosocial stress in older adults.

Long-term psychosocial stress is strongly associated with negative physical and mental health outcomes, as well as adverse health behaviours; however, little is known about the role that stress plays on the epigenome. One proposed mechanism by which stress affects DNA methylation is through health behaviours. We conducted an epigenome-wide association study (EWAS) of cumulative psychosocial stress (n = 2,689) from the Health and Retirement Study (mean age = 70.4 years), assessing DNA methylation (Illumina Infinium HumanMethylationEPIC Beadchip) at 789,656 CpG sites. For identified CpG sites, we conducted a formal mediation analysis to examine whether smoking, alcohol use, physical activity, and body mass index (BMI) mediate the relationship between stress and DNA methylation. Nine CpG sites were associated with psychosocial stress (all p < 9E-07; FDR q < 0.10). Additionally, health behaviours and/or BMI mediated 9.4% to 21.8% of the relationship between stress and methylation at eight of the nine CpGs. Several of the identified CpGs were in or near genes associated with cardiometabolic traits, psychosocial disorders, inflammation, and smoking. These findings support our hypothesis that psychosocial stress is associated with DNA methylation across the epigenome. Furthermore, specific health behaviours mediate only a modest percentage of this relationship, providing evidence that other mechanisms may link stress and DNA methylation.

Association analysis between an epigenetic alcohol risk score and blood pressure.

Background: Epigenome-wide association studies have revealed multiple DNA methylation sites (CpGs) associated with alcohol consumption, an important lifestyle risk factor for cardiovascular diseases. Results: We generated an alcohol consumption epigenetic risk score (ERS) based on previously reported 144 alcohol-associated CpGs and examined the association of the ERS with systolic blood pressure (SBP), diastolic blood pressure (DBP), and hypertension (HTN) in 3,898 Framingham Heart Study (FHS) participants. We found an association of alcohol intake with the ERS in the meta-analysis with 0.09 units higher ERS per drink consumed per day (p < 0.0001). Cross-sectional analyses in FHS revealed that a one-unit increment of the ERS was associated with 1.93 mm Hg higher SBP (p = 4.64E-07), 0.68 mm Hg higher DBP (p = 0.006), and an odds ratio of 1.78 for HTN (p < 2E-16). Meta-analysis of the cross-sectional association of the ERS with BP traits in eight independent external cohorts (n = 11,544) showed similar relationships with blood pressure levels, i.e., a one-unit increase in ERS was associated with 0.74 (p = 0.002) and 0.50 (p = 0.0006) mm Hg higher SBP and DBP, but could not confirm the association with hypertension. Longitudinal analyses in FHS (n = 3,260) and five independent external cohorts (n = 4,021) showed that the baseline ERS was not associated with a change in blood pressure over time or with incident HTN. Conclusions: Our findings provide proof-of-concept that utilizing an ERS is a useful approach to capture the recent health consequences of lifestyle behaviors such as alcohol consumption.

50,000 years of Evolutionary History of India: Insights from ~2,700 Whole Genome Sequences.

India has been underrepresented in whole genome sequencing studies. We generated 2,762 high coverage genomes from India-including individuals from most geographic regions, speakers of all major languages, and tribal and caste groups-providing a comprehensive survey of genetic variation in India. With these data, we reconstruct the evolutionary history of India through space and time at fine scales. We show that most Indians derive ancestry from three ancestral groups related to ancient Iranian farmers, Eurasian Steppe pastoralists and South Asian hunter-gatherers. We uncover a common source of Iranian-related ancestry from early Neolithic cultures of Central Asia into the ancestors of Ancestral South Indians (ASI), Ancestral North Indians (ANI), Austro-asiatic-related and East Asian-related groups in India. Following these admixtures, India experienced a major demographic shift towards endogamy, resulting in extensive homozygosity and identity-by-descent sharing among individuals. At deep time scales, Indians derive around 1-2% of their ancestry from gene flow from archaic hominins, Neanderthals and Denisovans. By assembling the surviving fragments of archaic ancestry in modern Indians, we recover ~1.5 Gb (or 50%) of the introgressing Neanderthal and ~0.6 Gb (or 20%) of the introgressing Denisovan genomes, more than any other previous archaic ancestry study. Moreover, Indians have the largest variation in Neanderthal ancestry, as well as the highest amount of population-specific Neanderthal segments among worldwide groups. Finally, we demonstrate that most of the genetic variation in Indians stems from a single major migration out of Africa that occurred around 50,000 years ago, with minimal contribution from earlier migration waves. Together, these analyses provide a detailed view of the population history of India and underscore the value of expanding genomic surveys to diverse groups outside Europe.