Asthma Susceptibility Gene ORMDL3 Promotes Autophagy in Human Bronchial Epithelium

The genome-wide association study (GWAS)-identified asthma susceptibility risk alleles on chromosome 17q21 increase the expression of ORMDL3 (ORMDL sphingolipid biosynthesis regulator 3) in lung tissue. Given the importance of epithelial integrity in asthma, we hypothesized that ORMDL3 directly impacted bronchial epithelial function. To determine whether and how ORMDL3 expression impacts the bronchial epithelium, in studies using both primary human bronchial epithelial cells and human bronchial epithelial cell line, 16HBE (16HBE14o-), we assessed the impact of ORMDL3 on autophagy. Studies included: autophagosome detection by electron microscopy, RFP-GFP-LC3B to assess autophagic activity, and Western blot analysis of autophagy-related proteins. Mechanistic assessments included immunoprecipitation assays, intracellular calcium mobilization assessments, and cell viability assays. Coexpression of ORMDL3 and autophagy-related genes was measured in primary human bronchial epithelial cells derived from 44 subjects. Overexpressing ORMDL3 demonstrated increased numbers of autophagosomes and increased levels of autophagy-related proteins LC3B, ATG3, ATG7, and ATG16L1. ORMDL3 overexpression promotes autophagy and subsequent cell death by impairing intracellular calcium mobilization through interacting with SERCA2. Strong correlation was observed between expression of ORMDL3 and autophagy-related genes in patient-derived bronchial epithelial cells. Increased ORMDL3 expression induces autophagy, possibly through interacting with SERCA2, thereby inhibiting intracellular calcium influx, and induces cell death, impairing bronchial epithelial function in asthma.

Early Changes in Immune Cell Count, Metabolism, and Function Following Sleeve Gastrectomy: A Prospective Human Study

OBJECTIVE

To characterize longitudinal changes in blood biomarkers, leukocyte composition, and gene expression following laparoscopic sleeve gastrectomy (LSG).

BACKGROUND

LSG is an effective treatment for obesity, leading to sustainable weight loss and improvements in obesity-related comorbidities and inflammatory profiles. However, the effects of LSG on immune function and metabolism remain uncertain.

METHODS

Prospective data were collected from 23 enrolled human subjects from a single institution. Parameters of weight, comorbidities, and trends in blood biomarkers and leukocyte subsets were observed from preoperative baseline to 1 year postsurgery in 3-month follow-up intervals. RNA sequencing was performed on pairs of whole blood samples from the first 6 subjects of the study (baseline and 3 months postsurgery) to identify genome-wide gene expression changes associated with undergoing LSG.

RESULTS

LSG led to a significant decrease in mean total body weight loss (18.1%) at 3 months and among diabetic subjects a reduction in hemoglobin A1c. Improvements in clinical inflammatory and hormonal biomarkers were demonstrated as early as 3 months after LSG. A reduction in neutrophil-lymphocyte ratio was observed, driven by a reduction in absolute neutrophil counts. Gene set enrichment analyses of differential whole blood gene expression demonstrated that after 3 months LSG induced transcriptomic changes not only in inflammatory cytokine pathways but also in several key metabolic pathways related to energy metabolism.

CONCLUSIONS

LSG induces significant changes in the composition and metabolism of immune cells as early as 3 months postoperatively. Further evaluation is required of bariatric surgery's effects on immunometabolism and the consequences for host defense and metabolic disease.

Early changes in immune cell metabolism and function are a hallmark of sleeve gastrectomy: a prospective human study

OBJECTIVE

To characterize longitudinal changes in blood biomarkers, leukocyte composition, and gene expression following laparoscopic sleeve gastrectomy (LSG).

BACKGROUND

LSG is an effective treatment for obesity, leading to sustainable weight loss and improvements in obesity-related co-morbidities and inflammatory profiles. However, the effects of LSG on immune function and metabolism remain uncertain.

METHODS

Prospective data was collected from 23 enrolled human subjects from a single institution. Parameters of weight, co-morbidities, and trends in blood biomarkers and leukocyte subsets were observed from pre-operative baseline to one year in three-month follow-up intervals. RNA-sequencing was performed on pairs of whole blood samples from the first six subjects of the study (baseline and three months post-surgery) to identify genome-wide gene expression changes associated with undergoing LSG.

RESULTS

LSG led to a significant decrease in mean total body weight loss (18.1%) at three months and among diabetic subjects a reduction in HbA1c. Improvements in clinical inflammatory and hormonal biomarkers were demonstrated as early as three months after LSG. A reduction in neutrophil-lymphocyte ratio was observed, driven by a reduction in absolute neutrophil counts. Gene set enrichment analyses of differential whole blood gene expression demonstrated that after three months, LSG induced transcriptomic changes not only in inflammatory cytokine pathways but also in several key metabolic pathways related to energy metabolism.

CONCLUSIONS

LSG induces significant changes in the composition and metabolism of immune cells as early as three months post-operatively. Further evaluation is required of bariatric surgery's effects on immunometabolism and consequences for host defense and metabolic disease.

Expression of SMARCD1 interacts with age in association with asthma control on inhaled corticosteroid therapy

Background

Global gene expression levels are known to be highly dependent upon gross demographic features including age, yet identification of age-related genomic indicators has yet to be comprehensively undertaken in a disease and treatment-specific context.

Methods

We used gene expression data from CD4+ lymphocytes in the Asthma BioRepository for Integrative Genomic Exploration (Asthma BRIDGE), an open-access collection of subjects participating in genetic studies of asthma with available gene expression data. Replication population participants were Puerto Rico islanders recruited as part of the ongoing Genes environments & Admixture in Latino Americans (GALA II), who provided nasal brushings for transcript sequencing. The main outcome measure was chronic asthma control as derived by questionnaires. Genomic associations were performed using regression of chronic asthma control score on gene expression with age in years as a covariate, including a multiplicative interaction term for gene expression times age.

Results

The SMARCD1 gene (SWI/SNF-related matrix-associated actin-dependent regulator of chromatin subfamily D member 1) interacted with age to influence chronic asthma control on inhaled corticosteroids, with a doubling of expression leading to an increase of 1.3 units of chronic asthma control per year (95% CI [0.86, 1.74], p = 6 × 10^− 9), suggesting worsening asthma control with increasing age. This result replicated in GALA II (p = 3.8 × 10^− 8). Cellular assays confirmed the role of SMARCD1 in glucocorticoid response in airway epithelial cells.

Conclusion

Focusing on age-dependent factors may help identify novel indicators of asthma medication response. Age appears to modulate the effect of SMARCD1 on asthma control with inhaled corticosteroids.

Pharmacogenomics and Placebo Response in a Randomized Clinical Trial in Asthma

Genetic variation may differentially modify drug and placebo treatment effects in randomized clinical trials. In asthma, although lung function and asthma control improvements are commonplace with placebo, pharmacogenomics of placebo vs. drug response remains unexamined. In a genomewide association study of subjective and objective outcomes with placebo treatment in Childhood Asthma Management Program of nedocromil/budesonide vs. placebo (N = 604), effect estimates for lead single nucleotide polymorphisms (SNPs) were compared across arms. The coughing/wheezing lead SNP, rs2392165 (β = 0.94; P = 1.10E-07) mapped to BBS9, a gene implicated in lung development that contains a lung function expression quantitative trait locus. The effect was attenuated with budesonide (P_interaction  = 1.48E-07), but not nedocromil (P_interaction  = 0.06). The lead forced vital capacity SNP, rs12930749 (β = -5.80; P = 1.47E-06), mapped to KIAA0556, a locus genomewide associated with respiratory diseases. The rs12930749 effect was attenuated with budesonide (P_interaction  = 1.32E-02) and nedocromil (P_interaction  = 1.09E-02). Pharmacogenomic analysis revealed differential effects with placebo and drug treatment that could potentially guide precision drug development in asthma.

Pathway analysis of a genome-wide gene by air pollution interaction study in asthmatic children

OBJECTIVES

We aimed to investigate the role of genetics in the respiratory response of asthmatic children to air pollution, with a genome-wide level analysis of gene by nitrogen dioxide (NO2) and carbon monoxide (CO) interaction on lung function and to identify biological pathways involved.

METHODS

We used a two-step method for fast linear mixed model computations for genome-wide association studies, exploring whether variants modify the longitudinal relationship between 4-month average pollution and post-bronchodilator FEV1 in 522 Caucasian and 88 African-American asthmatic children. Top hits were confirmed with classic linear mixed-effect models. We used the improved gene set enrichment analysis for GWAS (i-GSEA4GWAS) to identify plausible pathways.

RESULTS

Two SNPs near the EPHA3 (rs13090972 and rs958144) and one in TXNDC8 (rs7041938) showed significant interactions with NO2 in Caucasians but we did not replicate this locus in African-Americans. SNP-CO interactions did not reach genome-wide significance. The i-GSEA4GWAS showed a pathway linked to the HO-1/CO system to be associated with CO-related FEV1 changes. For NO2-related FEV1 responses, we identified pathways involved in cellular adhesion, oxidative stress, inflammation, and metabolic responses.

CONCLUSION

The host lung function response to long-term exposure to pollution is linked to genes involved in cellular adhesion, oxidative stress, inflammatory, and metabolic pathways.

Meta-analysis of genome-wide association studies for body fat distribution in 694 649 individuals of European ancestry

More than one in three adults worldwide is either overweight or obese. Epidemiological studies indicate that the location and distribution of excess fat, rather than general adiposity, are more informative for predicting risk of obesity sequelae, including cardiometabolic disease and cancer. We performed a genome-wide association study meta-analysis of body fat distribution, measured by waist-to-hip ratio (WHR) adjusted for body mass index (WHRadjBMI), and identified 463 signals in 346 loci. Heritability and variant effects were generally stronger in women than men, and we found approximately one-third of all signals to be sexually dimorphic. The 5% of individuals carrying the most WHRadjBMI-increasing alleles were 1.62 times more likely than the bottom 5% to have a WHR above the thresholds used for metabolic syndrome. These data, made publicly available, will inform the biology of body fat distribution and its relationship with disease.

Gene Coexpression Networks in Whole Blood Implicate Multiple Interrelated Molecular Pathways in Obesity in People with Asthma

OBJECTIVE

Asthmatic children who develop obesity through adolescence have poorer disease outcomes compared with those who do not. This study aimed to characterize the biology of childhood asthma complicated by adult obesity.

METHODS

Gene expression networks are powerful statistical tools for characterizing human disease that leverage the putative coregulatory relationships of genes to infer relevant biological pathways. Weighted gene coexpression network analysis of gene expression data was performed in whole blood from 514 adult asthmatic subjects. Then, module preservation and association replication analyses were performed in 418 subjects from two independent asthma cohorts (one pediatric and one adult).

RESULTS

A multivariate model was identified in which three gene coexpression network modules were associated with incident obesity in the discovery cohort (each P < 0.05). Two module memberships were enriched for genes in pathways related to platelets, integrins, extracellular matrix, smooth muscle, NF-κB signaling, and Hedgehog signaling. The network structures of each of the obesity modules were significantly preserved in both replication cohorts (permutation P = 9.999E-05). The corresponding module gene sets were significantly enriched for differential expression in subjects with obesity in both replication cohorts (each P < 0.05).

CONCLUSIONS

The gene coexpression network profiles thus implicate multiple interrelated pathways in the biology of an important endotype of asthma with obesity.

Novel eosinophilic gene expression networks associated with IgE in two distinct asthma populations

BACKGROUND

Asthma represents a significant public health burden; however, novel biological therapies targeting immunoglobulin E (IgE)-mediated pathways have widened clinical treatment options for the disease.

OBJECTIVE

In this study, we sought to identify gene transcripts and gene networks involved in the determination of serum IgE levels in people with asthma that can help inform the development of novel therapeutic agents.

METHODS

We analysed gene expression data from a cross-sectional study of 326 Costa Rican children with asthma, aged 6 to 12 years, from the Genetics of Asthma in Costa Rica Study and 610 young adults with asthma, aged 16 to 25 years, from the Childhood Asthma Management Program trial. We utilized differential gene expression analysis and performed weighted gene coexpression network analysis on 25 060 genes, to identify gene transcripts and network modules associated with total IgE, adjusting for age and gender. We used pathway enrichment analyses to identify key biological pathways underlying significant modules. We compared findings that replicated between both populations.

RESULTS

We identified 31 transcripts associated with total IgE that replicated between the two study cohorts. These results were notable for increased eosinophil-related transcripts (including IL5RA, CLC, SMPD3, CCL23 and CEBPE). Pathway enrichment identified the regulation of T cell tolerance as important in the determination of total IgE levels, supporting a key role for IDO1.

CONCLUSIONS AND CLINICAL RELEVANCE

These results provide robust evidence that biologically meaningful gene expression profiles (relating to eosinophilic and regulatory T cell pathways in particular) associated with total IgE levels can be identified in individuals diagnosed with asthma during childhood. These profiles and their constituent genes may represent novel therapeutic targets.

The eMERGE genotype set of 83,717 subjects imputed to ~40 million variants genome wide and association with the herpes zoster medical record phenotype

The Electronic Medical Records and Genomics (eMERGE) network is a network of medical centers with electronic medical records linked to existing biorepository samples for genomic discovery and genomic medicine research. The network sought to unify the genetic results from 78 Illumina and Affymetrix genotype array batches from 12 contributing medical centers for joint association analysis of 83,717 human participants. In this report, we describe the imputation of eMERGE results and methods to create the unified imputed merged set of genome‐wide variant genotype data. We imputed the data using the Michigan Imputation Server, which provides a missing single‐nucleotide variant genotype imputation service using the minimac3 imputation algorithm with the Haplotype Reference Consortium genotype reference set. We describe the quality control and filtering steps used in the generation of this data set and suggest generalizable quality thresholds for imputation and phenotype association studies. To test the merged imputed genotype set, we replicated a previously reported chromosome 6 HLA‐B herpes zoster (shingles) association and discovered a novel zoster‐associated loci in an epigenetic binding site near the terminus of chromosome 3 (3p29).

An Integrative Transcriptomic and Metabolomic Study of Lung Function in Children With Asthma

BACKGROUND

Single omic analyses have provided some insight into the basis of lung function in children with asthma, but the underlying biologic pathways are still poorly understood.

METHODS

Weighted gene coexpression network analysis (WGCNA) was used to identify modules of coregulated gene transcripts and metabolites in blood among 325 children with asthma from the Genetic Epidemiology of Asthma in Costa Rica study. The biology of modules associated with lung function as measured by FEV₁, the FEV₁/FVC ratio, bronchodilator response, and airway responsiveness to methacholine was explored. Significantly correlated gene-metabolite module pairs were then identified, and their constituent features were analyzed for biologic pathway enrichments.

RESULTS

WGCNA clustered 25,060 gene probes and 8,185 metabolite features into eight gene modules and eight metabolite modules, where four and six, respectively, were associated with lung function (P ≤ .05). The gene modules were enriched for immune, mitotic, and metabolic processes and asthma-associated microRNA targets. The metabolite modules were enriched for lipid and amino acid metabolism. Integration of correlated gene-metabolite modules expanded the single omic findings, linking the FEV₁/FVC ratio with ORMDL3 and dysregulated lipid metabolism. This finding was replicated in an independent population.

CONCLUSIONS

The results of this hypothesis-generating study suggest a mechanistic basis for multiple asthma genes, including ORMDL3, and a role for lipid metabolism. They demonstrate that integrating multiple omic technologies may provide a more informative picture of asthmatic lung function biology than single omic analyses.

Genetics and Genomics of Longitudinal Lung Function Patterns in Individuals with Asthma

RATIONALE

Patterns of longitudinal lung function growth and decline in childhood asthma have been shown to be important in determining risk for future respiratory ailments including chronic airway obstruction and chronic obstructive pulmonary disease.

OBJECTIVES

To determine the genetic underpinnings of lung function patterns in subjects with childhood asthma.

METHODS

We performed a genome-wide association study of 581 non-Hispanic white individuals with asthma that were previously classified by patterns of lung function growth and decline (normal growth, normal growth with early decline, reduced growth, and reduced growth with early decline). The strongest association was also measured in two additional cohorts: a small asthma cohort and a large chronic obstructive pulmonary disease metaanalysis cohort. Interaction between the genomic region encompassing the most strongly associated single-nucleotide polymorphism and nearby genes was assessed by two chromosome conformation capture assays.

MEASUREMENTS AND MAIN RESULTS

An intergenic single-nucleotide polymorphism (rs4445257) on chromosome 8 was strongly associated with the normal growth with early decline pattern compared with all other pattern groups (P = 6.7 × 10^-9; odds ratio, 2.8; 95% confidence interval, 2.0-4.0); replication analysis suggested this variant had opposite effects in normal growth with early decline and reduced growth with early decline pattern groups. Chromosome conformation capture experiments indicated a chromatin interaction between rs4445257 and the promoter of the distal CSMD3 gene.

CONCLUSIONS

Early decline in lung function after normal growth is associated with a genetic polymorphism that may also protect against early decline in reduced growth groups. Clinical trial registered with www.clinicaltrials.gov (NCT00000575).

Gene Expression Profiling in Blood Provides Reproducible Molecular Insights into Asthma Control

RATIONALE

Maintaining optimal symptom control remains the primary objective of asthma treatment. Better understanding of the biologic underpinnings of asthma control may lead to the development of improved clinical and pharmaceutical approaches.

OBJECTIVES

To identify molecular pathways and interrelated genes whose differential expression was associated with asthma control.

METHODS

We performed gene set enrichment analyses of asthma control in 1,170 adults with asthma, each with gene expression data derived from either whole blood (WB) or unstimulated CD4+ T lymphocytes (CD4), and a self-reported asthma control score representing either the preceding 6 months (chronic) or 7 days (acute). Our study comprised a discovery WB cohort (n = 245, chronic) and three independent, nonoverlapping replication cohorts: a second WB set (n = 448, acute) and two CD4 sets (n = 300, chronic; n = 77, acute).

MEASUREMENTS AND MAIN RESULTS

In the WB discovery cohort, we found significant overrepresentation of genes associated with asthma control in 1,106 gene sets from the Molecular Signatures Database (false discovery rate, <5%). Of these, 583 (53%) replicated in at least one replication cohort (false discovery rate, <25%). Suboptimal control was associated with signatures of eosinophilic and granulocytic inflammatory signals, whereas optimal control signatures were enriched for immature lymphocytic patterns. These signatures included two related biologic processes related to activation by TREM-1 (triggering receptor expressed on myeloid cells 1) and lipopolysaccharide.

CONCLUSIONS

Together, these results demonstrate the existence of specific, reproducible transcriptomic components in blood that vary with degree of asthma control and implicate a novel biologic target (TREM-1).

Integration of metabolomic and transcriptomic networks in pregnant women reveals biological pathways and predictive signatures associated with preeclampsia

INTRODUCTION

Preeclampsia is a leading cause of maternal and fetal mortality worldwide, yet its exact pathogenesis remains elusive.

OBJECTIVES

This study, nested within the Vitamin D Antenatal Asthma Reduction Trial (VDAART), aimed to develop integrated omics models of preeclampsia that have utility in both prediction and in the elucidation of underlying biological mechanisms.

METHODS

Metabolomic profiling was performed on first trimester plasma samples of 47 pregnant women from VDAART who subsequently developed preeclampsia and 62 controls with healthy pregnancies, using liquid-chromatography tandem mass-spectrometry. Metabolomic profiles were generated based on logistic regression models and assessed using Received Operator Characteristic Curve analysis. These profiles were compared to profiles from generated using third trimester samples. The first trimester metabolite profile was then integrated with a pre-existing transcriptomic profile using network methods.

RESULTS

In total, 72 (0.9%) metabolite features were associated (p<0.01) with preeclampsia after adjustment for maternal age, race, and gestational age. These features had moderate to good discriminatory ability; in ROC curve analyses a summary score based on these features displayed an area under the curve (AUC) of 0.794 (95%CI 0.700, 0.888). This profile retained the ability to distinguish preeclamptic from healthy pregnancies in the third trimester (AUC:0.762 (95% CI 0.663, 0.860)). Additionally, metabolite set enrichment analysis identified common pathways, including glycerophospholipid metabolism, at the two time-points. Integration with the transcriptomic signature refined these results suggesting a particular role for lipid imbalance, immune function and the circulatory system.

CONCLUSIONS

These findings suggest it is possible to develop a predictive metabolomic profile of preeclampsia. This profile is characterized by changes in lipid and amino acid metabolism and dysregulation of immune response and can be refined through interaction with transcriptomic data. However validation in larger and more diverse populations is required.

Whole genome prediction and heritability of childhood asthma phenotypes

Introduction

While whole genome prediction (WGP) methods have recently demonstrated successes in the prediction of complex genetic diseases, they have not yet been applied to asthma and related phenotypes. Longitudinal patterns of lung function differ between asthmatics, but these phenotypes have not been assessed for heritability or predictive ability. Herein, we assess the heritability and genetic predictability of asthma‐related phenotypes.

Methods

We applied several WGP methods to a well‐phenotyped cohort of 832 children with mild‐to‐moderate asthma from CAMP. We assessed narrow‐sense heritability and predictability for airway hyperresponsiveness, serum immunoglobulin E, blood eosinophil count, pre‐ and post‐bronchodilator forced expiratory volume in 1 sec (FEV₁), bronchodilator response, steroid responsiveness, and longitudinal patterns of lung function (normal growth, reduced growth, early decline, and their combinations). Prediction accuracy was evaluated using a training/testing set split of the cohort.

Results

We found that longitudinal lung function phenotypes demonstrated significant narrow‐sense heritability (reduced growth, 95%; normal growth with early decline, 55%). These same phenotypes also showed significant polygenic prediction (areas under the curve [AUCs] 56% to 62%). Including additional demographic covariates in the models increased prediction 4–8%, with reduced growth increasing from 62% to 66% AUC. We found that prediction with a genomic relatedness matrix was improved by filtering available SNPs based on chromatin evidence, and this result extended across cohorts.

Conclusions

Longitudinal reduced lung function growth displayed extremely high heritability. All phenotypes with significant heritability showed significant polygenic prediction. Using SNP‐prioritization increased prediction across cohorts. WGP methods show promise in predicting asthma‐related heritable traits.

Across-cohort QC analyses of GWAS summary statistics from complex traits

Genome-wide association studies (GWASs) have been successful in discovering SNP trait associations for many quantitative traits and common diseases. Typically, the effect sizes of SNP alleles are very small and this requires large genome-wide association meta-analyses (GWAMAs) to maximize statistical power. A trend towards ever-larger GWAMA is likely to continue, yet dealing with summary statistics from hundreds of cohorts increases logistical and quality control problems, including unknown sample overlap, and these can lead to both false positive and false negative findings. In this study, we propose four metrics and visualization tools for GWAMA, using summary statistics from cohort-level GWASs. We propose methods to examine the concordance between demographic information, and summary statistics and methods to investigate sample overlap. (I) We use the population genetics F_st statistic to verify the genetic origin of each cohort and their geographic location, and demonstrate using GWAMA data from the GIANT Consortium that geographic locations of cohorts can be recovered and outlier cohorts can be detected. (II) We conduct principal component analysis based on reported allele frequencies, and are able to recover the ancestral information for each cohort. (III) We propose a new statistic that uses the reported allelic effect sizes and their standard errors to identify significant sample overlap or heterogeneity between pairs of cohorts. (IV) To quantify unknown sample overlap across all pairs of cohorts, we propose a method that uses randomly generated genetic predictors that does not require the sharing of individual-level genotype data and does not breach individual privacy.

Patterns of Growth and Decline in Lung Function in Persistent Childhood Asthma

BACKGROUND

Tracking longitudinal measurements of growth and decline in lung function in patients with persistent childhood asthma may reveal links between asthma and subsequent chronic airflow obstruction.

METHODS

We classified children with asthma according to four characteristic patterns of lung-function growth and decline on the basis of graphs showing forced expiratory volume in 1 second (FEV1), representing spirometric measurements performed from childhood into adulthood. Risk factors associated with abnormal patterns were also examined. To define normal values, we used FEV1 values from participants in the National Health and Nutrition Examination Survey who did not have asthma.

RESULTS

Of the 684 study participants, 170 (25%) had a normal pattern of lung-function growth without early decline, and 514 (75%) had abnormal patterns: 176 (26%) had reduced growth and an early decline, 160 (23%) had reduced growth only, and 178 (26%) had normal growth and an early decline. Lower baseline values for FEV1, smaller bronchodilator response, airway hyperresponsiveness at baseline, and male sex were associated with reduced growth (P<0.001 for all comparisons). At the last spirometric measurement (mean [±SD] age, 26.0±1.8 years), 73 participants (11%) met Global Initiative for Chronic Obstructive Lung Disease spirometric criteria for lung-function impairment that was consistent with chronic obstructive pulmonary disease (COPD); these participants were more likely to have a reduced pattern of growth than a normal pattern (18% vs. 3%, P<0.001).

CONCLUSIONS

Childhood impairment of lung function and male sex were the most significant predictors of abnormal longitudinal patterns of lung-function growth and decline. Children with persistent asthma and reduced growth of lung function are at increased risk for fixed airflow obstruction and possibly COPD in early adulthood. (Funded by the Parker B. Francis Foundation and others; ClinicalTrials.gov number, NCT00000575.).

Expression Quantitative Trait Loci Information Improves Predictive Modeling of Disease Relevance of Non-Coding Genetic Variation

Disease-associated loci identified through genome-wide association studies (GWAS) frequently localize to non-coding sequence. We and others have demonstrated strong enrichment of such single nucleotide polymorphisms (SNPs) for expression quantitative trait loci (eQTLs), supporting an important role for regulatory genetic variation in complex disease pathogenesis. Herein we describe our initial efforts to develop a predictive model of disease-associated variants leveraging eQTL information. We first catalogued cis-acting eQTLs (SNPs within 100kb of target gene transcripts) by meta-analyzing four studies of three blood-derived tissues (n = 586). At a false discovery rate < 5%, we mapped eQTLs for 6,535 genes; these were enriched for disease-associated genes (P < 10⁻⁰⁴), particularly those related to immune diseases and metabolic traits. Based on eQTL information and other variant annotations (distance from target gene transcript, minor allele frequency, and chromatin state), we created multivariate logistic regression models to predict SNP membership in reported GWAS. The complete model revealed independent contributions of specific annotations as strong predictors, including evidence for an eQTL (odds ratio (OR) = 1.2–2.0, P < 10⁻¹¹) and the chromatin states of active promoters, different classes of strong or weak enhancers, or transcriptionally active regions (OR = 1.5–2.3, P < 10⁻¹¹). This complete prediction model including eQTL association information ultimately allowed for better discrimination of SNPs with higher probabilities of GWAS membership (6.3–10.0%, compared to 3.5% for a random SNP) than the other two models excluding eQTL information. This eQTL-based prediction model of disease relevance can help systematically prioritize non-coding GWAS SNPs for further functional characterization.

Stress and Bronchodilator Response in Children with Asthma

RATIONALE

Stress is associated with asthma morbidity in Puerto Ricans (PRs), who have reduced bronchodilator response (BDR).

OBJECTIVES

To examine whether stress and/or a gene regulating anxiety (ADCYAP1R1) is associated with BDR in PR and non-PR children with asthma.

METHODS

This was a cross-sectional study of stress and BDR (percent change in FEV1 after BD) in 234 PRs ages 9-14 years with asthma. We assessed child stress using the Checklist of Children's Distress Symptoms, and maternal stress using the Perceived Stress Scale. Replication analyses were conducted in two cohorts. Polymorphisms in ADCYAP1R1 were genotyped in our study and six replication studies. Multivariable models of stress and BDR were adjusted for age, sex, income, environmental tobacco smoke, and use of inhaled corticosteroids.

MEASUREMENTS AND MAIN RESULTS

High child stress was associated with reduced BDR in three cohorts. PR children who were highly stressed (upper quartile, Checklist of Children's Distress Symptoms) and whose mothers had high stress (upper quartile, Perceived Stress Scale) had a BDR that was 10.2% (95% confidence interval, 6.1-14.2%) lower than children who had neither high stress nor a highly stressed mother. A polymorphism in ADCYAP1R1 (rs34548976) was associated with reduced BDR. This single-nucleotide polymorphism is associated with reduced expression of the gene for the β2-adrenergic receptor (ADRB2) in CD4(+) lymphocytes of subjects with asthma, and it affects brain connectivity of the amygdala and the insula (a biomarker of anxiety).

CONCLUSIONS

High child stress and an ADCYAP1R1 single-nucleotide polymorphism are associated with reduced BDR in children with asthma. This is likely caused by down-regulation of ADRB2 in highly stressed children.

The metabolomics of asthma control: a promising link between genetics and disease

Short-acting β agonists (e.g., albuterol) are the most commonly used medications for asthma, a disease that affects over 300 million people in the world. Metabolomic profiling of asthmatics taking β agonists presents a new and promising resource for identifying the molecular determinants of asthma control. The objective is to identify novel genetic and biochemical predictors of asthma control using an integrative "omics" approach. We generated lipidomic data by liquid chromatography tandem mass spectrometry (LC-MS), - using plasma samples from 20 individuals with asthma. The outcome of interest was a binary indicator of asthma control defined by the use of albuterol inhalers in the preceding week. We integrated metabolomic data with genome-wide genotype, gene expression, and methylation data of this cohort to identify genomic and molecular indicators of asthma control. A Conditional Gaussian Bayesian Network (CGBN) was generated using the strongest predictors from each of these analyses. Integrative and metabolic pathway over-representation analyses (ORA) identified enrichment of known biological pathways within the strongest molecular determinants. Of the 64 metabolites measured, 32 had known identities. The CGBN model based on four SNPs (rs9522789, rs7147228, rs2701423, rs759582) and two metabolites-monoHETE_0863 and sphingosine-1-phosphate (S1P) could predict asthma control with an AUC of 95%. Integrative ORA identified 17 significantly enriched pathways related to cellular immune response, interferon signaling, and cytokine-related signaling, for which arachidonic acid, PGE2 and S1P, in addition to six genes (CHN1, PRKCE, GNA12, OASL, OAS1, and IFIT3) appeared to drive the pathway results. Of these predictors, S1P, GNA12, and PRKCE were enriched in the results from integrative and metabolic ORAs. Through an integrative analysis of metabolomic, genomic, and methylation data from a small cohort of asthmatics, we implicate altered metabolic pathways, related to sphingolipid metabolism, in asthma control. These results provide insight into the pathophysiology of asthma control.

Genome-wide expression profiles identify potential targets for gene-environment interactions in asthma severity

BACKGROUND

Gene-environment interaction studies using genome-wide association study data are often underpowered after adjustment for multiple comparisons. Differential gene expression in response to the exposure of interest can capture the most biologically relevant genes at the genome-wide level.

OBJECTIVE

We used differential genome-wide expression profiles from the Epidemiology of Home Allergens and Asthma birth cohort in response to Der f 1 allergen (sensitized vs nonsensitized) to inform a gene-environment study of dust mite exposure and asthma severity.

METHODS

Polymorphisms in differentially expressed genes were identified in genome-wide association study data from the Childhood Asthma Management Program, a clinical trial in childhood asthmatic patients. Home dust mite allergen levels (<10 or ≥10 μg/g dust) were assessed at baseline, and (≥1) severe asthma exacerbation (emergency department visit or hospitalization for asthma in the first trial year) served as the disease severity outcome. The Genetics of Asthma in Costa Rica Study and a Puerto Rico/Connecticut asthma cohort were used for replication.

RESULTS

IL9, IL5, and proteoglycan 2 expression (PRG2) was upregulated in Der f 1-stimulated PBMCs from dust mite-sensitized patients (adjusted P < .04). IL9 polymorphisms (rs11741137, rs2069885, and rs1859430) showed evidence for interaction with dust mite in the Childhood Asthma Management Program (P = .02 to .03), with replication in the Genetics of Asthma in Costa Rica Study (P = .04). Subjects with the dominant genotype for these IL9 polymorphisms were more likely to report a severe asthma exacerbation if exposed to increased dust mite levels.

CONCLUSIONS

Genome-wide differential gene expression in response to dust mite allergen identified IL9, a biologically plausible gene target that might interact with environmental dust mite to increase severe asthma exacerbations in children.

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, we conducted genome-wide association meta-analyses of waist and hip circumference-related traits in up to 224,459 individuals. We identified 49 loci (33 new) associated with waist-to-hip ratio adjusted for body mass index (WHRadjBMI) and an additional 19 loci newly associated with related waist and hip circumference measures (P<5×10⁻⁸). Twenty of the 49 WHRadjBMI loci showed significant sexual dimorphism, 19 of which displayed 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.

Defining the role of common variation in the genomic and biological architecture of adult human height

Using genome-wide data from 253,288 individuals, we identified 697 variants at genome-wide significance that together explained one-fifth of the heritability for adult height. By testing different numbers of variants in independent studies, we show that the most strongly associated ∼2,000, ∼3,700 and ∼9,500 SNPs explained ∼21%, ∼24% and ∼29% of phenotypic variance. Furthermore, all common variants together captured 60% of heritability. The 697 variants clustered in 423 loci were enriched for genes, pathways and tissue types known to be involved in growth and together implicated genes and pathways not highlighted in earlier efforts, such as signaling by fibroblast growth factors, WNT/β-catenin and chondroitin sulfate-related genes. We identified several genes and pathways not previously connected with human skeletal growth, including mTOR, osteoglycin and binding of hyaluronic acid. Our results indicate a genetic architecture for human height that is characterized by a very large but finite number (thousands) of causal variants.

Genome-wide interaction studies reveal sex-specific asthma risk alleles

Asthma is a complex disease with sex-specific differences in prevalence. Candidate gene studies have suggested that genotype-by-sex interaction effects on asthma risk exist, but this has not yet been explored at a genome-wide level. We aimed to identify sex-specific asthma risk alleles by performing a genome-wide scan for genotype-by-sex interactions in the ethnically diverse participants in the EVE Asthma Genetics Consortium. We performed male- and female-specific genome-wide association studies in 2653 male asthma cases, 2566 female asthma cases and 3830 non-asthma controls from European American, African American, African Caribbean and Latino populations. Association tests were conducted in each study sample, and the results were combined in ancestry-specific and cross-ancestry meta-analyses. Six sex-specific asthma risk loci had P-values < 1 × 10(-6), of which two were male specific and four were female specific; all were ancestry specific. The most significant sex-specific association in European Americans was at the interferon regulatory factor 1 (IRF1) locus on 5q31.1. We also identify a Latino female-specific association in RAP1GAP2. Both of these loci included single-nucleotide polymorphisms that are known expression quantitative trait loci and have been associated with asthma in independent studies. The IRF1 locus is a strong candidate region for male-specific asthma susceptibility due to the association and validation we demonstrate here, the known role of IRF1 in asthma-relevant immune pathways and prior reports of sex-specific differences in interferon responses.

Quality control and conduct of genome-wide association meta-analyses

Rigorous organization and quality control (QC) are necessary to facilitate successful genome-wide association meta-analyses (GWAMAs) of statistics aggregated across multiple genome-wide association studies. This protocol provides guidelines for (i) organizational aspects of GWAMAs, and for (ii) QC at the study file level, the meta-level across studies and the meta-analysis output level. Real-world examples highlight issues experienced and solutions developed by the GIANT Consortium that has conducted meta-analyses including data from 125 studies comprising more than 330,000 individuals. We provide a general protocol for conducting GWAMAs and carrying out QC to minimize errors and to guarantee maximum use of the data. We also include details for the use of a powerful and flexible software package called EasyQC. Precise timings will be greatly influenced by consortium size. For consortia of comparable size to the GIANT Consortium, this protocol takes a minimum of about 10 months to complete.

A meta-analysis of genome-wide association studies for adiponectin levels in East Asians identifies a novel locus near WDR11-FGFR2

Blood levels of adiponectin, an adipocyte-secreted protein correlated with metabolic and cardiovascular risks, are highly heritable. Genome-wide association (GWA) studies for adiponectin levels have identified 14 loci harboring variants associated with blood levels of adiponectin. To identify novel adiponectin-associated loci, particularly those of importance in East Asians, we conducted a meta-analysis of GWA studies for adiponectin in 7827 individuals, followed by two stages of replications in 4298 and 5954 additional individuals. We identified a novel adiponectin-associated locus on chromosome 10 near WDR11-FGFR2 (P = 3.0 × 10(-14)) and provided suggestive evidence for a locus on chromosome 12 near OR8S1-LALBA (P = 1.2 × 10(-7)). Of the adiponectin-associated loci previously described, we confirmed the association at CDH13 (P = 6.8 × 10(-165)), ADIPOQ (P = 1.8 × 10(-22)), PEPD (P = 3.6 × 10(-12)), CMIP (P = 2.1 × 10(-10)), ZNF664 (P = 2.3 × 10(-7)) and GPR109A (P = 7.4 × 10(-6)). Conditional analysis at ADIPOQ revealed a second signal with suggestive evidence of association only after conditioning on the lead SNP (Pinitial = 0.020; Pconditional = 7.0 × 10(-7)). We further confirmed the independence of two pairs of closely located loci (<2 Mb) on chromosome 16 at CMIP and CDH13, and on chromosome 12 at GPR109A and ZNF664. In addition, the newly identified signal near WDR11-FGFR2 exhibited evidence of association with triglycerides (P = 3.3 × 10(-4)), high density lipoprotein cholesterol (HDL-C, P = 4.9 × 10(-4)) and body mass index (BMI)-adjusted waist-hip ratio (P = 9.8 × 10(-3)). These findings improve our knowledge of the genetic basis of adiponectin variation, demonstrate the shared allelic architecture for adiponectin with lipids and central obesity and motivate further studies of underlying mechanisms.

Replication of LIN28B SNP association with age of menarche in young Filipino women

BACKGROUND

Age of menarche, or the timing of first menses in girls, is a physiological trait that shows substantial genetic heritability. Earlier age of menarche is associated with increased childhood adiposity and with adult risk of obesity and cardiovascular disease.

OBJECTIVES

We sought to further characterize the single nucleotide polymorphism (SNP) rs7759938 from the menarche locus LIN28B in 827 young Filipino women from the Cebu Longitudinal Health and Nutrition Survey (CLHNS).

METHODS

We tested rs7759938 for additive association with age of menarche and also tested whether childhood adiposity, as measured by body mass index (BMI) at age 8, mediated this relationship.

RESULTS

We observed nominal association of rs7759938 with age of menarche (β = -0.118 years, 95% confidence interval = (-0.216, -0.020), P = 0.019) with an effect direction consistent with the previous report. We also observed suggestive evidence that the effect of the SNP on age of menarche was independent of childhood BMI.

CONCLUSIONS

These data confirm the strongest gene reported in Europeans (LIN28B) as a contributor to age of menarche in an Asian population.

Genetic association with lipids in Filipinos: waist circumference modifies an APOA5 effect on triglyceride levels

Blood levels of lipoprotein cholesterol and triglycerides (TGs) are highly heritable and are major risk factors for cardiovascular disease (CVD). Approximately 100 lipid-associated loci have been identified in populations of European ancestry. We performed a genome-wide association study of lipid traits in 1,782 Filipino women from the Cebu Longitudinal Health and Nutrition Survey, and tested for evidence of interactions with waist circumference. We conducted additional association and interaction analyses in 1,719 of their young adult offspring. Genome-wide significant associations (P < 5 × 10⁻⁸) were detected at APOE for low density lipoprotein cholesterol and total cholesterol, and at APOA5 for TGs. Suggestive associations (P < 10⁻⁶) were detected at GCKR for TGs, and at CETP and TOM1 for high density lipoprotein cholesterol. Our data also supported the existence of allelic heterogeneity at APOA5, CETP, LIPC, and APOE. The secondary signal (Gly185Cys) at APOA5 exhibited a single nucleotide polymorphism (SNP)-by-waist circumference interaction affecting TGs (Pinteraction = 1.6 × 10⁻⁴), manifested by stronger SNP effects as waist circumference increased. These findings provide the first evidence that central obesity may accentuate the effect of the TG-increasing allele of the APOA5 signal, emphasizing that CVD risk could be reduced by central obesity control.

Sex-stratified genome-wide association studies including 270,000 individuals show sexual dimorphism in genetic loci for anthropometric traits

Given the anthropometric differences between men and women and previous evidence of sex-difference in genetic effects, we conducted a genome-wide search for sexually dimorphic associations with height, weight, body mass index, waist circumference, hip circumference, and waist-to-hip-ratio (133,723 individuals) and took forward 348 SNPs into follow-up (additional 137,052 individuals) in a total of 94 studies. Seven loci displayed significant sex-difference (FDR<5%), including four previously established (near GRB14/COBLL1, LYPLAL1/SLC30A10, VEGFA, ADAMTS9) and three novel anthropometric trait loci (near MAP3K1, HSD17B4, PPARG), all of which were genome-wide significant in women (P<5×10(-8)), but not in men. Sex-differences were apparent only for waist phenotypes, not for height, weight, BMI, or hip circumference. Moreover, we found no evidence for genetic effects with opposite directions in men versus women. The PPARG locus is of specific interest due to its role in diabetes genetics and therapy. Our results demonstrate the value of sex-specific GWAS to unravel the sexually dimorphic genetic underpinning of complex traits.

Genome-wide meta-analysis identifies 11 new loci for anthropometric traits and provides insights into genetic architecture

Approaches exploiting trait distribution extremes may be used to identify loci associated with common traits, but it is unknown whether these loci are generalizable to the broader population. In a genome-wide search for loci associated with the upper versus the lower 5th percentiles of body mass index, height and waist-to-hip ratio, as well as clinical classes of obesity, including up to 263,407 individuals of European ancestry, we identified 4 new loci (IGFBP4, H6PD, RSRC1 and PPP2R2A) influencing height detected in the distribution tails and 7 new loci (HNF4G, RPTOR, GNAT2, MRPS33P4, ADCY9, HS6ST3 and ZZZ3) for clinical classes of obesity. Further, we find a large overlap in genetic structure and the distribution of variants between traits based on extremes and the general population and little etiological heterogeneity between obesity subgroups.

Trans-ethnic fine-mapping of lipid loci identifies population-specific signals and allelic heterogeneity that increases the trait variance explained

Genome-wide association studies (GWAS) have identified ∼100 loci associated with blood lipid levels, but much of the trait heritability remains unexplained, and at most loci the identities of the trait-influencing variants remain unknown. We conducted a trans-ethnic fine-mapping study at 18, 22, and 18 GWAS loci on the Metabochip for their association with triglycerides (TG), high-density lipoprotein cholesterol (HDL-C), and low-density lipoprotein cholesterol (LDL-C), respectively, in individuals of African American (n = 6,832), East Asian (n = 9,449), and European (n = 10,829) ancestry. We aimed to identify the variants with strongest association at each locus, identify additional and population-specific signals, refine association signals, and assess the relative significance of previously described functional variants. Among the 58 loci, 33 exhibited evidence of association at P<1×10⁻⁴ in at least one ancestry group. Sequential conditional analyses revealed that ten, nine, and four loci in African Americans, Europeans, and East Asians, respectively, exhibited two or more signals. At these loci, accounting for all signals led to a 1.3- to 1.8-fold increase in the explained phenotypic variance compared to the strongest signals. Distinct signals across ancestry groups were identified at PCSK9 and APOA5. Trans-ethnic analyses narrowed the signals to smaller sets of variants at GCKR, PPP1R3B, ABO, LCAT, and ABCA1. Of 27 variants reported previously to have functional effects, 74% exhibited the strongest association at the respective signal. In conclusion, trans-ethnic high-density genotyping and analysis confirm the presence of allelic heterogeneity, allow the identification of population-specific variants, and limit the number of candidate SNPs for functional studies.

Lipid traits are heritable, but many of the DNA variants that influence lipid levels remain unknown. In a genomic region, more than one variant may affect gene expression or function, and the frequencies of these variants can differ across populations. Genotyping densely spaced variants in individuals with different ancestries may increase the chance of identifying variants that affect gene expression or function. We analyzed high-density genotyped variants for association with TG, HDL-C, and LDL-C in African Americans, East Asians, and Europeans. At several genomic regions, we provide evidence that two or more variants can influence lipid traits; across loci, these additional signals increase the proportion of trait variation that can be explained by genes. At some association signals shared across populations, combining data from individuals of different ancestries narrowed the set of likely functional variants. At PCSK9 and APOA5, the data suggest that different variants influence trait levels in different populations. Variants previously reported to alter gene expression or function frequently exhibited the strongest association at those signals. The multiple signals and population-specific characteristics of the loci described here may be shared by genetic loci for other complex traits.

A comprehensive SNP and indel imputability database

MOTIVATION

Genotype imputation has become an indispensible step in genome-wide association studies (GWAS). Imputation accuracy, directly influencing downstream analysis, has shown to be improved using re-sequencing-based reference panels; however, this comes at the cost of high computational burden due to the huge number of potentially imputable markers (tens of millions) discovered through sequencing a large number of individuals. Therefore, there is an increasing need for access to imputation quality information without actually conducting imputation. To facilitate this process, we have established a publicly available SNP and indel imputability database, aiming to provide direct access to imputation accuracy information for markers identified by the 1000 Genomes Project across four major populations and covering multiple GWAS genotyping platforms.

RESULTS

SNP and indel imputability information can be retrieved through a user-friendly interface by providing the ID(s) of the desired variant(s) or by specifying the desired genomic region. The query results can be refined by selecting relevant GWAS genotyping platform(s). This is the first database providing variant imputability information specific to each continental group and to each genotyping platform. In Filipino individuals from the Cebu Longitudinal Health and Nutrition Survey, our database can achieve an area under the receiver-operating characteristic curve of 0.97, 0.91, 0.88 and 0.79 for markers with minor allele frequency >5%, 3-5%, 1-3% and 0.5-1%, respectively. Specifically, by filtering out 48.6% of markers (corresponding to a reduction of up to 48.6% in computational costs for actual imputation) based on the imputability information in our database, we can remove 77%, 58%, 51% and 42% of the poorly imputed markers at the cost of only 0.3%, 0.8%, 1.5% and 4.6% of the well-imputed markers with minor allele frequency >5%, 3-5%, 1-3% and 0.5-1%, respectively.

AVAILABILITY

http://www.unc.edu/∼yunmli/imputability.html

Population-specific coding variant underlies genome-wide association with adiponectin level

Adiponectin is a protein hormone that can affect major metabolic processes including glucose regulation and fat metabolism. Our previous genome-wide association (GWA) study of circulating plasma adiponectin levels in Filipino women from the Cebu Longitudinal Health and Nutrition Survey (CLHNS) detected a 100 kb two-SNP haplotype at KNG1-ADIPOQ associated with reduced adiponectin (frequency = 0.050, P = 1.8 × 10(-25)). Subsequent genotyping of CLHNS young adult offspring detected an uncommon variant [minor allele frequency (MAF) = 0.025] located ~800 kb from ADIPOQ that showed strong association with lower adiponectin levels (P = 2.7 × 10(-15), n = 1695) and tagged a subset of KNG1-ADIPOQ haplotype carriers with even lower adiponectin levels. Sequencing of the ADIPOQ-coding region detected variant R221S (MAF = 0.015, P = 2.9 × 10(-69)), which explained 17.1% of the variance in adiponectin levels and largely accounted for the initial GWA signal in Filipinos. R221S was not present in 12 514 Europeans with previously sequenced exons. To explore the mechanism of this substitution, we re-measured adiponectin level in 20 R221S offspring carriers and 20 non-carriers using two alternative antibodies and determined that the presence of R221S resulted in artificially low quantification of adiponectin level using the original immunoassay. These data provide an example of an uncommon variant responsible for a GWA signal and demonstrate that genetic associations with phenotypes measured by antibody-based quantification methods can be affected by uncommon coding SNPs residing in the antibody target region.

Genome-wide association study of anthropometric traits and evidence of interactions with age and study year in Filipino women

Increased values of multiple adiposity-related anthropometric traits are important risk factors for many common complex diseases. We performed a genome-wide association (GWA) study for four quantitative traits related to body size and adiposity (BMI, weight, waist circumference, and height) in a cohort of 1,792 adult Filipino women from the Cebu Longitudinal Health and Nutrition Survey (CLHNS). This is the first GWA study of anthropometric traits in Filipinos, a population experiencing a rapid transition into a more obesogenic environment. In addition to identifying suggestive evidence of additional single-nucleotide polymorphism (SNP) association signals (P < 10(-5)), we replicated (P < 0.05, same direction of additive effect) associations previously reported in European populations of both BMI and weight with MC4R and FTO, of BMI with BDNF, and of height with EFEMP1, ZBTB38, and NPPC, but none with waist circumference. We also replicated loci reported in Japanese or Korean populations as associated with BMI (OTOL1) and height (HIST1H1PS2, C14orf145, GPC5). A difference in local linkage disequilibrium (LD) between European and Asian populations suggests a narrowed association region for BDNF, while still including a proposed functional nonsynonymous amino acid substitution variant (rs6265, Val66Met). Finally, we observed significant evidence (P < 0.0042) for age-by-genotype interactions influencing BMI for rs17782313 (MC4R) and rs9939609 (FTO), and for a study year-by-genotype interaction for rs4923461 (BDNF). Our results show that several genetic risk factors are associated with anthropometric traits in Filipinos and provide further insight into the effects of BDNF, FTO, and MC4R on BMI.

Genome-wide association study for adiponectin levels in Filipino women identifies CDH13 and a novel uncommon haplotype at KNG1-ADIPOQ

Adiponectin is an adipocyte-secreted protein involved in a variety of metabolic processes, including glucose regulation and fatty acid catabolism. We conducted a genome-wide association study to investigate the genetic loci associated with plasma adiponectin in 1776 unrelated Filipino women from the Cebu Longitudinal Health and Nutrition Survey (CLHNS). Our strongest signal for adiponectin mapped to the gene CDH13 (rs3865188, P ≤ 7.2 × 10(-16)), which encodes a receptor for high-molecular-weight forms of adiponectin. Strong association was also detected near the ADIPOQ gene (rs864265, P = 3.8 × 10(-9)) and at a novel signal 100 kb upstream near KNG1 (rs11924390, P = 7.6 × 10(-7)). All three signals were also observed in 1774 young adult CLHNS offspring and in combined analysis including all 3550 mothers and offspring samples (all P ≤ 1.6 × 10(-9)). An uncommon haplotype of rs11924390 and rs864265 (haplotype frequency = 0.050) was strongly associated with lower adiponectin compared with the most common C-G haplotype in both CLHNS mothers (P = 1.8 × 10(-25)) and offspring (P = 8.7 × 10(-32)). Comprehensive imputation of 2653 SNPs in a 2 Mb region using as reference combined CHB, JPT and CEU haplotypes from the 1000 Genomes Project revealed no variants that perfectly tagged this haplotype. Our findings provide the first genome-wide significant evidence of association with plasma adiponectin at the CDH13 locus and identify a novel uncommon KNG1-ADIPOQ haplotype strongly associated with adiponectin levels in Filipinos.

Biological, clinical and population relevance of 95 loci for blood lipids

Plasma concentrations of total cholesterol, low-density lipoprotein cholesterol, high-density lipoprotein cholesterol and triglycerides are among the most important risk factors for coronary artery disease (CAD) and are targets for therapeutic intervention. We screened the genome for common variants associated with plasma lipids in >100,000 individuals of European ancestry. Here we report 95 significantly associated loci (P < 5 x 10(-8)), with 59 showing genome-wide significant association with lipid traits for the first time. The newly reported associations include single nucleotide polymorphisms (SNPs) near known lipid regulators (for example, CYP7A1, NPC1L1 and SCARB1) as well as in scores of loci not previously implicated in lipoprotein metabolism. The 95 loci contribute not only to normal variation in lipid traits but also to extreme lipid phenotypes and have an impact on lipid traits in three non-European populations (East Asians, South Asians and African Americans). Our results identify several novel loci associated with plasma lipids that are also associated with CAD. Finally, we validated three of the novel genes-GALNT2, PPP1R3B and TTC39B-with experiments in mouse models. Taken together, our findings provide the foundation to develop a broader biological understanding of lipoprotein metabolism and to identify new therapeutic opportunities for the prevention of CAD.

Genome-wide association study of homocysteine levels in Filipinos provides evidence for CPS1 in women and a stronger MTHFR effect in young adults

Plasma homocysteine (Hcy) level is associated with cardiovascular disease and may play an etiologic role in vascular damage, a precursor for atherosclerosis. We performed a genome-wide association study for Hcy in 1786 unrelated Filipino women from the Cebu Longitudinal Health and Nutrition Survey (CLHNS). The most strongly associated single-nucleotide polymorphism (SNP) (rs7422339, P = 4.7 x 10(-13)) encodes Thr1405Asn in the gene CPS1 and explained 3.0% of variation in the Hcy level. The widely studied MTHFR C677T SNP (rs1801133) was also highly significant (P = 8.7 x 10(-10)) and explained 1.6% of the trait variation. We also genotyped these two SNPs in 1679 CLHNS young adult offspring. The MTHFR C677T SNP was strongly associated with Hcy (P = 1.9 x 10(-26)) and explained approximately 5.1% of the variation in the offspring. In contrast, the CPS1 variant was significant only in females (P = 0.11 in all; P = 0.0087 in females). Combined analysis of all samples confirmed that the MTHFR variant was more strongly associated with Hcy in the offspring (interaction P = 1.2 x 10(-5)). Furthermore, although there was evidence for a positive synergistic effect between the CPS1 and MTHFR SNPs in the offspring (interaction P = 0.0046), there was no significant evidence for an interaction in the mothers (P = 0.55). These data confirm a recent finding that CPS1 is a locus influencing Hcy levels in women and suggest that genetic effects on Hcy may differ across developmental stages.