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 under 18 years of age receiving care in the University of Michigan food allergy clinics. 148 patients undergoing peanut OIT (treatment group) were compared to the 880 patients avoiding 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<0.0001). The proportion of patients with private insurance was significantly higher in the treatment group compared to control group, 93.2% vs 82.2% (p=0.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 vs. the control group (0.51±0.18 vs 0.47±0.19) (p=0.015), while 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=0.020). CONCLUSION: 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 patients have equitable access to OIT.

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.

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.

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.

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.

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.

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.

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.

Tracking the role of Aire in immune tolerance to the eye with a TCR transgenic mouse model.

Roughly one-half of mice with partial defects in two immune tolerance pathways (AireGW/+Lyn-/- mice) spontaneously develop severe damage to their retinas due to T cell reactivity to Aire-regulated interphotoreceptor retinoid-binding protein (IRBP). Single-cell T cell receptor (TCR) sequencing of CD4+ T cells specific for a predominate epitope of IRBP showed a remarkable diversity of autoantigen-specific TCRs with greater clonal expansions in mice with disease. TCR transgenic mice made with an expanded IRBP-specific TCR (P2.U2) of intermediate affinity exhibited strong but incomplete negative selection of thymocytes. This negative selection was absent in IRBP-/- mice and greatly defective in AireGW/+ mice. Most P2.U2+/- mice and all P2.U.2+/-AireGW/+ mice rapidly developed inflammation of the retina and adjacent uvea (uveitis). Aire-dependent IRBP expression in the thymus also promoted Treg differentiation, but the niche for this fate determination was small, suggesting differences in antigen presentation leading to negative selection vs. thymic Treg differentiation and a stronger role for negative selection in preventing autoimmune disease in the retina.

Multi-omics and pathway analyses of genome-wide associations implicate regulation and immunity in verbal declarative memory performance.

BACKGROUND: Uncovering the functional relevance underlying verbal declarative memory (VDM) genome-wide association study (GWAS) results may facilitate the development of interventions to reduce age-related memory decline and dementia. METHODS: We performed multi-omics and pathway enrichment analyses of paragraph (PAR-dr) and word list (WL-dr) delayed recall GWAS from 29,076 older non-demented individuals of European descent. We assessed the relationship between single-variant associations and expression quantitative trait loci (eQTLs) in 44 tissues and methylation quantitative trait loci (meQTLs) in the hippocampus. We determined the relationship between gene associations and transcript levels in 53 tissues, annotation as immune genes, and regulation by transcription factors (TFs) and microRNAs. To identify significant pathways, gene set enrichment was tested in each cohort and meta-analyzed across cohorts. Analyses of differential expression in brain tissues were conducted for pathway component genes. RESULTS: The single-variant associations of VDM showed significant linkage disequilibrium (LD) with eQTLs across all tissues and meQTLs within the hippocampus. Stronger WL-dr gene associations correlated with reduced expression in four brain tissues, including the hippocampus. More robust PAR-dr and/or WL-dr gene associations were intricately linked with immunity and were influenced by 31 TFs and 2 microRNAs. Six pathways, including type I diabetes, exhibited significant associations with both PAR-dr and WL-dr. These pathways included fifteen MHC genes intricately linked to VDM performance, showing diverse expression patterns based on cognitive status in brain tissues. CONCLUSIONS: VDM genetic associations influence expression regulation via eQTLs and meQTLs. The involvement of TFs, microRNAs, MHC genes, and immune-related pathways contributes to VDM performance in older individuals.

Epigenetic Aging Is Associated With Measures of Midlife Muscle Volume and Attenuation in CARDIA Study.

BACKGROUND: GrimAge acceleration (GAA), an epigenetic marker that represents physiologic aging, is associated with age-related diseases including cancer and cardiovascular diseases. However, the associations between GAA and muscle mass and function are unknown. METHODS: We estimated measures of GAA in 1 118 Black and White participants from the Coronary Artery Risk Development in Young Adults (CARDIA) Study at exam years (Y) 15 (2000-2001) and 20 (2005-2006). Abdominal muscle composition was measured using CT scans at the Y25 (2010-2011) visit. We used multivariate regression models to examine associations of GAA estimates with muscle imaging measurements. RESULTS: In the CARDIA study, each 1-year higher GAA was associated with an average 1.1% (95% confidence interval [CI]: 0.6%, 1.5%) higher intermuscular adipose tissue (IMAT) volume for abdominal muscles. Each 1-year higher GAA was associated with an average -0.089 Hounsfield unit (HU; 95% CI: -0.146, -0.032) lower lean muscle attenuation and an average -0.049 HU (95% CI: -0.092, -0.007) lower IMAT attenuation for abdominal muscles. Stratified analyses showed that GAA was more strongly associated with higher abdominal muscle IMAT volume in females and significantly associated with lower lean muscle attenuation for White participants only. CONCLUSIONS: Higher GAA is associated with higher abdominal muscle IMAT volume and lower lean muscle attenuation in a midlife population.

Alcohol Use and Mortality Among Older Couples in the United States: Evidence of Individual and Partner Effects.

BACKGROUND AND OBJECTIVES: Spouses with concordant (i.e., similar) drinking behaviors often report better quality marriages and are married longer compared with those who report discordant drinking behaviors. Less is known regarding whether concordant or discordant patterns have implications for health, as couples grow older. The present study examined whether drinking patterns among older couples are associated with mortality over time. RESEARCH DESIGN AND METHODS: The Health and Retirement Study (HRS) is a nationally representative sample of individuals and their partners (married/cohabiting) over age 50 in the United States, in which participants completed surveys every 2 years. Participants included 4,656 married/cohabiting different-sex couples (9,312 individuals) who completed at least 3 waves of the HRS from 1996 to 2016. Participants reported whether they drank alcohol at all in the last 3 months, and if so, the average amount they drank per week. Mortality data were from 2016. RESULTS: Analyses revealed concordant drinking spouses (both indicated they drank in the last 3 months) survived longer than discordant drinking spouses (1 partner drinks and the other does not) and concordant nondrinking spouses. Analysis of average drinks per week showed a quadratic association with mortality such that light drinking predicted better survival rates among individuals and their partners compared with abstaining and heavy drinking. Further, similar levels of drinking in terms of the amount of drinking were associated with greater survival, particularly among wives. DISCUSSION AND IMPLICATIONS: This study moves the field forward by showing that survival varies as a function of one's own and one's partner's drinking.

Evaluation of the utility of rescans in the treatment of prostate and pelvic nodes with pencil beam scanning protons.

To determine the necessity of the first week CT simulation rescan of pencil beam scanning (PBS) prostate patients requiring treatment to the pelvic lymph nodes. Patients were treated on a prospective registry trial sponsored by the Proton Collaborative Group (PCG-NCT01255748). A total of 42 patients with high-risk prostate cancer requiring treatment to the pelvic lymph nodes were evaluated in a single calendar year. The cohort consisted of a mix of intact prostate and postprostatectomy patients. Most of the patients were treated with a simultaneous integrated boost (SIB) approach for the majority of the plan. The radiation prescriptions varied depending on whether the patient had an intact prostate or prostate bed. The plan geometry consisted of two lateral beams and a single field optimization (SFO) dosimetric matching technique using pencil beam scanning proton therapy. An in-house protocol was established wherein all high-risk prostate patients had at least 1 rescan evaluation performed during the first 5 ± 2 fractions, which was used to determine whether the nominal approved plan was robust to daily setup uncertainties and anatomical variations. If the evaluation failed clinical analysis, an adaptive replan was created. If 5% or more of the evaluated rescans resulted in a qualified adaptive plan, the planning technique would be considered insufficient. Of the 42 patients investigated, five (11.9%) required an adaptive plan. As it turned out, all five of these patients would have been rescanned within the first 5 fractions of treatment, independent of the established rescan protocol, due to a physician, dosimetrist, or therapist requesting a rescan to investigate specific areas of concern regarding setup or anatomic changes. Of the 5 adaptive plans, only one (2.4%) meets the criteria of a qualified adaptive plan. Our findings substantiated that this policy of a planned rescan with the 5th fraction was no longer necessary, the dosimetric technique had proven to be robust, and moving forward we will only perform these rescans if there is a significant issue with daily setups or observed changes in anatomy.

Type 2 Diabetes Modifies the Association of CAD Genomic Risk Variants With Subclinical Atherosclerosis.

BACKGROUND: Individuals with type 2 diabetes (T2D) have an increased risk of coronary artery disease (CAD), but questions remain about the underlying pathology. Identifying which CAD loci are modified by T2D in the development of subclinical atherosclerosis (coronary artery calcification [CAC], carotid intima-media thickness, or carotid plaque) may improve our understanding of the mechanisms leading to the increased CAD in T2D. METHODS: We compared the common and rare variant associations of known CAD loci from the literature on CAC, carotid intima-media thickness, and carotid plaque in up to 29 670 participants, including up to 24 157 normoglycemic controls and 5513 T2D cases leveraging whole-genome sequencing data from the Trans-Omics for Precision Medicine program. We included first-order T2D interaction terms in each model to determine whether CAD loci were modified by T2D. The genetic main and interaction effects were assessed using a joint test to determine whether a CAD variant, or gene-based rare variant set, was associated with the respective subclinical atherosclerosis measures and then further determined whether these loci had a significant interaction test. RESULTS: Using a Bonferroni-corrected significance threshold of P<1.6×10-4, we identified 3 genes (ATP1B1, ARVCF, and LIPG) associated with CAC and 2 genes (ABCG8 and EIF2B2) associated with carotid intima-media thickness and carotid plaque, respectively, through gene-based rare variant set analysis. Both ATP1B1 and ARVCF also had significantly different associations for CAC in T2D cases versus controls. No significant interaction tests were identified through the candidate single-variant analysis. CONCLUSIONS: These results highlight T2D as an important modifier of rare variant associations in CAD loci with CAC.

Circulating metabolites may illustrate relationship of alcohol consumption with cardiovascular disease.

BACKGROUND: Metabolite signatures of long-term alcohol consumption are lacking. To better understand the molecular basis linking alcohol drinking and cardiovascular disease (CVD), we investigated circulating metabolites associated with long-term alcohol consumption and examined whether these metabolites were associated with incident CVD. METHODS: Cumulative average alcohol consumption (g/day) was derived from the total consumption of beer, wine, and liquor on average of 19 years in 2428 Framingham Heart Study Offspring participants (mean age 56 years, 52% women). We used linear mixed models to investigate the associations of alcohol consumption with 211 log-transformed plasma metabolites, adjusting for age, sex, batch, smoking, diet, physical activity, BMI, and familial relationship. Cox models were used to test the association of alcohol-related metabolite scores with fatal and nonfatal incident CVD (myocardial infarction, coronary heart disease, stroke, and heart failure). RESULTS: We identified 60 metabolites associated with cumulative average alcohol consumption (p < 0.05/211 ≈ 0.00024). For example, 1 g/day increase of alcohol consumption was associated with higher levels of cholesteryl esters (e.g., CE 16:1, beta = 0.023 ± 0.002, p = 6.3e - 45) and phosphatidylcholine (e.g., PC 32:1, beta = 0.021 ± 0.002, p = 3.1e - 38). Survival analysis identified that 10 alcohol-associated metabolites were also associated with a differential CVD risk after adjusting for age, sex, and batch. Further, we built two alcohol consumption weighted metabolite scores using these 10 metabolites and showed that, with adjustment age, sex, batch, and common CVD risk factors, the two scores had comparable but opposite associations with incident CVD, hazard ratio 1.11 (95% CI = [1.02, 1.21], p = 0.02) vs 0.88 (95% CI = [0.78, 0.98], p = 0.02). CONCLUSIONS: We identified 60 long-term alcohol consumption-associated metabolites. The association analysis with incident CVD suggests a complex metabolic basis between alcohol consumption and CVD.

Methods for mediation analysis with high-dimensional DNA methylation data: Possible choices and comparisons.

Epigenetic researchers often evaluate DNA methylation as a potential mediator of the effect of social/environmental exposures on a health outcome. Modern statistical methods for jointly evaluating many mediators have not been widely adopted. We compare seven methods for high-dimensional mediation analysis with continuous outcomes through both diverse simulations and analysis of DNAm data from a large multi-ethnic cohort in the United States, while providing an R package for their seamless implementation and adoption. Among the considered choices, the best-performing methods for detecting active mediators in simulations are the Bayesian sparse linear mixed model (BSLMM) and high-dimensional mediation analysis (HDMA); while the preferred methods for estimating the global mediation effect are high-dimensional linear mediation analysis (HILMA) and principal component mediation analysis (PCMA). We provide guidelines for epigenetic researchers on choosing the best method in practice and offer suggestions for future methodological development.

Association Between Whole Blood-Derived Mitochondrial DNA Copy Number, Low-Density Lipoprotein Cholesterol, and Cardiovascular Disease Risk.

Background The relationship between mitochondrial DNA copy number (mtDNA CN) and cardiovascular disease remains elusive. Methods and Results We performed cross-sectional and prospective association analyses of blood-derived mtDNA CN and cardiovascular disease outcomes in 27 316 participants in 8 cohorts of multiple racial and ethnic groups with whole-genome sequencing. We also performed Mendelian randomization to explore causal relationships of mtDNA CN with coronary heart disease (CHD) and cardiometabolic risk factors (obesity, diabetes, hypertension, and hyperlipidemia). P<0.01 was used for significance. We validated most of the previously reported associations between mtDNA CN and cardiovascular disease outcomes. For example, 1-SD unit lower level of mtDNA CN was associated with 1.08 (95% CI, 1.04-1.12; P<0.001) times the hazard for developing incident CHD, adjusting for covariates. Mendelian randomization analyses showed no causal effect from a lower level of mtDNA CN to a higher CHD risk (ß=0.091; P=0.11) or in the reverse direction (ß=-0.012; P=0.076). Additional bidirectional Mendelian randomization analyses revealed that low-density lipoprotein cholesterol had a causal effect on mtDNA CN (ß=-0.084; P<0.001), but the reverse direction was not significant (P=0.059). No causal associations were observed between mtDNA CN and obesity, diabetes, and hypertension, in either direction. Multivariable Mendelian randomization analyses showed no causal effect of CHD on mtDNA CN, controlling for low-density lipoprotein cholesterol level (P=0.52), whereas there was a strong direct causal effect of higher low-density lipoprotein cholesterol on lower mtDNA CN, adjusting for CHD status (ß=-0.092; P<0.001). Conclusions Our findings indicate that high low-density lipoprotein cholesterol may underlie the complex relationships between mtDNA CN and vascular atherosclerosis.

Mosaic chromosomal alterations in blood across ancestries using whole-genome sequencing.

Megabase-scale mosaic chromosomal alterations (mCAs) in blood are prognostic markers for a host of human diseases. Here, to gain a better understanding of mCA rates in genetically diverse populations, we analyzed whole-genome sequencing data from 67,390 individuals from the National Heart, Lung, and Blood Institute Trans-Omics for Precision Medicine program. We observed higher sensitivity with whole-genome sequencing data, compared with array-based data, in uncovering mCAs at low mutant cell fractions and found that individuals of European ancestry have the highest rates of autosomal mCAs and the lowest rates of chromosome X mCAs, compared with individuals of African or Hispanic ancestry. Although further studies in diverse populations will be needed to replicate our findings, we report three loci associated with loss of chromosome X, associations between autosomal mCAs and rare variants in DCPS, ADM17, PPP1R16B and TET2 and ancestry-specific variants in ATM and MPL with mCAs in cis.

Multi-ancestry genome-wide study identifies effector genes and druggable pathways for coronary artery calcification.

Coronary artery calcification (CAC), a measure of subclinical atherosclerosis, predicts future symptomatic coronary artery disease (CAD). Identifying genetic risk factors for CAC may point to new therapeutic avenues for prevention. Currently, there are only four known risk loci for CAC identified from genome-wide association studies (GWAS) in the general population. Here we conducted the largest multi-ancestry GWAS meta-analysis of CAC to date, which comprised 26,909 individuals of European ancestry and 8,867 individuals of African ancestry. We identified 11 independent risk loci, of which eight were new for CAC and five had not been reported for CAD. These new CAC loci are related to bone mineralization, phosphate catabolism and hormone metabolic pathways. Several new loci harbor candidate causal genes supported by multiple lines of functional evidence and are regulators of smooth muscle cell-mediated calcification ex vivo and in vitro. Together, these findings help refine the genetic architecture of CAC and extend our understanding of the biological and potential druggable pathways underlying CAC.