DNA methylation near MAD1L1, KDM2B, and SOCS3 mediates the effect of socioeconomic status on elevated body mass index in African American adults.

Obesity and poverty disproportionally affect African American persons. Epigenetic mechanisms could partially explain the association between socioeconomic disadvantage and body mass index (BMI). We examined the extent to which epigenetic mechanisms mediate the effect of socioeconomic status (SES) on BMI. Using data from African American adults from the Atherosclerosis Risk in Communities (ARIC) Study (n = 2664, mean age = 57 years), education, income, and occupation were used to create a composite SES score at visit 1 (1987-1989). We conducted two methylation-wide association analyses to identify associations between SES (visit 1), BMI and cytosine-phosphate-guanine (CpG) sites measured at a subsequent visit (1990-1995). We then utilized structural equation modeling (SEM) to test whether identified sites mediated the association between earlier SES and BMI in sex-stratified models adjusted for demographic and risk factor covariates. Independent replication and meta-analyses were conducted using the Jackson Heart Study (JHS, n = 874, mean age 51 years, 2000-2004). Three CpG sites near MAD1L1, KDM2B, and SOCS3 (cg05095590, cg1370865, and cg18181703) were suggestively associated (P-value < 1.3×10-5) in ARIC and at array-wide significance (P-value < 1.3×10-7) in a combined meta-analysis of ARIC with JHS. SEM of these three sites revealed significant indirect effects in females (P-value < 5.8×10-3), each mediating 7%-20% of the total effect of SES on BMI. Nominally significant indirect effects were observed for two sites near MAD1L1 and KDM2B in males (P-value < 3.4×10-2), mediating -17 and -22% of the SES-BMI effect. These results provide further evidence that epigenetic modifications may be a potential pathway through which SES may "get under the skin" and contribute to downstream health disparities.

Predicted gene expression in ancestrally diverse populations leads to discovery of susceptibility loci for lifestyle and cardiometabolic traits.

One mechanism by which genetic factors influence complex traits and diseases is altering gene expression. Direct measurement of gene expression in relevant tissues is rarely tenable; however, genetically regulated gene expression (GReX) can be estimated using prediction models derived from large multi-omic datasets. These approaches have led to the discovery of many gene-trait associations, but whether models derived from predominantly European ancestry (EA) reference panels can map novel associations in ancestrally diverse populations remains unclear. We applied PrediXcan to impute GReX in 51,520 ancestrally diverse Population Architecture using Genomics and Epidemiology (PAGE) participants (35% African American, 45% Hispanic/Latino, 10% Asian, and 7% Hawaiian) across 25 key cardiometabolic traits and relevant tissues to identify 102 novel associations. We then compared associations in PAGE to those in a random subset of 50,000 White British participants from UK Biobank (UKBB50k) for height and body mass index (BMI). We identified 517 associations across 47 tissues in PAGE but not UKBB50k, demonstrating the importance of diverse samples in identifying trait-associated GReX. We observed that variants used in PrediXcan models were either more or less differentiated across continental-level populations than matched-control variants depending on the specific population reflecting sampling bias. Additionally, variants from identified genes specific to either PAGE or UKBB50k analyses were more ancestrally differentiated than those in genes detected in both analyses, underlining the value of population-specific discoveries. This suggests that while EA-derived transcriptome imputation models can identify new associations in non-EA populations, models derived from closely matched reference panels may yield further insights. Our findings call for more diversity in reference datasets of tissue-specific gene expression.

Genetic analyses of diverse populations improves discovery for complex traits.

Genome-wide association studies (GWAS) have laid the foundation for investigations into the biology of complex traits, drug development and clinical guidelines. However, the majority of discovery efforts are based on data from populations of European ancestry1-3. In light of the differential genetic architecture that is known to exist between populations, bias in representation can exacerbate existing disease and healthcare disparities. Critical variants may be missed if they have a low frequency or are completely absent in European populations, especially as the field shifts its attention towards rare variants, which are more likely to be population-specific4-10. Additionally, effect sizes and their derived risk prediction scores derived in one population may not accurately extrapolate to other populations11,12. Here we demonstrate the value of diverse, multi-ethnic participants in large-scale genomic studies. The Population Architecture using Genomics and Epidemiology (PAGE) study conducted a GWAS of 26 clinical and behavioural phenotypes in 49,839 non-European individuals. Using strategies tailored for analysis of multi-ethnic and admixed populations, we describe a framework for analysing diverse populations, identify 27 novel loci and 38 secondary signals at known loci, as well as replicate 1,444 GWAS catalogue associations across these traits. Our data show evidence of effect-size heterogeneity across ancestries for published GWAS associations, substantial benefits for fine-mapping using diverse cohorts and insights into clinical implications. In the United States-where minority populations have a disproportionately higher burden of chronic conditions13-the lack of representation of diverse populations in genetic research will result in inequitable access to precision medicine for those with the highest burden of disease. We strongly advocate for continued, large genome-wide efforts in diverse populations to maximize genetic discovery and reduce health disparities.

Exome sequencing of 1190 non-syndromic clubfoot cases reveals HOXD12 as a novel disease gene.

BACKGROUND: Clubfoot, presenting as a rigid inward and downward turning of the foot, is one of the most common congenital musculoskeletal anomalies. The aetiology of clubfoot is poorly understood and variants in known clubfoot disease genes account for only a small portion of the heritability. METHODS: Exome sequence data were generated from 1190 non-syndromic clubfoot cases and their family members from multiple ethnicities. Ultra-rare variant burden analysis was performed comparing 857 unrelated clubfoot cases with European ancestry with two independent ethnicity-matched control groups (1043 in-house and 56 885 gnomAD controls). Additional variants in prioritised genes were identified in a larger cohort, including probands with non-European ancestry. Segregation analysis was performed in multiplex families when available. RESULTS: Rare variants in 29 genes were enriched in clubfoot cases, including PITX1 (a known clubfoot disease gene), HOXD12, COL12A1, COL9A3 and LMX1B. In addition, rare variants in posterior HOX genes (HOX9-13) were enriched overall in clubfoot cases. In total, variants in these genes were present in 8.4% (100/1190) of clubfoot cases with both European and non-European ancestry. Among these, 3 are de novo and 22 show variable penetrance, including 4 HOXD12 variants that segregate with clubfoot. CONCLUSION: We report HOXD12 as a novel clubfoot disease gene and demonstrate a phenotypic expansion of known disease genes (myopathy gene COL12A1, Ehlers-Danlos syndrome gene COL9A3 and nail-patella syndrome gene LMX1B) to include isolated clubfoot.

Discovery and fine-mapping of height loci via high-density imputation of GWASs in individuals of African ancestry.

Although many loci have been associated with height in European ancestry populations, very few have been identified in African ancestry individuals. Furthermore, many of the known loci have yet to be generalized to and fine-mapped within a large-scale African ancestry sample. We performed sex-combined and sex-stratified meta-analyses in up to 52,764 individuals with height and genome-wide genotyping data from the African Ancestry Anthropometry Genetics Consortium (AAAGC). We additionally combined our African ancestry meta-analysis results with published European genome-wide association study (GWAS) data. In the African ancestry analyses, we identified three novel loci (SLC4A3, NCOA2, ECD/FAM149B1) in sex-combined results and two loci (CRB1, KLF6) in women only. In the African plus European sex-combined GWAS, we identified an additional three novel loci (RCCD1, G6PC3, CEP95) which were equally driven by AAAGC and European results. Among 39 genome-wide significant signals at known loci, conditioning index SNPs from European studies identified 20 secondary signals. Two of the 20 new secondary signals and none of the 8 novel loci had minor allele frequencies (MAF) < 5%. Of 802 known European height signals, 643 displayed directionally consistent associations with height, of which 205 were nominally significant (p < 0.05) in the African ancestry sex-combined sample. Furthermore, 148 of 241 loci contained ≤20 variants in the credible sets that jointly account for 99% of the posterior probability of driving the associations. In summary, trans-ethnic meta-analyses revealed novel signals and further improved fine-mapping of putative causal variants in loci shared between African and European ancestry populations.

Pan-ancestry exome-wide association analyses of COVID-19 outcomes in 586,157 individuals.

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) causes coronavirus disease 2019 (COVID-19), a respiratory illness that can result in hospitalization or death. We used exome sequence data to investigate associations between rare genetic variants and seven COVID-19 outcomes in 586,157 individuals, including 20,952 with COVID-19. After accounting for multiple testing, we did not identify any clear associations with rare variants either exome wide or when specifically focusing on (1) 13 interferon pathway genes in which rare deleterious variants have been reported in individuals with severe COVID-19, (2) 281 genes located in susceptibility loci identified by the COVID-19 Host Genetics Initiative, or (3) 32 additional genes of immunologic relevance and/or therapeutic potential. Our analyses indicate there are no significant associations with rare protein-coding variants with detectable effect sizes at our current sample sizes. Analyses will be updated as additional data become available, and results are publicly available through the Regeneron Genetics Center COVID-19 Results Browser.

Genome-wide association study of body fat distribution traits in Hispanics/Latinos from the HCHS/SOL.

Central obesity is a leading health concern with a great burden carried by ethnic minority populations, especially Hispanics/Latinos. Genetic factors contribute to the obesity burden overall and to inter-population differences. We aimed to identify the loci associated with central adiposity measured as waist-to-hip ratio (WHR), waist circumference (WC) and hip circumference (HIP) adjusted for body mass index (adjBMI) by using the Hispanic Community Health Study/Study of Latinos (HCHS/SOL); determine if differences in associations differ by background group within HCHS/SOL and determine whether previously reported associations generalize to HCHS/SOL. Our analyses included 7472 women and 5200 men of mainland (Mexican, Central and South American) and Caribbean (Puerto Rican, Cuban and Dominican) background residing in the USA. We performed genome-wide association analyses stratified and combined across sexes using linear mixed-model regression. We identified 16 variants for waist-to-hip ratio adjusted for body mass index (WHRadjBMI), 22 for waist circumference adjusted for body mass index (WCadjBMI) and 28 for hip circumference adjusted for body mass index (HIPadjBMI), which reached suggestive significance (P < 1 × 10-6). Many loci exhibited differences in strength of associations by ethnic background and sex. We brought a total of 66 variants forward for validation in cohorts (N = 34 161) with participants of Hispanic/Latino, African and European descent. We confirmed four novel loci (P < 0.05 and consistent direction of effect, and P < 5 × 10-8 after meta-analysis), including two for WHRadjBMI (rs13301996, rs79478137); one for WCadjBMI (rs3168072) and one for HIPadjBMI (rs28692724). Also, we generalized previously reported associations to HCHS/SOL, (8 for WHRadjBMI, 10 for WCadjBMI and 12 for HIPadjBMI). Our study highlights the importance of large-scale genomic studies in ancestrally diverse Hispanic/Latino populations for identifying and characterizing central obesity susceptibility that may be ancestry-specific.

Assessing efficiency of fine-mapping obesity-associated variants through leveraging ancestry architecture and functional annotation using PAGE and UKBB cohorts.

Inadequate representation of non-European ancestry populations in genome-wide association studies (GWAS) has limited opportunities to isolate functional variants. Fine-mapping in multi-ancestry populations should improve the efficiency of prioritizing variants for functional interrogation. To evaluate this hypothesis, we leveraged ancestry architecture to perform comparative GWAS and fine-mapping of obesity-related phenotypes in European ancestry populations from the UK Biobank (UKBB) and multi-ancestry samples from the Population Architecture for Genetic Epidemiology (PAGE) consortium with comparable sample sizes. In the investigated regions with genome-wide significant associations for obesity-related traits, fine-mapping in our ancestrally diverse sample led to 95% and 99% credible sets (CS) with fewer variants than in the European ancestry sample. Lead fine-mapped variants in PAGE regions had higher average coding scores, and higher average posterior probabilities for causality compared to UKBB. Importantly, 99% CS in PAGE loci contained strong expression quantitative trait loci (eQTLs) in adipose tissues or harbored more variants in tighter linkage disequilibrium (LD) with eQTLs. Leveraging ancestrally diverse populations with heterogeneous ancestry architectures, coupled with functional annotation, increased fine-mapping efficiency and performance, and reduced the set of candidate variants for consideration for future functional studies. Significant overlap in genetic causal variants across populations suggests generalizability of genetic mechanisms underpinning obesity-related traits across populations.

Rare and low-frequency coding variants alter human adult height.

Height is a highly heritable, classic polygenic trait with approximately 700 common associated variants identified through genome-wide association studies so far. Here, we report 83 height-associated coding variants with lower minor-allele frequencies (in the range of 0.1-4.8%) and effects of up to 2 centimetres per allele (such as those in IHH, STC2, AR and CRISPLD2), greater than ten times the average effect of common variants. In functional follow-up studies, rare height-increasing alleles of STC2 (giving an increase of 1-2 centimetres per allele) compromised proteolytic inhibition of PAPP-A and increased cleavage of IGFBP-4 in vitro, resulting in higher bioavailability of insulin-like growth factors. These 83 height-associated variants overlap genes that are mutated in monogenic growth disorders and highlight new biological candidates (such as ADAMTS3, IL11RA and NOX4) and pathways (such as proteoglycan and glycosaminoglycan synthesis) involved in growth. Our results demonstrate that sufficiently large sample sizes can uncover rare and low-frequency variants of moderate-to-large effect associated with polygenic human phenotypes, and that these variants implicate relevant genes and pathways.

Pulmonary Function and Blood DNA Methylation: A Multiancestry Epigenome-Wide Association Meta-analysis.

Rationale: Methylation integrates factors present at birth and modifiable across the lifespan that can influence pulmonary function. Studies are limited in scope and replication. Objectives: To conduct large-scale epigenome-wide meta-analyses of blood DNA methylation and pulmonary function. Methods: Twelve cohorts analyzed associations of methylation at cytosine-phosphate-guanine probes (CpGs), using Illumina 450K or EPIC/850K arrays, with FEV1, FVC, and FEV1/FVC. We performed multiancestry epigenome-wide meta-analyses (total of 17,503 individuals; 14,761 European, 2,549 African, and 193 Hispanic/Latino ancestries) and interpreted results using integrative epigenomics. Measurements and Main Results: We identified 1,267 CpGs (1,042 genes) differentially methylated (false discovery rate, <0.025) in relation to FEV1, FVC, or FEV1/FVC, including 1,240 novel and 73 also related to chronic obstructive pulmonary disease (1,787 cases). We found 294 CpGs unique to European or African ancestry and 395 CpGs unique to never or ever smokers. The majority of significant CpGs correlated with nearby gene expression in blood. Findings were enriched in key regulatory elements for gene function, including accessible chromatin elements, in both blood and lung. Sixty-nine implicated genes are targets of investigational or approved drugs. One example novel gene highlighted by integrative epigenomic and druggable target analysis is TNFRSF4. Mendelian randomization and colocalization analyses suggest that epigenome-wide association study signals capture causal regulatory genomic loci. Conclusions: We identified numerous novel loci differentially methylated in relation to pulmonary function; few were detected in large genome-wide association studies. Integrative analyses highlight functional relevance and potential therapeutic targets. This comprehensive discovery of potentially modifiable, novel lung function loci expands knowledge gained from genetic studies, providing insights into lung pathogenesis.

Infant Growth Trajectories and Lipid Levels in Adolescence: Evidence From a Chilean Infancy Cohort.

Growth in early infancy is hypothesized to affect chronic disease risk factors later in life. To date, most reports draw on European-ancestry cohorts with few repeated observations in early infancy. We investigated the association between infant growth before 6 months and lipid levels in adolescents in a Hispanic/Latino cohort. We characterized infant growth from birth to 5 months in male (n = 311) and female (n = 285) infants from the Santiago Longitudinal Study (1991-1996) using 3 metrics: weight (kg), length (cm), and weight-for-length (g/cm). Superimposition by translation and rotation (SITAR) and latent growth mixture models (LGMMs) were used to estimate the association between infant growth characteristics and lipid levels at age 17 years. We found a positive relationship between the SITAR length velocity parameter before 6 months of age and high-density lipoprotein cholesterol levels in adolescence (11.5, 95% confidence interval; 3.4, 19.5), indicating higher high-density lipoprotein cholesterol levels occurring with faster length growth. The strongest associations from the LGMMs were between higher low-density lipoprotein cholesterol and slower weight-for-length growth, following a pattern of associations between slower growth and adverse lipid profiles. Further research in this window of time can confirm the association between early infant growth as an exposure and adolescent cardiovascular disease risk factors.

Genetic determinants of metabolic biomarkers and their associations with cardiometabolic traits in Hispanic/Latino adolescents.

BACKGROUND: Metabolic regulation plays a significant role in energy homeostasis, and adolescence is a crucial life stage for the development of cardiometabolic disease (CMD). This study aims to investigate the genetic determinants of metabolic biomarkers-adiponectin, leptin, ghrelin, and orexin-and their associations with CMD risk factors. METHODS: We characterized the genetic determinants of the biomarkers among Hispanic/Latino adolescents of the Santiago Longitudinal Study (SLS) and identified the cumulative effects of genetic variants on adiponectin and leptin using biomarker polygenic risk scores (PRS). We further investigated the direct and indirect effect of the biomarker PRS on downstream body fat percent (BF%) and glycemic traits using structural equation modeling. RESULTS: We identified putatively novel genetic variants associated with the metabolic biomarkers. A substantial amount of biomarker variance was explained by SLS-specific PRS, and the prediction was improved by including the putatively novel loci. Fasting blood insulin and insulin resistance were associated with PRS for adiponectin, leptin, and ghrelin, and BF% was associated with PRS for adiponectin and leptin. We found evidence of substantial mediation of these associations by the biomarker levels. CONCLUSIONS: The genetic underpinnings of metabolic biomarkers can affect the early development of CMD, partly mediated by the biomarkers. IMPACT: This study characterized the genetic underpinnings of four metabolic hormones and investigated their potential influence on adiposity and insulin biology among Hispanic/Latino adolescents. Fasting blood insulin and insulin resistance were associated with polygenic risk score (PRS) for adiponectin, leptin, and ghrelin, with evidence of some degree of mediation by the biomarker levels. Body fat percent (BF%) was also associated with PRS for adiponectin and leptin. This provides important insight on biological mechanisms underlying early metabolic dysfunction and reveals candidates for prevention efforts. Our findings also highlight the importance of ancestrally diverse populations to facilitate valid studies of the genetic architecture of metabolic biomarker levels.

Gaseous air pollutants and DNA methylation in a methylome-wide association study of an ethnically and environmentally diverse population of U.S. adults.

Epigenetic mechanisms may underlie air pollution-health outcome associations. We estimated gaseous air pollutant-DNA methylation (DNAm) associations using twelve subpopulations within Women's Health Initiative (WHI) and Atherosclerosis Risk in Communities (ARIC) cohorts (n = 8397; mean age 61.3 years; 83% female; 46% African-American, 46% European-American, 8% Hispanic/Latino). We used geocoded participant address-specific mean ambient carbon monoxide (CO), nitrogen oxides (NO2; NOx), ozone (O3), and sulfur dioxide (SO2) concentrations estimated over the 2-, 7-, 28-, and 365-day periods before collection of blood samples used to generate Illumina 450 k array leukocyte DNAm measurements. We estimated methylome-wide, subpopulation- and race/ethnicity-stratified pollutant-DNAm associations in multi-level, linear mixed-effects models adjusted for sociodemographic, behavioral, meteorological, and technical covariates. We combined stratum-specific estimates in inverse variance-weighted meta-analyses and characterized significant associations (false discovery rate; FDR<0.05) at Cytosine-phosphate-Guanine (CpG) sites without among-strata heterogeneity (PCochran's Q > 0.05). We attempted replication in the Cooperative Health Research in Region of Augsburg (KORA) study and Normative Aging Study (NAS). We observed a -0.3 (95% CI: -0.4, -0.2) unit decrease in percent DNAm per interquartile range (IQR, 7.3 ppb) increase in 28-day mean NO2 concentration at cg01885635 (chromosome 3; regulatory region 290 bp upstream from ZNF621; FDR = 0.03). At intragenic sites cg21849932 (chromosome 20; LIME1; intron 3) and cg05353869 (chromosome 11; KLHL35; exon 2), we observed a -0.3 (95% CI: -0.4, -0.2) unit decrease (FDR = 0.04) and a 1.2 (95% CI: 0.7, 1.7) unit increase (FDR = 0.04), respectively, in percent DNAm per IQR (17.6 ppb) increase in 7-day mean ozone concentration. Results were not fully replicated in KORA and NAS. We identified three CpG sites potentially susceptible to gaseous air pollution-induced DNAm changes near genes relevant for cardiovascular and lung disease. Further harmonized investigations with a range of gaseous pollutants and averaging durations are needed to determine the effect of gaseous air pollutants on DNA methylation and ultimately gene expression.

Dynamic relationships between body fat and circulating adipokine levels from adolescence to young adulthood: The Santiago Longitudinal Study.

BACKGROUND AND AIMS: Adipose tissue secretes adipokines such as adiponectin and leptin, playing important roles in energy metabolism. The longitudinal associations between such adipokines and body fat accumulation have not been established, especially during adolescence and young adulthood and in diverse populations. The study aims to assess the longitudinal association between body fat measured with dual X-ray absorptiometry and plasma adipokines from adolescence to young adulthood. METHODS AND RESULTS: Among Hispanic/Latino participants (N = 537) aged 16.8 (SD: 0.3) years of the Santiago Longitudinal Study, we implemented structural equation modeling to estimate the sex-specific associations between adiposity (body fat percent (BF%) and proportion of trunk fat (PTF)) and adipokines (adiponectin and leptin levels) during adolescence (16 y) and these values after 6 years of follow-up (22 y). In addition, we further investigated whether the associations differed by baseline insulin resistance (IR) status. We found evidence for associations between 16 y BF% and 22 y leptin levels (ß (SE): 0.58 (0.06) for females; 0.53 (0.05) for males), between 16 y PTF and 22 y adiponectin levels (ß (SE): -0.31 (0.06) for females; -0.18 (0.06) for males) and between 16 y adiponectin levels and 22 y BF% (ß (SE): 0.12 (0.04) for both females and males). CONCLUSION: We observed dynamic relationships between adiposity and adipokines levels from late adolescence to young adulthood in a Hispanic/Latino population further demonstrating the importance of this period of the life course in the development of obesity.

Patient reported data integration for management of Eosinophilic Esophagitis.

Eosinophilic Esophagitis (EoE) is an esophageal allergic inflammatory disorder triggered by food proteins. Symptoms of EoE are variable within and between individuals. Presenting symptoms may include dysphagia, food bolus impaction, dyspepsia, or more subtle symptoms such as feeding disorders, regurgitation sensation, or nausea. The development and validation of a pediatric EoE patient self-reported and parent proxy-reported outcome symptom scoring tool was created by Franciosi et al. published in BMJ 2011, titled the Pediatric Eosinophilic Esophagitis Symptom Score (PEESS™ v2.0). To date, its use is largely for research purposes. We propose to evaluate the implementation of the PEESS™ v2.0 in a prospective interventional controlled clinical practice. The study included 620 patients over an 18-month period. Surveys were delivered and administered digitally every month through the MyGeisinger.org Patient Portal. Our analysis demonstrated symptom severity and symptom frequency scores significantly improved over time. However, counter to our hypothesis, patients who completed the PEESS™v2.0 ultimately had higher EoE-related health care utilization of office visits and endoscopies compared with those who did not complete the PEESS™v2.0. This could be related to greater awareness of disease activity and/or increased willingness to seek care. Our study, in the context of mobile health tool and patient-reported outcome trends, represents an opportunity for improved disease monitoring at-home within the field of eosinophilic gastrointestinal diseases.

New genetic loci link adipose and insulin biology to body fat distribution.

Body fat distribution is a heritable trait and a well-established predictor of adverse metabolic outcomes, independent of overall adiposity. To increase our understanding of the genetic basis of body fat distribution and its molecular links to cardiometabolic traits, here we conduct genome-wide association meta-analyses of traits related to waist and hip circumferences in up to 224,459 individuals. We identify 49 loci (33 new) associated with waist-to-hip ratio adjusted for body mass index (BMI), and an additional 19 loci newly associated with related waist and hip circumference measures (P < 5 × 10(-8)). In total, 20 of the 49 waist-to-hip ratio adjusted for BMI loci show significant sexual dimorphism, 19 of which display a stronger effect in women. The identified loci were enriched for genes expressed in adipose tissue and for putative regulatory elements in adipocytes. Pathway analyses implicated adipogenesis, angiogenesis, transcriptional regulation and insulin resistance as processes affecting fat distribution, providing insight into potential pathophysiological mechanisms.

Genetic studies of body mass index yield new insights for obesity biology.

Obesity is heritable and predisposes to many diseases. To understand the genetic basis of obesity better, here we conduct a genome-wide association study and Metabochip meta-analysis of body mass index (BMI), a measure commonly used to define obesity and assess adiposity, in up to 339,224 individuals. This analysis identifies 97 BMI-associated loci (P < 5 × 10(-8)), 56 of which are novel. Five loci demonstrate clear evidence of several independent association signals, and many loci have significant effects on other metabolic phenotypes. The 97 loci account for ∼2.7% of BMI variation, and genome-wide estimates suggest that common variation accounts for >20% of BMI variation. Pathway analyses provide strong support for a role of the central nervous system in obesity susceptibility and implicate new genes and pathways, including those related to synaptic function, glutamate signalling, insulin secretion/action, energy metabolism, lipid biology and adipogenesis.

Epigenetically mediated electrocardiographic manifestations of sub-chronic exposures to ambient particulate matter air pollution in the Women's Health Initiative and Atherosclerosis Risk in Communities Study.

BACKGROUND: Short-duration exposure to ambient particulate matter (PM) air pollution is associated with cardiac autonomic dysfunction and prolonged ventricular repolarization. However, associations with sub-chronic exposures to coarser particulates are relatively poorly characterized as are molecular mechanisms underlying their potential relationships with cardiovascular disease. MATERIALS AND METHODS: We estimated associations between monthly mean concentrations of PM < 10 µm and 2.5-10 µm in diameter (PM10; PM2.5-10) with time-domain measures of heart rate variability (HRV) and QT interval duration (QT) among U.S. women and men in the Women's Health Initiative and Atherosclerosis Risk in Communities Study (nHRV = 82,107; nQT = 76,711). Then we examined mediation of the PM-HRV and PM-QT associations by DNA methylation (DNAm) at three Cytosine-phosphate-Guanine (CpG) sites (cg19004594, cg24102420, cg12124767) with known sensitivity to monthly mean PM concentrations in a subset of the participants (nHRV = 7,169; nQT = 6,895). After multiply imputing missing PM, electrocardiographic and covariable data, we estimated associations using attrition-weighted, linear, mixed, longitudinal models adjusting for sociodemographic, behavioral, meteorological, and clinical characteristics. We assessed mediation by estimating the proportions of PM-HRV and PM-QT associations mediated by DNAm. RESULTS: We found little evidence of PM-HRV association, PM-QT association, or mediation by DNAm. CONCLUSIONS: The findings suggest that among racially/ethnically and environmentally diverse U.S. populations, sub-chronic exposures to coarser particulates may not exert appreciable, epigenetically mediated effects on cardiac autonomic function or ventricular repolarization. Further investigation in better-powered studies is warranted, with additional focus on shorter duration exposures to finer particulates and non-electrocardiographic outcomes among relatively susceptible populations.

Genome-wide meta-analysis of macronutrient intake of 91,114 European ancestry participants from the cohorts for heart and aging research in genomic epidemiology consortium.

Macronutrient intake, the proportion of calories consumed from carbohydrate, fat, and protein, is an important risk factor for metabolic diseases with significant familial aggregation. Previous studies have identified two genetic loci for macronutrient intake, but incomplete coverage of genetic variation and modest sample sizes have hindered the discovery of additional loci. Here, we expanded the genetic landscape of macronutrient intake, identifying 12 suggestively significant loci (P < 1 × 10-6) associated with intake of any macronutrient in 91,114 European ancestry participants. Four loci replicated and reached genome-wide significance in a combined meta-analysis including 123,659 European descent participants, unraveling two novel loci; a common variant in RARB locus for carbohydrate intake and a rare variant in DRAM1 locus for protein intake, and corroborating earlier FGF21 and FTO findings. In additional analysis of 144,770 participants from the UK Biobank, all identified associations from the two-stage analysis were confirmed except for DRAM1. Identified loci might have implications in brain and adipose tissue biology and have clinical impact in obesity-related phenotypes. Our findings provide new insight into biological functions related to macronutrient intake.

A genome-wide association study of polycystic ovary syndrome identified from electronic health records.

BACKGROUND: Polycystic ovary syndrome is the most common endocrine disorder affecting women of reproductive age. A number of criteria have been developed for clinical diagnosis of polycystic ovary syndrome, with the Rotterdam criteria being the most inclusive. Evidence suggests that polycystic ovary syndrome is significantly heritable, and previous studies have identified genetic variants associated with polycystic ovary syndrome diagnosed using different criteria. The widely adopted electronic health record system provides an opportunity to identify patients with polycystic ovary syndrome using the Rotterdam criteria for genetic studies. OBJECTIVE: To identify novel associated genetic variants under the same phenotype definition, we extracted polycystic ovary syndrome cases and unaffected controls based on the Rotterdam criteria from the electronic health records and performed a discovery-validation genome-wide association study. STUDY DESIGN: We developed a polycystic ovary syndrome phenotyping algorithm on the basis of the Rotterdam criteria and applied it to 3 electronic health record-linked biobanks to identify cases and controls for genetic study. In the discovery phase, we performed an individual genome-wide association study using the Geisinger MyCode and the Electronic Medical Records and Genomics cohorts, which were then meta-analyzed. We attempted validation of the significant association loci (P<1×10-6) in the BioVU cohort. All association analyses used logistic regression, assuming an additive genetic model, and adjusted for principal components to control for population stratification. An inverse-variance fixed-effect model was adopted for meta-analysis. In addition, we examined the top variants to evaluate their associations with each criterion in the phenotyping algorithm. We used the STRING database to characterize protein-protein interaction network. RESULTS: Using the same algorithm based on the Rotterdam criteria, we identified 2995 patients with polycystic ovary syndrome and 53,599 population controls in total (2742 cases and 51,438 controls from the discovery phase; 253 cases and 2161 controls in the validation phase). We identified 1 novel genome-wide significant variant rs17186366 (odds ratio [OR]=1.37 [1.23, 1.54], P=2.8×10-8) located near SOD2. In addition, 2 loci with suggestive association were also identified: rs113168128 (OR=1.72 [1.42, 2.10], P=5.2×10-8), an intronic variant of ERBB4 that is independent from the previously published variants, and rs144248326 (OR=2.13 [1.52, 2.86], P=8.45×10-7), a novel intronic variant in WWTR1. In the further association tests of the top 3 single-nucleotide polymorphisms with each criterion in the polycystic ovary syndrome algorithm, we found that rs17186366 (SOD2) was associated with polycystic ovaries and hyperandrogenism, whereas rs11316812 (ERBB4) and rs144248326 (WWTR1) were mainly associated with oligomenorrhea or infertility. We also validated the previously reported association with DENND1A1. Using the STRING database to characterize protein-protein interactions, we found both ERBB4 and WWTR1 can interact with YAP1, which has been previously associated with polycystic ovary syndrome. CONCLUSION: Through a discovery-validation genome-wide association study on polycystic ovary syndrome identified from electronic health records using an algorithm based on Rotterdam criteria, we identified and validated a novel genome-wide significant association with a variant near SOD2. We also identified a novel independent variant within ERBB4 and a suggestive association with WWTR1. With previously identified polycystic ovary syndrome gene YAP1, the ERBB4-YAP1-WWTR1 network suggests involvement of the epidermal growth factor receptor and the Hippo pathway in the multifactorial etiology of polycystic ovary syndrome.