Selection, optimization and validation of ten chronic disease polygenic risk scores for clinical implementation in diverse US populations.

Polygenic risk scores (PRSs) have improved in predictive performance, but several challenges remain to be addressed before PRSs can be implemented in the clinic, including reduced predictive performance of PRSs in diverse populations, and the interpretation and communication of genetic results to both providers and patients. To address these challenges, the National Human Genome Research Institute-funded Electronic Medical Records and Genomics (eMERGE) Network has developed a framework and pipeline for return of a PRS-based genome-informed risk assessment to 25,000 diverse adults and children as part of a clinical study. From an initial list of 23 conditions, ten were selected for implementation based on PRS performance, medical actionability and potential clinical utility, including cardiometabolic diseases and cancer. Standardized metrics were considered in the selection process, with additional consideration given to strength of evidence in African and Hispanic populations. We then developed a pipeline for clinical PRS implementation (score transfer to a clinical laboratory, validation and verification of score performance), and used genetic ancestry to calibrate PRS mean and variance, utilizing genetically diverse data from 13,475 participants of the All of Us Research Program cohort to train and test model parameters. Finally, we created a framework for regulatory compliance and developed a PRS clinical report for return to providers and for inclusion in an additional genome-informed risk assessment. The initial experience from eMERGE can inform the approach needed to implement PRS-based testing in diverse clinical settings.

Selection, optimization, and validation of ten chronic disease polygenic risk scores for clinical implementation in diverse populations.

Polygenic risk scores (PRS) have improved in predictive performance supporting their use in clinical practice. Reduced predictive performance of PRS in diverse populations can exacerbate existing health disparities. The NHGRI-funded eMERGE Network is returning a PRS-based genome-informed risk assessment to 25,000 diverse adults and children. We assessed PRS performance, medical actionability, and potential clinical utility for 23 conditions. Standardized metrics were considered in the selection process with additional consideration given to strength of evidence in African and Hispanic populations. Ten conditions were selected with a range of high-risk thresholds: atrial fibrillation, breast cancer, chronic kidney disease, coronary heart disease, hypercholesterolemia, prostate cancer, asthma, type 1 diabetes, obesity, and type 2 diabetes. We developed a pipeline for clinical PRS implementation, used genetic ancestry to calibrate PRS mean and variance, created a framework for regulatory compliance, and developed a PRS clinical report. eMERGE's experience informs the infrastructure needed to implement PRS-based implementation in diverse clinical settings.

More appendicular lean mass relative to body mass index is associated with lower incident diabetes in middle-aged adults in the CARDIA study.

BACKGROUND AND AIMS: Although lower lean mass is associated with greater diabetes prevalence in cross-sectional studies, prospective data specifically in middle-aged Black and White adults are lacking. Relative appendicular lean mass (ALM), such as ALM adjusted for body mass index (BMI), is important to consider since muscle mass is associated with overall body size. We investigated whether ALM/BMI is associated with incident type 2 diabetes in the Coronary Artery Risk Development in Young Adults study. METHODS AND RESULTS: 1893 middle-aged adults (55% women) were included. ALM was measured by DXA in 2005-06. Incident type 2 diabetes was defined in 2010-11 or 2015-16 as fasting glucose ≥7 mmol/L (126 mg/dL), 2-h glucose on OGTT ≥11.1 mmol/L (200 mg/dL) (2010-11 only), HbA1C ≥48 mmol/mol (6.5%) (2010-11 only), or glucose-lowering medications. Cox regression models with sex stratification were performed. In men and women, ALM/BMI was 1.07 ± 0.14 (mean ± SD) and 0.73 ± 0.12, respectively. Seventy men (8.2%) and 71 women (6.8%) developed type 2 diabetes. Per sex-specific SD higher ALM/BMI, unadjusted diabetes risk was lower by 21% in men [HR 0.79 (0.62-0.99), p = 0.04] and 29% in women [HR 0.71 (0.55-0.91), p = 0.008]. After adjusting for age, race, smoking, education, physical activity, and waist circumference, the association was no longer significant. Adjustment for waist circumference accounted for the attenuation in men. CONCLUSION: Although more appendicular lean mass relative to BMI is associated with lower incident type 2 diabetes in middle-aged men and women over 10 years, its effect may be through other metabolic risk factors such as waist circumference, which is a correlate of abdominal fat mass.

Association between muscle mass and diabetes prevalence independent of body fat distribution in adults under 50 years old.

BACKGROUND/OBJECTIVES: Although relatively less muscle mass has been associated with greater diabetes prevalence, whether there is an association between muscle mass and diabetes prevalence independent of body fat distribution is unknown. The objective was to determine whether less skeletal muscle mass is associated with greater diabetes prevalence in young men and women independent of body fat distribution. SUBJECTS/METHODS: One thousand seven hundred and sixty-four adults, aged 20-49 years old, from the United States National Health and Nutrition Examination Survey (2005-2006). Body composition, including appendicular lean mass (ALM), was measured by dual-energy x-ray absorptiometry. Diabetes was defined as fasting blood glucose ≥7 mmol/l, 2-h blood glucose ≥11.1 mmol/l on 75 g OGTT, HbA1c ≥ 48 mmol/mol (6.5%), use of diabetes medications, or self-reported diagnosis of diabetes. RESULTS: The odds of diabetes were 1.31 times higher in men [OR 1.31 (1.18-1.45), p = 0.0001], and 1.24 times higher in women [OR 1.24 (1.05-1.46), p = 0.01], per percent decrease in ALM/weight after controlling for age, race, height, smoking, and education. After additionally controlling for android/gynoid fat, the odds of diabetes were 1.20 times higher per percent decrease in ALM/weight in men [OR 1.20 (1.04-1.37), p = 0.01]; an inverse association was also observed in women, albeit was not statistically significant [OR 1.08 (0.90-1.30), p = 0.42]. CONCLUSIONS: Less muscle mass was associated with greater diabetes prevalence independent of body fat distribution in young men. The association was not statistically significant in women after controlling for android and gynoid adiposity. Low muscle mass could be a causal factor in the development of type 2 diabetes or a correlated marker of higher metabolic risk.

A MUC5B Gene Polymorphism, rs35705950-T, Confers Protective Effects Against COVID-19 Hospitalization but Not Severe Disease or Mortality.

Rationale: A common MUC5B gene polymorphism, rs35705950-T, is associated with idiopathic pulmonary fibrosis (IPF), but its role in severe acute respiratory syndrome coronavirus 2 infection and disease severity is unclear. Objectives: To assess whether rs35705950-T confers differential risk for clinical outcomes associated with coronavirus disease (COVID-19) infection among participants in the Million Veteran Program (MVP). Methods: The MUC5B rs35705950-T allele was directly genotyped among MVP participants; clinical events and comorbidities were extracted from the electronic health records. Associations between the incidence or severity of COVID-19 and rs35705950-T were analyzed within each ancestry group in the MVP followed by transancestry meta-analysis. Replication and joint meta-analysis were conducted using summary statistics from the COVID-19 Host Genetics Initiative (HGI). Sensitivity analyses with adjustment for additional covariates (body mass index, Charlson comorbidity index, smoking, asbestosis, rheumatoid arthritis with interstitial lung disease, and IPF) and associations with post-COVID-19 pneumonia were performed in MVP subjects. Measurements and Main Results: The rs35705950-T allele was associated with fewer COVID-19 hospitalizations in transancestry meta-analyses within the MVP (Ncases = 4,325; Ncontrols = 507,640; OR = 0.89 [0.82-0.97]; P = 6.86 × 10-3) and joint meta-analyses with the HGI (Ncases = 13,320; Ncontrols = 1,508,841; OR, 0.90 [0.86-0.95]; P = 8.99 × 10-5). The rs35705950-T allele was not associated with reduced COVID-19 positivity in transancestry meta-analysis within the MVP (Ncases = 19,168/Ncontrols = 492,854; OR, 0.98 [0.95-1.01]; P = 0.06) but was nominally significant (P < 0.05) in the joint meta-analysis with the HGI (Ncases = 44,820; Ncontrols = 1,775,827; OR, 0.97 [0.95-1.00]; P = 0.03). Associations were not observed with severe outcomes or mortality. Among individuals of European ancestry in the MVP, rs35705950-T was associated with fewer post-COVID-19 pneumonia events (OR, 0.82 [0.72-0.93]; P = 0.001). Conclusions: The MUC5B variant rs35705950-T may confer protection in COVID-19 hospitalizations.

Understanding the Link Between Obesity and Severe COVID-19 Outcomes: Causal Mediation by Systemic Inflammatory Response.

BACKGROUND: Obesity is an established risk factor for severe COVID-19 outcomes. The mechanistic underpinnings of this association are not well-understood. OBJECTIVE: To evaluate the mediating role of systemic inflammation in obesity-associated COVID-19 outcomes. METHODS: This hospital-based, observational study included 3828 SARS-CoV-2-infected patients who were hospitalized February to May 2020 at Massachusetts General Hospital (MGH) or Columbia University Irving Medical Center/New York Presbyterian Hospital (CUIMC/NYP). We use mediation analysis to evaluate whether peak inflammatory biomarkers (C-reactive protein [CRP], erythrocyte sedimentation rate [ESR], D-dimer, ferritin, white blood cell count and interleukin-6) are in the causal pathway between obesity (BMI ≥ 30) and mechanical ventilation or death within 28 days of presentation to care. RESULTS: In the MGH cohort (n = 1202), obesity was associated with greater likelihood of ventilation or death (OR = 1.73; 95% CI = [1.25, 2.41]; P = 0.001) and higher peak CRP (P < 0.001) compared with nonobese patients. The estimated proportion of the association between obesity and ventilation or death mediated by CRP was 0.49 (P < 0.001). Evidence of mediation was more pronounced in patients < 65 years (proportion mediated = 0.52 [P < 0.001] vs 0.44 [P = 0.180]). Findings were more moderate but consistent for peak ESR. Mediation by other inflammatory markers was not supported. Results were replicated in CUIMC/NYP cohort (n = 2626). CONCLUSION: Findings support systemic inflammatory pathways in obesity-associated severe COVID-19 disease, particularly in patients < 65 years, captured by CRP and ESR. Contextualized in clinical trial findings, these results reveal therapeutic opportunity to target systemic inflammatory pathways and monitor interventions in high-risk subgroups and particularly obese patients.

Validating a non-invasive, ALT-based non-alcoholic fatty liver phenotype in the million veteran program.

BACKGROUND & AIMS: Given ongoing challenges in non-invasive non-alcoholic liver disease (NAFLD) diagnosis, we sought to validate an ALT-based NAFLD phenotype using measures readily available in electronic health records (EHRs) and population-based studies by leveraging the clinical and genetic data in the Million Veteran Program (MVP), a multi-ethnic mega-biobank of US Veterans. METHODS: MVP participants with alanine aminotransferases (ALT) >40 units/L for men and >30 units/L for women without other causes of liver disease were compared to controls with normal ALT. Genetic variants spanning eight NAFLD risk or ALT-associated loci (LYPLAL1, GCKR, HSD17B13, TRIB1, PPP1R3B, ERLIN1, TM6SF2, PNPLA3) were tested for NAFLD associations with sensitivity analyses adjusting for metabolic risk factors and alcohol consumption. A manual EHR review assessed performance characteristics of the NAFLD phenotype with imaging and biopsy data as gold standards. Genetic associations with advanced fibrosis were explored using FIB4, NAFLD Fibrosis Score and platelet counts. RESULTS: Among 322,259 MVP participants, 19% met non-invasive criteria for NAFLD. Trans-ethnic meta-analysis replicated associations with previously reported genetic variants in all but LYPLAL1 and GCKR loci (P<6x10-3), without attenuation when adjusted for metabolic risk factors and alcohol consumption. At the previously reported LYPLAL1 locus, the established genetic variant did not appear to be associated with NAFLD, however the regional association plot showed a significant association with NAFLD 279kb downstream. In the EHR validation, the ALT-based NAFLD phenotype yielded a positive predictive value 0.89 and 0.84 for liver biopsy and abdominal imaging, respectively (inter-rater reliability (Cohen's kappa = 0.98)). HSD17B13 and PNPLA3 loci were associated with advanced fibrosis. CONCLUSIONS: We validate a simple, non-invasive ALT-based NAFLD phenotype using EHR data by leveraging previously established NAFLD risk-associated genetic polymorphisms.

Diabetes Prevalence and Its Relationship With Education, Wealth, and BMI in 29 Low- and Middle-Income Countries.

OBJECTIVE: Diabetes is a rapidly growing health problem in low- and middle-income countries (LMICs), but empirical data on its prevalence and relationship to socioeconomic status are scarce. We estimated diabetes prevalence and the subset with undiagnosed diabetes in 29 LMICs and evaluated the relationship of education, household wealth, and BMI with diabetes risk. RESEARCH DESIGN AND METHODS: We pooled individual-level data from 29 nationally representative surveys conducted between 2008 and 2016, totaling 588,574 participants aged ≥25 years. Diabetes prevalence and the subset with undiagnosed diabetes was calculated overall and by country, World Bank income group (WBIG), and geographic region. Multivariable Poisson regression models were used to estimate relative risk (RR). RESULTS: Overall, prevalence of diabetes in 29 LMICs was 7.5% (95% CI 7.1-8.0) and of undiagnosed diabetes 4.9% (4.6-5.3). Diabetes prevalence increased with increasing WBIG: countries with low-income economies (LICs) 6.7% (5.5-8.1), lower-middle-income economies (LMIs) 7.1% (6.6-7.6), and upper-middle-income economies (UMIs) 8.2% (7.5-9.0). Compared with no formal education, greater educational attainment was associated with an increased risk of diabetes across WBIGs, after adjusting for BMI (LICs RR 1.47 [95% CI 1.22-1.78], LMIs 1.14 [1.06-1.23], and UMIs 1.28 [1.02-1.61]). CONCLUSIONS: Among 29 LMICs, diabetes prevalence was substantial and increased with increasing WBIG. In contrast to the association seen in high-income countries, diabetes risk was highest among those with greater educational attainment, independent of BMI. LMICs included in this analysis may be at an advanced stage in the nutrition transition but with no reversal in the socioeconomic gradient of diabetes risk.

Epigenome-Wide Association Study of Incident Type 2 Diabetes in a British Population: EPIC-Norfolk Study.

Epigenetic changes may contribute substantially to risks of diseases of aging. Previous studies reported seven methylation variable positions (MVPs) robustly associated with incident type 2 diabetes mellitus (T2DM). However, their causal roles in T2DM are unclear. In an incident T2DM case-cohort study nested within the population-based European Prospective Investigation into Cancer and Nutrition (EPIC)-Norfolk cohort, we used whole blood DNA collected at baseline, up to 11 years before T2DM onset, to investigate the role of methylation in the etiology of T2DM. We identified 15 novel MVPs with robust associations with incident T2DM and robustly confirmed three MVPs identified previously (near to TXNIP, ABCG1, and SREBF1). All 18 MVPs showed directionally consistent associations with incident and prevalent T2DM in independent studies. Further conditional analyses suggested that the identified epigenetic signals appear related to T2DM via glucose and obesity-related pathways acting before the collection of baseline samples. We integrated genome-wide genetic data to identify methylation-associated quantitative trait loci robustly associated with 16 of the 18 MVPs and found one MVP, cg00574958 at CPT1A, with a possible direct causal role in T2DM. None of the implicated genes were previously highlighted by genetic association studies, suggesting that DNA methylation studies may reveal novel biological mechanisms involved in tissue responses to glycemia.

Cardiovascular Risk Factors and Dehydroepiandrosterone Sulfate Among Latinos in the Boston Puerto Rican Health Study.

Low blood dehydroepiandrosterone sulfate (DHEAS) levels have strong positive associations with stroke and coronary heart disease. However, it is unclear whether DHEAS is independently associated with cardiovascular risk factors. Therefore, we examined the association between cardiovascular risk factors and DHEAS concentration among a high-risk population of Latinos (Puerto Ricans aged 45 to 75 years at baseline) in a cross-sectional analysis of the Boston Puerto Rican Health Study. Of eligible participants, 72% completed baseline interviews and provided blood samples. Complete data were available for 1355 participants. Associations between cardiovascular risk factors (age, sex, total cholesterol, high-density lipid cholesterol, triglycerides, and glucose) and log-transformed DHEAS (µg/dL) were assessed. In robust multivariable regression analyses, DHEAS was significantly inversely associated with age (ß = -12.4; 95% CI: -15.2, -9.7; per 5 years), being female (vs. male) (ß = -46; 95% CI: -55.3, -36.6), and plasma triglyceride concentration (ß = -0.2; 95% CI: -0.3, -0.1; per 10 mg/dL) and was positively associated with total cholesterol and plasma glucose levels (ß = 1.8; 95% CI: 0.6, 3 and ß = 0.2; 95% CI: 0.04, 0.3, respectively, per 10 mg/dL) after adjustment for smoking, alcohol, and physical activity and for postmenopausal hormone use in women. Estimates were unchanged after adjustment for measures of chronic disease and inflammation. Women exhibited a stronger age-related decline in DHEAS and a positive association with glucose in contrast to findings among men (P interaction < 0.05). In conclusion, in this large study of Latinos with a heavy cardiovascular risk factor burden, we observed significant associations between cardiovascular disease (CVD) risk factors and DHEAS, with variations by sex. These findings improve our understanding of the role DHEAS may play in CVD etiology.

Association of Methylation Signals With Incident Coronary Heart Disease in an Epigenome-Wide Assessment of Circulating Tumor Necrosis Factor α.

Importance: Tumor necrosis factor α (TNF-α) is a proinflammatory cytokine with manifold consequences for mammalian pathophysiology, including cardiovascular disease. A deeper understanding of TNF-α biology may enhance treatment precision. Objective: To conduct an epigenome-wide analysis of blood-derived DNA methylation and TNF-α levels and to assess the clinical relevance of findings. Design, Setting, and Participants: This meta-analysis assessed epigenome-wide associations in circulating TNF-α concentrations from 5 cohort studies and 1 interventional trial, with replication in 3 additional cohort studies. Follow-up analyses investigated associations of identified methylation loci with gene expression and incident coronary heart disease; this meta-analysis included 11 461 participants who experienced 1895 coronary events. Exposures: Circulating TNF-α concentration. Main Outcomes and Measures: DNA methylation at approximately 450 000 loci, neighboring DNA sequence variation, gene expression, and incident coronary heart disease. Results: The discovery cohort included 4794 participants, and the replication study included 816 participants (overall mean [SD] age, 60.7 [8.5] years). In the discovery stage, circulating TNF-α levels were associated with methylation of 7 cytosine-phosphate-guanine (CpG) sites, 3 of which were located in or near DTX3L-PARP9 at cg00959259 (ß [SE] = -0.01 [0.003]; P = 7.36 × 10-8), cg08122652 (ß [SE] = -0.008 [0.002]; P = 2.24 × 10-7), and cg22930808(ß [SE] = -0.01 [0.002]; P = 6.92 × 10-8); NLRC5 at cg16411857 (ß [SE] = -0.01 [0.002]; P = 2.14 × 10-13) and cg07839457 (ß [SE] = -0.02 [0.003]; P = 6.31 × 10-10); or ABO, at cg13683939 (ß [SE] = 0.04 [0.008]; P = 1.42 × 10-7) and cg24267699 (ß [SE] = -0.009 [0.002]; P = 1.67 × 10-7), after accounting for multiple testing. Of these, negative associations between TNF-α concentration and methylation of 2 loci in NLRC5 and 1 in DTX3L-14 PARP9 were replicated. Replicated TNF-α-linked CpG sites were associated with 9% to 19% decreased risk of incident coronary heart disease per 10% higher methylation per CpG site (cg16411857: hazard ratio [HR], 0.86; 95% CI, 0.78-1.95; P = .003; cg07839457: HR, 0.89; 95% CI, 0.80-0.94; P = 3.1 × 10-5; cg00959259: HR, 0.91; 95% CI, 0.84-0.97; P = .002; cg08122652: HR, 0.81; 95% CI, 0.74-0.89; P = 2.0 × 10-5). Conclusions and Relevance: We identified and replicated novel epigenetic correlates of circulating TNF-α concentration in blood samples and linked these loci to coronary heart disease risk, opening opportunities for validation and therapeutic applications.

Ambient air pollution, adipokines, and glucose homeostasis: The Framingham Heart Study.

OBJECTIVE: To examine associations of proximity to major roadways, sustained exposure to fine particulate matter (PM2.5), and acute exposure to ambient air pollutants with adipokines and measures of glucose homeostasis among participants living in the northeastern United States. METHODS: We included 5958 participants from the Framingham Offspring cohort examination cycle 7 (1998-2001) and 8 (2005-2008) and Third Generation cohort examination cycle 1 (2002-2005) and 2 (2008-2011), who did not have type 2 diabetes at the time of examination visit. We calculated 2003 annual average PM2.5 at participants' home address, residential distance to the nearest major roadway, and daily PM2.5, black carbon (BC), sulfate, nitrogen oxides (NOx), and ozone concentrations. We used linear mixed effects models for fasting glucose, insulin, and homeostasis model assessment of insulin resistance (HOMA-IR) which were measured up to twice, and used linear regression models for adiponectin, resistin, leptin, and hemoglobin A1c (HbA1c) which were measured only once, adjusting for demographics, socioeconomic position, lifestyle, time, and seasonality. RESULTS: The mean age was 51years and 55% were women. Participants who lived 64m (25th percentile) from a major roadway had 0.28% (95% CI: 0.05%, 0.51%) higher fasting plasma glucose than participants who lived 413m (75th percentile) away, and the association appeared to be driven by participants who lived within 50m from a major roadway. Higher exposures to 3- to 7-day moving averages of BC and NOx were associated with higher glucose whereas the associations for ozone were negative. The associations otherwise were generally null and did not differ by median age, sex, educational attainment, obesity status, or prediabetes status. CONCLUSIONS: Living closer to a major roadway or acute exposure to traffic-related air pollutants were associated with dysregulated glucose homeostasis but not with adipokines among participants from the Framingham Offspring and Third Generation cohorts.

Determinants of Smoking and Quitting in HIV-Infected Individuals.

BACKGROUND: Cigarette smoking is widespread among HIV-infected patients, who confront increased risk of smoking-related co-morbidities. The effects of HIV infection and HIV-related variables on smoking and smoking cessation are incompletely understood. We investigated the correlates of smoking and quitting in an HIV-infected cohort using a validated natural language processor to determine smoking status. METHOD: We developed and validated an algorithm using natural language processing (NLP) to ascertain smoking status from electronic health record data. The algorithm was applied to records for a cohort of 3487 HIV-infected from a large health care system in Boston, USA, and 9446 uninfected control patients matched 3:1 on age, gender, race and clinical encounters. NLP was used to identify and classify smoking-related portions of free-text notes. These classifications were combined into patient-year smoking status and used to classify patients as ever versus never smokers and current smokers versus non-smokers. Generalized linear models were used to assess associations of HIV with 3 outcomes, ever smoking, current smoking, and current smoking in analyses limited to ever smokers (persistent smoking), while adjusting for demographics, cardiovascular risk factors, and psychiatric illness. Analyses were repeated within the HIV cohort, with the addition of CD4 cell count and HIV viral load to assess associations of these HIV-related factors with the smoking outcomes. RESULTS: Using the natural language processing algorithm to assign annual smoking status yielded sensitivity of 92.4, specificity of 86.2, and AUC of 0.89 (95% confidence interval [CI] 0.88-0.91). Ever and current smoking were more common in HIV-infected patients than controls (54% vs. 44% and 42% vs. 30%, respectively, both P<0.001). In multivariate models HIV was independently associated with ever smoking (adjusted rate ratio [ARR] 1.18, 95% CI 1.13-1.24, P <0.001), current smoking (ARR 1.33, 95% CI 1.25-1.40, P<0.001), and persistent smoking (ARR 1.11, 95% CI 1.07-1.15, P<0.001). Within the HIV cohort, having a detectable HIV RNA was significantly associated with all three smoking outcomes. CONCLUSIONS: HIV was independently associated with both smoking and not quitting smoking, using a novel algorithm to ascertain smoking status from electronic health record data and accounting for multiple confounding clinical factors. Further research is needed to identify HIV-related barriers to smoking cessation and develop aggressive interventions specific to HIV-infected patients.

General Framework for Meta-Analysis of Haplotype Association Tests.

For complex traits, most associated single nucleotide variants (SNV) discovered to date have a small effect, and detection of association is only possible with large sample sizes. Because of patient confidentiality concerns, it is often not possible to pool genetic data from multiple cohorts, and meta-analysis has emerged as the method of choice to combine results from multiple studies. Many meta-analysis methods are available for single SNV analyses. As new approaches allow the capture of low frequency and rare genetic variation, it is of interest to jointly consider multiple variants to improve power. However, for the analysis of haplotypes formed by multiple SNVs, meta-analysis remains a challenge, because different haplotypes may be observed across studies. We propose a two-stage meta-analysis approach to combine haplotype analysis results. In the first stage, each cohort estimate haplotype effect sizes in a regression framework, accounting for relatedness among observations if appropriate. For the second stage, we use a multivariate generalized least square meta-analysis approach to combine haplotype effect estimates from multiple cohorts. Haplotype-specific association tests and a global test of independence between haplotypes and traits are obtained within our framework. We demonstrate through simulation studies that we control the type-I error rate, and our approach is more powerful than inverse variance weighted meta-analysis of single SNV analysis when haplotype effects are present. We replicate a published haplotype association between fasting glucose-associated locus (G6PC2) and fasting glucose in seven studies from the Cohorts for Heart and Aging Research in Genomic Epidemiology Consortium and we provide more precise haplotype effect estimates.

Transancestral fine-mapping of four type 2 diabetes susceptibility loci highlights potential causal regulatory mechanisms.

To gain insight into potential regulatory mechanisms through which the effects of variants at four established type 2 diabetes (T2D) susceptibility loci (CDKAL1, CDKN2A-B, IGF2BP2 and KCNQ1) are mediated, we undertook transancestral fine-mapping in 22 086 cases and 42 539 controls of East Asian, European, South Asian, African American and Mexican American descent. Through high-density imputation and conditional analyses, we identified seven distinct association signals at these four loci, each with allelic effects on T2D susceptibility that were homogenous across ancestry groups. By leveraging differences in the structure of linkage disequilibrium between diverse populations, and increased sample size, we localised the variants most likely to drive each distinct association signal. We demonstrated that integration of these genetic fine-mapping data with genomic annotation can highlight potential causal regulatory elements in T2D-relevant tissues. These analyses provide insight into the mechanisms through which T2D association signals are mediated, and suggest future routes to understanding the biology of specific disease susceptibility loci.

Type 2 Diabetes Genetic Predisposition, Obesity, and All-Cause Mortality Risk in the U.S.: A Multiethnic Analysis.

OBJECTIVE: Type 2 diabetes (T2D) is associated with increased mortality in ethnically diverse populations, although the extent to which this association is genetically determined is unknown. We sought to determine whether T2D-related genetic variants predicted all-cause mortality, even after accounting for BMI, in the Third National Health and Nutrition Examination Survey. RESEARCH DESIGN AND METHODS: We modeled mortality risk using a genetic risk score (GRS) from a weighted sum of risk alleles at 38 T2D-related single nucleotide polymorphisms. In age-, sex-, and BMI-adjusted logistic regression models, accounting for the complex survey design, we tested the association with mortality in 6,501 participants. We repeated the analysis within ethnicities (2,528 non-Hispanic white [NHW], 1,979 non-Hispanic black [NHB], and 1,994 Mexican American [MA]) and within BMI categories (<25, 25-30, and ≥30 kg/m(2)). Significance was set at P < 0.05. RESULTS: Over 17 years, 1,556 participants died. GRS was associated with mortality risk (OR 1.04 [95% CI 1.00-1.07] per T2D-associated risk allele, P = 0.05). Within ethnicities, GRS was positively associated with mortality risk in NHW and NHB, but not in MA (0.95 [0.90-1.01], P = 0.07). The negative trend in MA was largely driven by those with BMI <25 kg/m(2) (0.91 [0.82-1.00]). In NHW, the positive association was strongest among those with BMI ≥30 kg/m(2) (1.07 [1.02-1.12]). CONCLUSIONS: In the U.S., a higher T2D genetic risk was associated with increased mortality risk, especially among obese NHW. The underlying genetic basis for mortality likely involves complex interactions with factors related to ethnicity, T2D, and body weight.

Integromic analysis of genetic variation and gene expression identifies networks for cardiovascular disease phenotypes.

BACKGROUND: Cardiovascular disease (CVD) reflects a highly coordinated complex of traits. Although genome-wide association studies have reported numerous single nucleotide polymorphisms (SNPs) to be associated with CVD, the role of most of these variants in disease processes remains unknown. METHODS AND RESULTS: We built a CVD network using 1512 SNPs associated with 21 CVD traits in genome-wide association studies (at P≤5×10(-8)) and cross-linked different traits by virtue of their shared SNP associations. We then explored whole blood gene expression in relation to these SNPs in 5257 participants in the Framingham Heart Study. At a false discovery rate <0.05, we identified 370 cis-expression quantitative trait loci (eQTLs; SNPs associated with altered expression of nearby genes) and 44 trans-eQTLs (SNPs associated with altered expression of remote genes). The eQTL network revealed 13 CVD-related modules. Searching for association of eQTL genes with CVD risk factors (lipids, blood pressure, fasting blood glucose, and body mass index) in the same individuals, we found examples in which the expression of eQTL genes was significantly associated with these CVD phenotypes. In addition, mediation tests suggested that a subset of SNPs previously associated with CVD phenotypes in genome-wide association studies may exert their function by altering expression of eQTL genes (eg, LDLR and PCSK7), which in turn may promote interindividual variation in phenotypes. CONCLUSIONS: Using a network approach to analyze CVD traits, we identified complex networks of SNP-phenotype and SNP-transcript connections. Integrating the CVD network with phenotypic data, we identified biological pathways that may provide insights into potential drug targets for treatment or prevention of CVD.

Electronic medical record cancer incidence over six years comparing new users of glargine with new users of NPH insulin.

BACKGROUND: Recent studies suggested that insulin glargine use could be associated with increased risk of cancer. We compared the incidence of cancer in new users of glargine versus new users of NPH in a longitudinal clinical cohort with diabetes for up to 6 years. METHODS AND FINDINGS: From all patients who had been regularly followed at Massachusetts General Hospital from 1/01/2005 to 12/31/2010, 3,680 patients who had a medication record for glargine or NPH usage were obtained from the electronic medical record (EMR). From those we selected 539 new glargine users (age: 60.1±13.6 years, BMI: 32.7±7.5 kg/m2) and 343 new NPH users (61.5±14.1 years, 32.7±8.3 kg/m2) who had no prevalent cancer during 19 months prior to glargine or NPH initiation. All incident cancer cases were ascertained from the EMR requiring at least 2 ICD-9 codes within a 2 month period. Insulin exposure time and cumulative dose were validated. The statistical analysis compared the rates of cancer in new glargine vs. new NPH users while on treatment, adjusted for the propensity to receive one or the other insulin. There were 26 and 28 new cancer cases in new glargine and new NPH users for 1559 and 1126 person-years follow-up, respectively. There were no differences in the propensity-adjusted clinical characteristics between groups. The adjusted hazard ratio for the cancer incidence comparing glargine vs. NPH use was 0.65 (95% CI: 0.36-1.19). CONCLUSIONS: Insulin glargine is not associated with development of cancers when compared with NPH in this longitudinal and carefully retrieved EMR data.

Association of levels of fasting glucose and insulin with rare variants at the chromosome 11p11.2-MADD locus: Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium Targeted Sequencing Study.

BACKGROUND: Common variation at the 11p11.2 locus, encompassing MADD, ACP2, NR1H3, MYBPC3, and SPI1, has been associated in genome-wide association studies with fasting glucose and insulin (FI). In the Cohorts for Heart and Aging Research in Genomic Epidemiology Targeted Sequencing Study, we sequenced 5 gene regions at 11p11.2 to identify rare, potentially functional variants influencing fasting glucose or FI levels. METHODS AND RESULTS: Sequencing (mean depth, 38×) across 16.1 kb in 3566 individuals without diabetes mellitus identified 653 variants, 79.9% of which were rare (minor allele frequency <1%) and novel. We analyzed rare variants in 5 gene regions with FI or fasting glucose using the sequence kernel association test. At NR1H3, 53 rare variants were jointly associated with FI (P=2.73×10(-3)); of these, 7 were predicted to have regulatory function and showed association with FI (P=1.28×10(-3)). Conditioning on 2 previously associated variants at MADD (rs7944584, rs10838687) did not attenuate this association, suggesting that there are >2 independent signals at 11p11.2. One predicted regulatory variant, chr11:47227430 (hg18; minor allele frequency=0.00068), contributed 20.6% to the overall sequence kernel association test score at NR1H3, lies in intron 2 of NR1H3, and is a predicted binding site for forkhead box A1 (FOXA1), a transcription factor associated with insulin regulation. In human HepG2 hepatoma cells, the rare chr11:47227430 A allele disrupted FOXA1 binding and reduced FOXA1-dependent transcriptional activity. CONCLUSIONS: Sequencing at 11p11.2-NR1H3 identified rare variation associated with FI. One variant, chr11:47227430, seems to be functional, with the rare A allele reducing transcription factor FOXA1 binding and FOXA1-dependent transcriptional activity.

Magnesium intake is inversely associated with coronary artery calcification: the Framingham Heart Study.

OBJECTIVES: The aim of this study was to examine whether magnesium intake is associated with coronary artery calcification (CAC) and abdominal aortic calcification (AAC). BACKGROUND: Animal and cell studies suggest that magnesium may prevent calcification within atherosclerotic plaques underlying cardiovascular disease. Little is known about the association of magnesium intake and atherosclerotic calcification in humans. METHODS: We examined cross-sectional associations of self-reported total (dietary and supplemental) magnesium intake estimated by food frequency questionnaire with CAC and AAC in participants of the Framingham Heart Study who were free of cardiovascular disease and underwent Multi-Detector Computed Tomography (MDCT) of the heart and abdomen (n = 2,695; age: 53 ± 11 years), using multivariate-adjusted Tobit regression. CAC and AAC were quantified using modified Agatston scores (AS). Models were adjusted for age, sex, body mass index, smoking status, systolic blood pressure, fasting insulin, total-to-high-density lipoprotein cholesterol ratio, use of hormone replacement therapy (women only), menopausal status (women only), treatment for hyperlipidemia, hypertension, cardiovascular disease prevention, or diabetes, as well as self-reported intake of calcium, vitamins D and K, saturated fat, fiber, alcohol, and energy. Secondary analyses included logistic regressions of CAC and AAC outcomes as cut-points (AS >0 and AS ≥90th percentile for age and sex), as well as sex-stratified analyses. RESULTS: In fully adjusted models, a 50-mg/day increment in self-reported total magnesium intake was associated with 22% lower CAC (p < 0.001) and 12% lower AAC (p = 0.07). Consistent with these observations, the odds of having any CAC were 58% lower (p trend: <0.001) and any AAC were 34% lower (p trend: 0.01), in those with the highest compared to those with the lowest magnesium intake. Stronger inverse associations were observed in women than in men. CONCLUSIONS: In community-dwelling participants free of cardiovascular disease, self-reported magnesium intake was inversely associated with arterial calcification, which may play a contributing role in magnesium's protective associations in stroke and fatal coronary heart disease.