Meta-analysis of genome-wide association studies in East Asian-ancestry populations identifies four new loci for body mass index

A variant in KCNQ1 gene predicts metabolic syndrome among northern urban Han Chinese women

BACKGROUND

Previous studies have reported that the potassium voltage-gated channel subfamily Q member 1 (KCNQ1) gene is associated with diabetes in both European and Asian population. This study aims to find a predictable single nucleotide polymorphism (SNP) to predict the risk of metabolic syndrome (MetS) through investigating the association of SNP in KCNQ1 gene with MetS in Han Chinese women of northern urban area.

METHODS

Six SNPs were selected and genotyped in 1381 unrelated women aged 21 and above, who have had physical check-up in Shandong Provincial Qianfoshan Hospital. Cox proportional model was conducted to access the association between SNPs and MetS.

RESULTS

Sixty one women developed MetS between 2010 and 2015 during the 3055 person-year of follow-up. The cumulative incidence density was 19.964/1000 person-year. The SNP rs163182 was associated with MetS both in the additive genetic model (RR = 1.658, 95% CI: 1.144-2.402) and in the recessive genetic model (RR = 2.461, 95% CI: 1.347-4.496). It remained significant after adjustment. This relationship was also observed in MetS components (BMI and SBP).

CONCLUSION

A novel association between rs163182 and MetS was found in this study, which can predict the occurrence of MetS among northern urban Han Chinese women. More investigations are needed to be done to assess the possible pathway in which KCNQ1 gene affects MetS.

A Large Multiethnic Genome-Wide Association Study of Adult Body Mass Index Identifies Novel Loci

Body mass index (BMI), a proxy measure for obesity, is determined by both environmental (including ethnicity, age, and sex) and genetic factors, with > 400 BMI-associated loci identified to date. However, the impact, interplay, and underlying biological mechanisms among BMI, environment, genetics, and ancestry are not completely understood. To further examine these relationships, we utilized 427,509 calendar year-averaged BMI measurements from 100,418 adults from the single large multiethnic Genetic Epidemiology Research on Adult Health and Aging (GERA) cohort. We observed substantial independent ancestry and nationality differences, including ancestry principal component interactions and nonlinear effects. To increase the list of BMI-associated variants before assessing other differences, we conducted a genome-wide association study (GWAS) in GERA, with replication in the Genetic Investigation of Anthropomorphic Traits (GIANT) consortium combined with the UK Biobank (UKB), followed by GWAS in GERA combined with GIANT, with replication in the UKB. We discovered 30 novel independent BMI loci (P < 5.0 × 10^-8) that replicated. We then assessed the proportion of BMI variance explained by sex in the UKB using previously identified loci compared to previously and newly identified loci and found slight increases: from 3.0 to 3.3% for males and from 2.7 to 3.0% for females. Further, the variance explained by previously and newly identified variants decreased with increasing age in the GERA and UKB cohorts, echoed in the variance explained by the entire genome, which also showed gene-age interaction effects. Finally, we conducted a tissue expression QTL enrichment analysis, which revealed that GWAS BMI-associated variants were enriched in the cerebellum, consistent with prior work in humans and mice.

Birthweight, time-varying adiposity growth and early menarche in girls: A Mendelian randomisation and mediation analysis

OBJECTIVE

To explore the causal effect of time-varying z-BMI growth on early menarche using Mendelian randomisation (MR); to identify critical adiposity predictors of early menarche; to compare the effects of birthweight and time-varying z-BMI growth as mediators of the path from genes to early menarche using mediation analysis.

METHODS

We used data from the Taiwan Children Health Study with 21 obesity-related single-nucleotide polymorphisms (SNPs) to yield genetic (instrumental variable)IVs for adiposity. Children with available data on genotyping, birthweight, adiposity, and menarcheal age were included.

RESULTS

In MR analyses, results based on the time-varying z-BMI growth show more statistical power and capture more information of adiposity growth (p=0.01) than those based on single point z-BMI (p=0.02). Among adiposity measures, critical predictors of early menarche are fat free mass (RR=1.33, 95% CI 1.07-1.65) and waist/height ratio (RR=1.27, 95% CI 1.03-1.56). Other potential predictors of early menarche are sum of skinfold (RR=1.24, 95% CI 1.03-1.48) and total body fat (RR=1.20, 95% CI 1.05-1.38). In both one-mediation and multi-mediation analyses, time-varying z-BMI growth in the prepubertal years plays a crucial mediator in the pathway from the genes to early menarche.

CONCLUSIONS

This study discovered that greater prepubertal adiposity growth is a crucial mediator in the path from genes to early menarche. For girls with genes positively associated with obesity; and/or of lower birthweight, a strategy to prevent childhood adiposity should be implemented in order to avoid early menarche development.

Trans-Ethnic Polygenic Analysis Supports Genetic Overlaps of Lumbar Disc Degeneration With Height, Body Mass Index, and Bone Mineral Density

Lumbar disc degeneration (LDD) is age-related break-down in the fibrocartilaginous joints between lumbar vertebrae. It is a major cause of low back pain and is conventionally assessed by magnetic resonance imaging (MRI). Like most other complex traits, LDD is likely polygenic and influenced by both genetic and environmental factors. However, genome-wide association studies (GWASs) of LDD have uncovered few susceptibility loci due to the limited sample size. Previous epidemiology studies of LDD also reported multiple heritable risk factors, including height, body mass index (BMI), bone mineral density (BMD), lipid levels, etc. Genetics can help elucidate causality between traits and suggest loci with pleiotropic effects. One such approach is polygenic score (PGS) which summarizes the effect of multiple variants by the summation of alleles weighted by estimated effects from GWAS. To investigate genetic overlaps of LDD and related heritable risk factors, we calculated the PGS of height, BMI, BMD and lipid levels in a Chinese population-based cohort with spine MRI examination and a Japanese case-control cohort of lumbar disc herniation (LDH) requiring surgery. Because most large-scale GWASs were done in European populations, PGS of corresponding traits were created using weights from European GWASs. We calibrated their prediction performance in independent Chinese samples, then tested associations with MRI-derived LDD scores and LDH affection status. The PGS of height, BMI, BMD and lipid levels were strongly associated with respective phenotypes in Chinese, but phenotype variances explained were lower than in Europeans which would reduce the power to detect genetic overlaps. Despite of this, the PGS of BMI and lumbar spine BMD were significantly associated with LDD scores; and the PGS of height was associated with the increased the liability of LDH. Furthermore, linkage disequilibrium score regression suggested that, osteoarthritis, another degenerative disorder that shares common features with LDD, also showed genetic correlations with height, BMI and BMD. The findings suggest a common key contribution of biomechanical stress to the pathogenesis of LDD and will direct the future search for pleiotropic genes.

Causal relationships between adiposity and childhood asthma: bi-directional Mendelian Randomization analysis

BACKGROUND/OBJECTIVES

Obesity and asthma are common chronic diseases and have been reported to be mutually causative. We investigated the causal direction of the relationship between adiposity and asthma using genetic markers as instrumental variables (IVs) in bi-directional Mendelian randomization (MR) analysis.

SUBJECTS/METHODS

We used data from the Taiwan Children Health Study with 24 body mass index (BMI)-single-nucleotide polymorphisms (SNPs, combined into a weighted allelic score) and 16 asthma-SNPs (combined into two weighted allelic scores, separately for asthma inflammatory and antioxidative genes) to yield genetic IVs for adiposity and asthma, respectively.

RESULTS

The weighted allele score for BMI was strongly associated with adiposity (p = 2 × 10^-16) and active asthma (p = 0.03). The two-stage least square regression risk ratio (RR) for the effect of BMI on asthma was 1.04 (95% confidence interval: 1.00-1.07, p = 0.03). Although the weighted asthma genetic scores were significantly associated with asthma (p = 8.4 × 10^-3), no association was seen for genetically instrumented asthma with BMI using MR. Central obesity was the most accurate predictor of asthma. Adiposity showed higher causal effects on asthma in boys and children with non-atopic asthma. Sensitivity analysis for MR revealed no directional genetic pleiotropy effects. The causal effect RRs of BMI on asthma were 1.04, 1.08, and 1.03 for inverse-variance weighted, MR-Egger regression (slope), and weighted median methods, respectively, all in accordance with the MR estimates.

CONCLUSIONS

High adiposity may lead to asthma, whereas the effects of asthma on adiposity accumulation are likely to be small.

Comprehensive review and annotation of susceptibility SNPs associated with obesity-related traits

We aimed to summarize the results of genetic association studies for obesity and provide a comprehensive annotation of all susceptibility single nucleotide polymorphisms (SNPs). A total of 72 studies were summarized, resulting in 90,361 susceptibility SNPs (738 index SNPs and 89,623 linkage disequilibrium SNPs). Over 90% of the susceptibility SNPs are located in non-coding regions, and it is challenging to understand their functional significance. Therefore, we annotated these SNPs by using various functional databases. We identified 24,623 functional SNPs, including 4 nonsense SNPs, 479 missense SNPs, 399 untranslated region SNPs which might affect microRNA binding, 262 promoter and 5,492 enhancer SNPs which might affect transcription factor binding, 7 splicing sites, 76 SNPs which might affect gene methylation levels, 1,839 SNPs under natural selection and 17,351 SNPs which might modify histone binding. Expression quantitative trait loci analyses for functional SNPs identified 98 target genes, including 69 protein coding genes, 27 long non-coding RNAs and 3 processed transcripts. The percentage of protein coding genes that could be correlated with obesity-related pathways directly or through gene-gene interaction is 75.36 (52/69). Our results may serve as an encyclopaedia of obesity susceptibility SNPs and offer guide for functional experiments.

Identification of rs7350481 at chromosome 11q23.3 as a novel susceptibility locus for metabolic syndrome in Japanese individuals by an exome-wide association study

We have performed exome-wide association studies to identify genetic variants that influence body mass index or confer susceptibility to obesity or metabolic syndrome in Japanese. The exome-wide association study for body mass index included 12,890 subjects, and those for obesity and metabolic syndrome included 12,968 subjects (3954 individuals with obesity, 9014 controls) and 6817 subjects (3998 individuals with MetS, 2819 controls), respectively. Exome-wide association studies were performed with Illumina HumanExome-12 DNA Analysis BeadChip or Infinium Exome-24 BeadChip arrays. The relation of genotypes of single nucleotide polymorphisms to body mass index was examined by linear regression analysis, and that of allele frequencies of single nucleotide polymorphisms to obesity or metabolic syndrome was evaluated with Fisher's exact test. The exome-wide association studies identified six, 11, and 40 single nucleotide polymorphisms as being significantly associated with body mass index, obesity (P <1.21 × 10⁻⁶), or metabolic syndrome (P <1.20 × 10⁻⁶), respectively. Subsequent multivariable logistic regression analysis with adjustment for age and sex revealed that three and five single nucleotide polymorphisms were related (P < 0.05) to obesity or metabolic syndrome, respectively, with one of these latter polymorphisms—rs7350481 (C/T) at chromosome 11q23.3—also being significantly (P < 3.13 × 10⁻⁴) associated with metabolic syndrome. The polymorphism rs7350481 may thus be a novel susceptibility locus for metabolic syndrome in Japanese. In addition, single nucleotide polymorphisms in three genes (CROT, TSC1, RIN3) and at four loci (ANKK1, ZNF804B, CSRNP3, 17p11.2) were implicated as candidate determinants of obesity and metabolic syndrome, respectively.

A Common Variation in the Caveolin 1 Gene Is Associated with High Serum Triglycerides and Metabolic Syndrome in an Admixed Latin American Population

Background: The caveolin 1 (CAV1) gene has been associated with metabolic traits in animal models and human cohorts. Recently, a prevalent variant in CAV1 has been found to be related to metabolic syndrome in Hispanics living in North America. Since Hispanics represent an admixed population at high risk for cardiovascular diseases, in this study a Latin American population with a similar genetic background was assessed.

Objective: To analyze a genetic association between CAV1 and metabolic traits in an admixed Latin American population.

Methods: A cross-sectional study was carried out with adults from the Colombian Caribbean Coast, selected in urban clusters and work places through a stratified sampling to include diverse ages and socioeconomic groups. Blood pressure and waist circumference were registered. Serum concentrations of glucose, triglycerides, and high-density lipoprotein cholesterol were measured from an 8-hr fasting whole-blood sample. Two previously analyzed CAV1 single nucleotide polymorphisms were genotyped (rs926198 and rs11773845). A logistic regression model was applied to estimate the associations. An admixture adjustment was performed through a Bayesian model.

Results: A total of 605 subjects were included. rs11773845 was associated with hypertriglyceridemia [odds ratio (OR) = 1.33, p = 0.001] and the metabolic syndrome (OR = 1.53, p = 0.02). When admixture adjustment was performed these genetic associations preserved their statistical significance. There were no significant associations between rs926198 and metabolic traits.

Conclusions: The CAV1 variation rs11773845 was found to be consistently associated with high serum triglycerides and the metabolic syndrome. This is the first report of a relationship between CAV1 variants and serum triglycerides in Latin America.

Identification of polymorphisms in 12q24.1, ACAD10, and BRAP as novel genetic determinants of blood pressure in Japanese by exome-wide association studies

We performed exome-wide association studies to identify genetic variants that influence systolic or diastolic blood pressure or confer susceptibility to hypertension in Japanese. The exome-wide association studies were performed with the use of Illumina HumanExome-12 DNA Analysis BeadChip or Infinium Exome-24 BeadChip arrays and with 14,678 subjects, including 8215 individuals with hypertension and 6463 controls. The relation of genotypes of 41,843 single nucleotide polymorphisms to systolic or diastolic blood pressure was examined by linear regression analysis. After Bonferroni's correction, 44 and eight polymorphisms were significantly (P < 1.19 × 10⁻⁶) associated with systolic or diastolic blood pressure, respectively, with six polymorphisms (rs12229654, rs671, rs11066015, rs2074356, rs3782886, rs11066280) being associated with both systolic and diastolic blood pressure. Examination of the relation of allele frequencies to hypertension with Fisher's exact test revealed that 100 of the 41,843 single nucleotide polymorphisms were significantly (P < 1.19 × 10⁻⁶) associated with hypertension. Subsequent multivariable logistic regression analysis with adjustment for age and sex showed that five polymorphisms (rs150854849, rs202069030, rs139012426, rs12229654, rs76974938) were significantly (P < 1.25 × 10⁻⁴) associated with hypertension. The polymorphism rs12229654 was thus associated with both systolic and diastolic blood pressure and with hypertension. Six polymorphisms (rs12229654 at 12q24.1, rs671 of ALDH2, rs11066015 of ACAD10, rs2074356 and rs11066280 of HECTD4, and rs3782886 of BRAP) were found to be associated with both systolic and diastolic blood pressure, with those at 12q24.1 or in ACAD10 or BRAP being novel determinants of blood pressure in Japanese.

Genetic factors associated with risk of metabolic syndrome and hepatocellular carcinoma

Although the metabolic syndrome is a commonplace topic, its potential threats to public health is a problem that cannot be neglected. As the living conditions improved significantly over the past few years, the morbidity of metabolic syndrome has also steadily risen, and the onset age is becoming younger. The hepatocellular carcinoma (HCC), is one of the most prevalent life-threatening human cancers worldwide, incidence of which is also on the rise, gradually occupied the top of the list associated with metabolic syndrome related complication. Despite the advanced improvement of HCC management, the lifestyle, environmental factors, obesity, hepatitis B virus (HBV) infection have been recognized as risk factors for the development of liver cancer. In recent years, genetic studies, especially the genome-wide association studies (GWASs) were widely performed, a new era of the human genome research was created, which has significantly promoted the study of complex disease genetics. These progresses have contributed to the discovery of abundant number of genomic loci convincingly linked with complex metabolic feature and HCC. In this review, we briefly summarize the association between metabolic syndrome and HCC, focusing on the genetic factors contributed to metabolic syndrome and HCC.

Supporting precision medicine by data mining across multi-disciplines: an integrative approach for generating comprehensive linkages between single nucleotide variants (SNVs) and drug-binding sites

Motivation

Genetic variants in drug targets and metabolizing enzymes often have important functional implications, including altering the efficacy and toxicity of drugs. Identifying single nucleotide variants (SNVs) that contribute to differences in drug response and understanding their underlying mechanisms are fundamental to successful implementation of the precision medicine model. This work reports an effort to collect, classify and analyze SNVs that may affect the optimal response to currently approved drugs.

Results

An integrated approach was taken involving data mining across multiple information resources including databases containing drugs, drug targets, chemical structures, protein-ligand structure complexes, genetic and clinical variations as well as protein sequence alignment tools. We obtained 2640 SNVs of interest, most of which occur rarely in populations (minor allele frequency < 0.01). Clinical significance of only 9.56% of the SNVs is known in ClinVar, although 79.02% are predicted as deleterious. The examples here demonstrate that even if the mapped SNVs predicted as deleterious may not result in significant structural modifications, they can plausibly modify the protein-drug interactions, affecting selectivity and drug-binding affinity. Our analysis identifies potentially deleterious SNVs present on drug-binding residues that are relevant for further studies in the context of precision medicine.

Availability and Implementation

Data are available from Supplementary information file.

Contact

yanli.wang@nih.gov.

Supplementary information

Supplementary Tables S1-S5 are available at Bioinformatics online.

Obesity genetics: insights from the Pakistani population

The Pakistani population is extensively diverse, indicating a genetic admixture of European and Central/West Asian migrants with indigenous South Asian gene pools. Pakistanis are organized in different ethnicities/castes based on cultural, linguistic and geographical origin. While Pakistan is facing a rapid nutritional transition, the rising prevalence of obesity is driving a growing burden of health complications and mortality. This represents a unique opportunity for the research community to study the interplay between obesogenic environmental changes and obesity predisposing genes in the time frame of one generation. This review recapitulates the ancestral origins of Pakistani population, the societal determinants of the rise in obesity and its governmental management. We describe the contribution of syndromic, monogenic non-syndromic and polygenic obesity genes identified in the Pakistani population. We then discuss the utility of gene identification approaches based on large consanguineous families and original gene × environment interaction study designs in discovering new obesity genes and causal pathways. Elucidation of the genetic basis of obesity in the Pakistani population may result in improved methods of obesity prevention and treatment globally.

Trans-ethnic analysis of metabochip data identifies two new loci associated with BMI

OBJECTIVE

Body mass index (BMI) is commonly used to assess obesity, which is associated with numerous diseases and negative health outcomes. BMI has been shown to be a heritable, polygenic trait, with close to 100 loci previously identified and replicated in multiple populations. We aim to replicate known BMI loci and identify novel associations in a trans-ethnic study population.

SUBJECTS

Using eligible participants from the Population Architecture using Genomics and Epidemiology consortium, we conducted a trans-ethnic meta-analysis of 102 514 African Americans, Hispanics, Asian/Native Hawaiian, Native Americans and European Americans. Participants were genotyped on over 200 000 SNPs on the Illumina Metabochip custom array, or imputed into the 1000 Genomes Project (Phase I). Linear regression of the natural log of BMI, adjusting for age, sex, study site (if applicable), and ancestry principal components, was conducted for each race/ethnicity within each study cohort. Race/ethnicity-specific, and combined meta-analyses used fixed-effects models.

RESULTS

We replicated 15 of 21 BMI loci included on the Metabochip, and identified two novel BMI loci at 1q41 (rs2820436) and 2q31.1 (rs10930502) at the Metabochip-wide significance threshold (P<2.5 × 10-7). Bioinformatic functional investigation of SNPs at these loci suggests a possible impact on pathways that regulate metabolism and adipose tissue.

CONCLUSION

Conducting studies in genetically diverse populations continues to be a valuable strategy for replicating known loci and uncovering novel BMI associations.

Genetics of obesity: what genetic association studies have taught us about the biology of obesity and its complications

Genome-wide association studies (GWAS) for BMI, waist-to-hip ratio, and other adiposity traits have identified more than 300 single-nucleotide polymorphisms (SNPs). Although there is reason to hope that these discoveries will eventually lead to new preventive and therapeutic agents for obesity, this will take time because such developments require detailed mechanistic understanding of how an SNP influences phenotype (and this information is largely unavailable). Fortunately, absence of functional information has not prevented GWAS findings from providing insights into the biology of obesity. Genes near loci regulating total body mass are enriched for expression in the CNS, whereas genes for fat distribution are enriched in adipose tissue itself. Gene by environment and lifestyle interaction analyses have revealed that our increasingly obesogenic environment might be amplifying genetic risk for obesity, yet those at highest risk could mitigate this risk by increasing physical activity and possibly by avoiding specific dietary components. GWAS findings have also been used in mendelian randomisation analyses probing the causal association between obesity and its many putative complications. In supporting a causal association of obesity with diabetes, coronary heart disease, specific cancers, and other conditions, these analyses have clinical relevance in identifying which outcomes could be preventable through weight loss interventions.

A genome-wide association study in the Japanese population identifies the 12q24 locus for habitual coffee consumption: The J-MICC Study

Coffee is one of the most widely consumed beverages worldwide, and its role in human health has received much attention. Although genome-wide association studies (GWASs) have investigated genetic variants associated with coffee consumption in European populations, no such study has yet been conducted in an Asian population. Here, we conducted a GWAS to identify common genetic variations that affected coffee consumption in a Japanese population of 11,261 participants recruited as a part of the Japan Multi-Institutional Collaborative Cohort (J-MICC) study. Coffee consumption was collected using a self-administered questionnaire, and converted from categories to cups/day. In the discovery stage (n = 6,312), we found 2 independent loci (12q24.12–13 and 5q33.3) that met suggestive significance (P < 1 × 10⁻⁶). In the replication stage (n = 4,949), the lead variant for the 12q24.12–13 locus (rs2074356) was significantly associated with habitual coffee consumption (P = 2.2 × 10⁻⁶), whereas the lead variant for the 5q33.3 locus (rs1957553) was not (P = 0.53). A meta-analysis of the discovery and replication populations, and the combined analysis using all subjects, revealed that rs2074356 achieved genome-wide significance (P = 2.2 × 10⁻¹⁶ for a meta-analysis). These findings indicate that the 12q24.12-13 locus is associated with coffee consumption among a Japanese population.

Nutrient Sensing, Signaling and Ageing: The Role of IGF-1 and mTOR in Ageing and Age-Related Disease

Nutrient signaling through insulin/IGF-1 was the first pathway demonstrated to regulate ageing and age-related disease in model organisms. Pharmacological or dietary interventions targeting nutrient signaling pathways have been shown to robustly attenuate ageing in many organisms. Caloric restriction, the most widely studied longevity promoting intervention, works through multiple nutrient signaling pathways, while inhibition of mTOR through treatment with rapamycin reproducibly delays ageing and disease through specific inhibition of the mTOR complexes. Although the benefits of reduced insulin/IGF-1 in lifespan and health are well documented in model organisms, defining the precise role of the IGF-1 in human ageing and age-related disease has proven more difficult. Association studies provide some insight but also reveal paradoxes. Low serum IGF-1 predicts longevity, but IGF-1 decreases with age and IGF-1 therapy benefits some of age-related pathologies. Circulating IGF-1 has been associated both positively and negatively with risk of age-related diseases in humans, and in some cases both activation and inhibition of IGF-1 signaling have provided benefit in animal models of the same diseases. Interventions designed modulate the nutrient sensing signaling pathways positively or negatively are already available for clinical use, highlighting the need for a clear understanding of the role of nutrient signaling in ageing and age-related disease. This chapter examines data from model organisms and human genetic association studies, with a special emphasis on IGF-1 and mTOR, and discusses potential models for resolving the paradoxes surrounding IGF-1 data.

Ethnic and population differences in the genetic predisposition to human obesity

Obesity rates have escalated to the point of a global pandemic with varying prevalence across ethnic groups. These differences are partially explained by lifestyle factors in addition to genetic predisposition to obesity. This review provides a comprehensive examination of the ethnic differences in the genetic architecture of obesity. Using examples from evolution, heritability, admixture, monogenic and polygenic studies of obesity, we provide explanations for ethnic differences in the prevalence of obesity. The debate over definitions of race and ethnicity, the advantages and limitations of multi-ethnic studies and future directions of research are also discussed. Multi-ethnic studies have great potential to provide a better understanding of ethnic differences in the prevalence of obesity that may result in more targeted and personalized obesity treatments.

Protein-altering variants associated with body mass index implicate pathways that control energy intake and expenditure in obesity

Genome-wide association studies (GWAS) have identified >250 loci for body mass index (BMI), implicating pathways related to neuronal biology. Most GWAS loci represent clusters of common, noncoding variants from which pinpointing causal genes remains challenging. Here we combined data from 718,734 individuals to discover rare and low-frequency (minor allele frequency (MAF) < 5%) coding variants associated with BMI. We identified 14 coding variants in 13 genes, of which 8 variants were in genes (ZBTB7B, ACHE, RAPGEF3, RAB21, ZFHX3, ENTPD6, ZFR2 and ZNF169) newly implicated in human obesity, 2 variants were in genes (MC4R and KSR2) previously observed to be mutated in extreme obesity and 2 variants were in GIPR. The effect sizes of rare variants are ~10 times larger than those of common variants, with the largest effect observed in carriers of an MC4R mutation introducing a stop codon (p.Tyr35Ter, MAF = 0.01%), who weighed ~7 kg more than non-carriers. Pathway analyses based on the variants associated with BMI confirm enrichment of neuronal genes and provide new evidence for adipocyte and energy expenditure biology, widening the potential of genetically supported therapeutic targets in obesity.

GWAS of the electrocardiographic QT interval in Hispanics/Latinos generalizes previously identified loci and identifies population-specific signals

QT interval prolongation is a heritable risk factor for ventricular arrhythmias and can predispose to sudden death. Most genome-wide association studies (GWAS) of QT were performed in European ancestral populations, leaving other groups uncharacterized. Herein we present the first QT GWAS of Hispanic/Latinos using data on 15,997 participants from four studies. Study-specific summary results of the association between 1000 Genomes Project (1000G) imputed SNPs and electrocardiographically measured QT were combined using fixed-effects meta-analysis. We identified 41 genome-wide significant SNPs that mapped to 13 previously identified QT loci. Conditional analyses distinguished six secondary signals at NOS1AP (n = 2), ATP1B1 (n = 2), SCN5A (n = 1), and KCNQ1 (n = 1). Comparison of linkage disequilibrium patterns between the 13 lead SNPs and six secondary signals with previously reported index SNPs in 1000G super populations suggested that the SCN5A and KCNE1 lead SNPs were potentially novel and population-specific. Finally, of the 42 suggestively associated loci, AJAP1 was suggestively associated with QT in a prior East Asian GWAS; in contrast BVES and CAP2 murine knockouts caused cardiac conduction defects. Our results indicate that whereas the same loci influence QT across populations, population-specific variation exists, motivating future trans-ethnic and ancestrally diverse QT GWAS.

The contribution of recently identified adult BMI risk loci to paediatric obesity in a Singaporean Chinese childhood dataset

INTRODUCTION

Recent genome-wide association studies have identified 103 adult obesity risk loci; however, it is unclear if these findings are relevant to East-Asian childhood body mass index (BMI) levels.

METHODS AND RESULTS

We evaluated for paediatric obesity associations at these risk loci utilizing genome-wide data from Chinese childhood subjects in the Singapore Cohort study Of the Risk factors for Myopia study (N = 1006). A weighted gene-risk score of all adult obesity risk loci in the Singapore Cohort study Of the Risk factors for Myopia study showed strong associations with BMI at age 9 (p-value = 3.40 × 10^-12 ) and 4-year average BMI (age 9 to 12, p-value = 6.67 × 10^-8 ). Directionally consistent nominal associations for 15 index single nucleotide polymorphisms (SNPs) (p-value < 0.05) were observed. Pathway analysis with genes from these 15 replicating loci revealed over-representation for the G-protein-coupled receptor (GPCR)-mediated integration of entero-endocrine signalling pathway exemplified by L-cell (adjusted p-value = 0.018). Evaluations of birth weight to modify the effects of BMI risk SNPs in paediatric obesity did not reveal significant interactions, and these SNPs were generally not associated with birth weight.

CONCLUSIONS

At least some common adult BMI risk variants predispose to paediatric obesity risk in East-Asians.

Population differentiation in allele frequencies of obesity-associated SNPs

BACKGROUND

Obesity is emerging as a global health problem, with more than one-third of the world's adult population being overweight or obese. In this study, we investigated worldwide population differentiation in allele frequencies of obesity-associated SNPs (single nucleotide polymorphisms).

RESULTS

We collected a total of 225 obesity-associated SNPs from a public database. Their population-level allele frequencies were derived based on the genotype data from 1000 Genomes Project (phase 3). We used hypergeometric model to assess whether the effect allele at a given SNP is significantly enriched or depleted in each of the 26 populations surveyed in the 1000 Genomes Project with respect to the overall pooled population. Our results indicate that 195 out of 225 SNPs (86.7%) possess effect alleles significantly enriched or depleted in at least one of the 26 populations. Populations within the same continental group exhibit similar allele enrichment/depletion patterns whereas inter-continental populations show distinct patterns. Among the 225 SNPs, 15 SNPs cluster in the first intron region of the FTO gene, which is a major gene associated with body-mass index (BMI) and fat mass. African populations exhibit much smaller blocks of LD (linkage disequilibrium) among these15 SNPs while European and Asian populations have larger blocks. To estimate the cumulative effect of all variants associated with obesity, we developed the personal composite genetic risk score for obesity. Our results indicate that the East Asian populations have the lowest averages of the composite risk scores, whereas three European populations have the highest averages. In addition, the population-level average of composite genetic risk scores is significantly correlated (R2 = 0.35, P = 0.0060) with obesity prevalence.

CONCLUSIONS

We have detected substantial population differentiation in allele frequencies of obesity-associated SNPs. The results will help elucidate the genetic basis which may contribute to population disparities in obesity prevalence.

Epigenetic reprogramming in metabolic disorders: nutritional factors and beyond

Environmental factors (e.g., malnutrition and physical inactivity) contribute largely to metabolic disorders including obesity, type 2 diabetes, cardiometabolic disease and nonalcoholic fatty liver diseases. The abnormalities in metabolic activity and pathways have been increasingly associated with altered DNA methylation, histone modification and noncoding RNAs, whereas lifestyle interventions targeting diet and physical activity can reverse the epigenetic and metabolic changes. Here we review recent evidence primarily from human studies that links DNA methylation reprogramming to metabolic derangements or improvements, with a focus on cross-tissue (e.g., the liver, skeletal muscle, pancreas, adipose tissue and blood samples) epigenetic markers, mechanistic mediators of the epigenetic reprogramming, and the potential of using epigenetic traits to predict disease risk and intervention response. The challenges in epigenetic studies addressing the mechanisms of metabolic diseases and future directions are also discussed and prospected.

Convergence between biological, behavioural and genetic determinants of obesity

Multiple biological, behavioural and genetic determinants or correlates of obesity have been identified to date. Genome-wide association studies (GWAS) have contributed to the identification of more than 100 obesity-associated genetic variants, but their roles in causal processes leading to obesity remain largely unknown. Most variants are likely to have tissue-specific regulatory roles through joint contributions to biological pathways and networks, through changes in gene expression that influence quantitative traits, or through the regulation of the epigenome. The recent availability of large-scale functional genomics resources provides an opportunity to re-examine obesity GWAS data to begin elucidating the function of genetic variants. Interrogation of knockout mouse phenotype resources provides a further avenue to test for evidence of convergence between genetic variation and biological or behavioural determinants of obesity.

The role of established East Asian obesity-related loci on pediatric leptin levels highlights a neuronal influence on body weight regulation in Chinese children and adolescents: the BCAMS study

Genome-wide association studies have identified multiple variants associated with adult obesity, mostly in European-ancestry populations. We aimed to systematically assess the contribution of key loci, which had been previously shown to be associated in East Asian adults, to childhood obesity, related adipokine profiles and metabolic traits in a Chinese pediatric population. Twelve single-nucleotide polymorphisms (SNPs) plus metabolic profiles and levels of five adipokines (leptin, adiponectin, resistin, fibroblast growth factor 21 and retinol binding protein 4) were evaluated in 3,506 Chinese children and adolescents aged 6-18. After correction for multiple comparisons, six of these SNPs were robustly associated with childhood obesity: FTO-rs1558902 (P=5.6×10⁻⁵), MC4R-rs2331841 (P=4.4×10⁻⁴), GNPDA2-rs16858082 (P = 3.4×10⁻⁴), PCSK1-rs261967 (P = 0.001), SEC16B-rs516636 (P = 0.004) and MAP2K5-rs4776970 (P = 0.004), with odds ratios ranging from 1.211 to 1.421; while ITIH4-rs2535633 and BDNF-rs2030323 yielded nominal association with the same trait (P < 0.05). Moreover, the risk alleles of six SNPs displayed significant (P < 0.004) or nominal (P < 0.05) association with leptin levels, namely at in/near PCSK1, MC4R, FTO, MAP2K5, GNPDA2 and BDNF plus their cumulative genetic score yielded stronger association with increased leptin levels (P = 6.2×10⁻¹¹). Our results reveal that key obesity-associated loci previously reported in Europeans, but also associated with East Asian adults, are also associated with obesity and/or metabolic quantitative traits in Chinese children. These associations coincide with six brain-expressed loci that correlate with leptin levels, thus may point to an important neuronal influence on body weight regulation in the pediatric setting.

Interaction between obesity and the Hypoxia Inducible Factor 3 Alpha Subunit rs3826795 polymorphism in relation with plasma alanine aminotransferase

BACKGROUND

Hypoxia Inducible Factor 3 Alpha Subunit (HIF3A) DNA has been demonstrated to be associated with obesity in the methylation level, and it also has a Body Mass Index (BMI)-independent association with plasma alanine aminotransferase (ALT). However, the relation among obesity, plasma ALT, HIF3A polymorphism and methylation remains unclear. This study aims to identify the association between HIF3A polymorphism and plasma ALT, and further to determine whether the effect of HIF3A polymorphism on ALT could be modified by obesity or mediated by DNA methylation.

METHODS

The HIF3A rs3826795 polymorphism was genotyped in a case-control study including 2030 Chinese children aged 7-18 years (705 obese cases and 1325 non-obese controls). Furthermore, the HIF3A DNA methylation of the peripheral blood was measured in 110 severely obese children and 110 age- and gender- matched normal-weight controls.

RESULTS

There was no overall association between the HIF3A rs3826795 polymorphism and ALT. A significant interaction between obesity and rs3826795 in relation with ALT was found (P _inter = 0.042), with rs3826795 G-allele number elevating ALT significantly only in obese children (β' = 0.075, P = 0.037), but not in non-obese children (β' = -0.009, P = 0.741). Additionally, a mediation effect of HIF3A methylation was found in the association between the HIF3A rs3826795 polymorphism and ALT among obese children (β' = 0.242, P = 0.014).

CONCLUSION

This is the first study to report the interaction between obesity and HIF3A gene in relation with ALT, and also to reveal a mediation effect among the HIF3A polymorphism, methylation and ALT. This study provides new evidence to the function of HIF3A gene, which would be helpful for future risk assessment and personalized treatment of liver diseases.

Association between age at menarche and cardiovascular disease: A systematic review on risk and potential mechanisms

Age at menarche (AAM) has been reported to be associated with the risk of cardiovascular disease (CVD), but the shape of and the mechanisms behind this association remain unclear. We reviewed the data on the association between AAM and different subtypes of CVD, and used shared genetic loci to identify possible mechanisms underlying this association using shared genetic association. We searched the databases of PubMed, Web of Science and Embase through to April 2017. We included articles with any clinically manifest CVD endpoint and for any ethnicity. We identified single nucleotide polymorphisms (SNPs) for AAM in genome-wide association studies (GWAS) in Caucasians through PubMed and HuGE Navigator, and searched whether these SNPs or any of their proxies were associated with any CVD-related trait. Eight studies in Caucasian populations reported an inverse linear relation between AAM and CVD risk, whereas one large study reported a significant U-shaped relation between them. Data from Asian populations were contradictory and inconclusive. In total, 122 AAM SNPs were identified at a genome-wide significance level (p<5×10^-8). Of those, 18 were also associated with various CVD-related traits, primarily body mass index (BMI), obesity, and height. In conclusion, early AAM and possibly also late AAM increase the risk of CVD in Caucasian populations. Weight and height may be part of the mechanism underlying the relation between AAM and CVD risk in Caucasians. Data on other ethnicities are too limited for meaningful analysis and conclusions.

Genetically predicted high body mass index is associated with increased gastric cancer risk

Epidemiological studies have linked body mass index (BMI) with risk of gastrointestinal cancers. However, for gastric cancer, the relationship is more controversial. In particular, it is unclear whether the observed association is due to confounding or bias inherent in conventional observational studies. To investigate whether BMI is causally associated with gastric cancer risk, we applied Mendelian randomization using individual-level data from 2631 gastric cancer cases and 4373 cancer-free controls. We derived a weighted genetic risk score (wGRS) using 37 BMI-associated genetic variants as an instrumental variable. We used logistic regression to estimate odds ratios (OR) and 95% confidence intervals (CI) for associations between genetically predicted BMI and gastric cancer risk. We observed that higher genetically determined BMI was associated with increased gastric cancer risk (per standard deviation (SD) increase in the wGRS: OR=1.07, 95% CI: 1.02-1.13, P=4.94 × 10^-3). Compared with individuals in the bottom tertile of the BMI wGRS, those in the top tertile had 1.14-fold (95% CI: 1.01-1.29) increased risk of developing gastric cancer. Sensitivity analyses using alternative causal inference measures demonstrated consistent association. Our study indicated that genetically high BMI was associated with increased gastric cancer risk, suggesting that high BMI may have a causal role in the etiology of gastric cancer.

The ALDH2 gene rs671 polymorphism is not associated with essential hypertension

Essential hypertension (EH) is a worldwide problem. Acetaldehyde dehydrogenase 2 (ALDH2) gene has been suggested to be correlated with EH. However, the results are inconsistent. This study aimed to investigate the associations of ALDH2 rs671 polymorphism with EH in a Chinese Han population in Shanghai. Genotype of ALDH2 rs671 was analyzed in 1923 EH patients and 1115 control subjects. We found no association between ALDH2 rs671 and EH risk or EH-related quantitative blood chemistry values. Furthermore, a meta-analysis was performed and the summary results from 11220 patients and 8339 control subjects were consistent with our findings. These results indicated that rs671 of ALDH2 may not associate with the risk of EH.

Trans-ethnic fine-mapping of genetic loci for body mass index in the diverse ancestral populations of the Population Architecture using Genomics and Epidemiology (PAGE) Study reveals evidence for multiple signals at established loci

Most body mass index (BMI) genetic loci have been identified in studies of primarily European ancestries. The effect of these loci in other racial/ethnic groups is less clear. Thus, we aimed to characterize the generalizability of 170 established BMI variants, or their proxies, to diverse US populations and trans-ethnically fine-map 36 BMI loci using a sample of >102,000 adults of African, Hispanic/Latino, Asian, European and American Indian/Alaskan Native descent from the Population Architecture using Genomics and Epidemiology Study. We performed linear regression of the natural log of BMI (18.5-70 kg/m2) on the additive single nucleotide polymorphisms (SNPs) at BMI loci on the MetaboChip (Illumina, Inc.), adjusting for age, sex, population stratification, study site, or relatedness. We then performed fixed-effect meta-analyses and a Bayesian trans-ethnic meta-analysis to empirically cluster by allele frequency differences. Finally, we approximated conditional and joint associations to test for the presence of secondary signals. We noted directional consistency with the previously reported risk alleles beyond what would have been expected by chance (binomial p < 0.05). Nearly, a quarter of the previously described BMI index SNPs and 29 of 36 densely-genotyped BMI loci on the MetaboChip replicated/generalized in trans-ethnic analyses. We observed multiple signals at nine loci, including the description of seven loci with novel multiple signals. This study supports the generalization of most common genetic loci to diverse ancestral populations and emphasizes the importance of dense multiethnic genomic data in refining the functional variation at genetic loci of interest and describing several loci with multiple underlying genetic variants.

Recent Advances in Human Genetics and Epigenetics of Adiposity: Pathway to Precision Medicine?

Obesity is a heritable trait that contributes to substantial global morbidity and mortality. Here, we summarize findings from the past decade of genetic and epigenetic research focused on unravelling the underpinnings of adiposity. More than 140 genetic regions now are known to influence adiposity traits. The genetics of general adiposity, as measured by body mass index, and that of abdominal obesity, as measured by waist-to-hip ratio, have distinct biological backgrounds. Gene expression associated with general adiposity is enriched in the nervous system. In contrast, genes associated with abdominal adiposity function in adipose tissue. Recent population-based epigenetic analyses have highlighted additional distinct loci. We discuss how associated genetic variants can lead to understanding causal mechanisms, and to disentangling reverse causation in epigenetic analyses. Discoveries emerging from population genomics are identifying new disease markers and potential novel drug targets to better define and combat obesity and related diseases.

Discovery and fine-mapping of adiposity loci using high density imputation of genome-wide association studies in individuals of African ancestry: African Ancestry Anthropometry Genetics Consortium

Genome-wide association studies (GWAS) have identified >300 loci associated with measures of adiposity including body mass index (BMI) and waist-to-hip ratio (adjusted for BMI, WHR_adjBMI), but few have been identified through screening of the African ancestry genomes. We performed large scale meta-analyses and replications in up to 52,895 individuals for BMI and up to 23,095 individuals for WHR_adjBMI from the African Ancestry Anthropometry Genetics Consortium (AAAGC) using 1000 Genomes phase 1 imputed GWAS to improve coverage of both common and low frequency variants in the low linkage disequilibrium African ancestry genomes. In the sex-combined analyses, we identified one novel locus (TCF7L2/HABP2) for WHR_adjBMI and eight previously established loci at P < 5×10⁻⁸: seven for BMI, and one for WHR_adjBMI in African ancestry individuals. An additional novel locus (SPRYD7/DLEU2) was identified for WHR_adjBMI when combined with European GWAS. In the sex-stratified analyses, we identified three novel loci for BMI (INTS10/LPL and MLC1 in men, IRX4/IRX2 in women) and four for WHR_adjBMI (SSX2IP, CASC8, PDE3B and ZDHHC1/HSD11B2 in women) in individuals of African ancestry or both African and European ancestry. For four of the novel variants, the minor allele frequency was low (<5%). In the trans-ethnic fine mapping of 47 BMI loci and 27 WHR_adjBMI loci that were locus-wide significant (P < 0.05 adjusted for effective number of variants per locus) from the African ancestry sex-combined and sex-stratified analyses, 26 BMI loci and 17 WHR_adjBMI loci contained ≤ 20 variants in the credible sets that jointly account for 99% posterior probability of driving the associations. The lead variants in 13 of these loci had a high probability of being causal. As compared to our previous HapMap imputed GWAS for BMI and WHR_adjBMI including up to 71,412 and 27,350 African ancestry individuals, respectively, our results suggest that 1000 Genomes imputation showed modest improvement in identifying GWAS loci including low frequency variants. Trans-ethnic meta-analyses further improved fine mapping of putative causal variants in loci shared between the African and European ancestry populations.

Genome-wide association studies (GWAS) have identified >300 genetic regions that influence body size and shape as measured by body mass index (BMI) and waist-to-hip ratio (WHR), respectively, but few have been identified in populations of African ancestry. We conducted large scale high coverage GWAS and replication of these traits in 52,895 and 23,095 individuals of African ancestry, respectively, followed by additional replication in European populations. We identified 10 genome-wide significant loci in all individuals, and an additional seven loci by analyzing men and women separately. We combined African and European ancestry GWAS and were able to narrow down 43 out of 74 African ancestry associated genetic regions to contain small number of putative causal variants. Our results highlight the improvement of applying high density genome coverage and combining multiple ancestries in the identification and refinement of location of genetic regions associated with adiposity traits.

Susceptibility loci for metabolic syndrome and metabolic components identified in Han Chinese: a multi-stage genome-wide association study

Metabolic syndrome (MetS), a cluster of metabolic disturbances that increase the risk for cardiovascular disease and diabetes, was because of genetic susceptibility and environmental risk factors. To identify the genetic variants associated with MetS and metabolic components, we conducted a genome-wide association study followed by replications in totally 12,720 participants from the north, north-eastern and eastern China. In combined analyses, independent of the top known signal at rs651821 on APOA5, we newly identified a secondary triglyceride-associated signal at rs180326 on BUD13 (P_combined = 2.4 × 10^-8 ). Notably, by an integrated analysis of the genotypes and the serum levels of APOA5, BUD13 and triglyceride, we observed that BUD13 was another potential mediator, besides APOA5, of the association between rs651821 and serum triglyceride. rs671 (ALDH2), an east Asian-specific common variant, was found to be associated with MetS (P_combined = 9.7 × 10^-22 ) in Han Chinese. The effects of rs671 on metabolic components were more prominent in drinkers than in non-drinkers. The replicated loci provided information on the genetic basis and mechanisms of MetS and metabolic components in Han Chinese.

Stably transfected adherent cancer cell models with decreased expression of 5'-nucleotidase cN-II

The 5'-nucleotidase cN-II has been shown to be associated with the sensitivity to nucleoside analogues, the survival of cytarabine treated leukemia patients and to cell proliferation. Due to the lack of relevant cell models for solid tumors, we developed four cell lines with low cN-II expression and characterized them concerning their in vitro sensitivity to cancer drugs and their intracellular nucleotide pools. All four cell models had an important decrease of cN-II expression but did not show modified sensitivity, cell proliferation or nucleotide pools. Our cell models will be important for the study of the role of cN-II in human cancer cells.

The transcriptional architecture of phenotypic dimorphism

The profound differences in gene expression between the sexes are increasingly used to study the molecular basis of sexual dimorphism, sexual selection and sexual conflict. Studies of transcriptional architecture, based on comparisons of gene expression, have also been implemented for a wide variety of other intra-specific polymorphisms. These efforts are based on key assumptions regarding the relationship between transcriptional architecture, phenotypic variation and the target of selection. Some of these assumptions are better supported by available evidence than others. In all cases, the evidence is largely circumstantial, leaving considerable gaps in our understanding of the relationship between transcriptional and phenotypic dimorphism.

The genetic overlap between mood disorders and cardiometabolic diseases: a systematic review of genome wide and candidate gene studies

Meta-analyses of genome-wide association studies (meta-GWASs) and candidate gene studies have identified genetic variants associated with cardiovascular diseases, metabolic diseases and mood disorders. Although previous efforts were successful for individual disease conditions (single disease), limited information exists on shared genetic risk between these disorders. This article presents a detailed review and analysis of cardiometabolic diseases risk (CMD-R) genes that are also associated with mood disorders. First, we reviewed meta-GWASs published until January 2016, for the diseases 'type 2 diabetes, coronary artery disease, hypertension' and/or for the risk factors 'blood pressure, obesity, plasma lipid levels, insulin and glucose related traits'. We then searched the literature for published associations of these CMD-R genes with mood disorders. We considered studies that reported a significant association of at least one of the CMD-R genes and 'depression' or 'depressive disorder' or 'depressive symptoms' or 'bipolar disorder' or 'lithium treatment response in bipolar disorder', or 'serotonin reuptake inhibitors treatment response in major depression'. Our review revealed 24 potential pleiotropic genes that are likely to be shared between mood disorders and CMD-Rs. These genes include MTHFR, CACNA1D, CACNB2, GNAS, ADRB1, NCAN, REST, FTO, POMC, BDNF, CREB, ITIH4, LEP, GSK3B, SLC18A1, TLR4, PPP1R1B, APOE, CRY2, HTR1A, ADRA2A, TCF7L2, MTNR1B and IGF1. A pathway analysis of these genes revealed significant pathways: corticotrophin-releasing hormone signaling, AMPK signaling, cAMP-mediated or G-protein coupled receptor signaling, axonal guidance signaling, serotonin or dopamine receptors signaling, dopamine-DARPP32 feedback in cAMP signaling, circadian rhythm signaling and leptin signaling. Our review provides insights into the shared biological mechanisms of mood disorders and cardiometabolic diseases.

Age at menarche and age at natural menopause in East Asian women: a genome-wide association study

Age at menarche (AM) and age at natural menopause (ANM) are complex traits with a high heritability. Abnormal timing of menarche or menopause is associated with a reduced span of fertility and risk for several age-related diseases including breast, endometrial and ovarian cancer, cardiovascular disease, and osteoporosis. To identify novel genetic loci for AM or ANM in East Asian women and to replicate previously identified loci primarily in women of European ancestry by genome-wide association studies (GWASs), we conducted a two-stage GWAS. Stage I aimed to discover promising novel AM and ANM loci using GWAS data of 8073 women from Shanghai, China. The Stage II replication study used the data from another Chinese GWAS (n = 1230 for AM and n = 1458 for ANM), a Korean GWAS (n = 4215 for AM and n = 1739 for ANM), and de novo genotyping of 2877 additional Chinese women. Previous GWAS-identified loci for AM and ANM were also evaluated. We identified two suggestive menarcheal age loci tagged by rs79195475 at 10q21.3 (beta = -0.118 years, P = 3.4 × 10^-6) and rs1023935 at 4p15.1 (beta = -0.145 years, P = 4.9 × 10^-6) and one menopausal age locus tagged by rs3818134 at 22q12.2 (beta = -0.276 years, P = 8.8 × 10^-6). These suggestive loci warrant a further validation in independent populations. Although limited by low statistical power, we replicated 19 of the 98 menarche loci and 5 of the 20 menopause loci previously identified in women of European ancestry in East Asian women, suggesting a shared genetic architecture for these two traits across populations.

Joint association analysis of a binary and a quantitative trait in family samples

In recent years, improved genotyping and sequencing technologies have enabled the discovery of new loci associated with various diseases or traits. For instance, by testing the association with each single-nucleotide variant (SNV) separately, genome-wide association studies (GWAS) have achieved tremendous success in identifying SNVs associated with specific traits. However, little is known about the common genetic basis of multiple traits owing to lack of efficient methods. With the use of extended quasi-likelihood, a Wald test has been proposed to perform a bivariate analysis of a continuous and a binary trait in unrelated samples. However, owing to its low computational efficiency, it has not been implemented in real applications to large-scale genetic studies. In this paper, we propose an efficient bivariate robust score test for two traits, one continuous and one binary, based on extended generalized estimating equations. Our approach is applicable to both family-based and unrelated study designs and can be extended to test the association of multiple traits. Our simulation studies demonstrate the type-I error rate of our approach is well controlled in all minor allele frequency (MAF) scenarios, with MAF ranging from 1 to 30%, and the method is more powerful in certain MAF scenarios than univariate testing with correction for multiple testing. Because of the computational advantage of score tests, our approach is readily applicable to GWAS or sequencing studies. Finally, we present a real application to uncover genetic variants associated with body mass index and type-2 diabetes in the Framingham Heart Study.

A GWAS follow-up of obesity-related SNPs in SYPL2 reveals sex-specific association with hip circumference

Objective

A novel single‐nucleotide polymorphism (SNP) associated with morbid obesity was recently identified by exome sequencing. The purpose of this study was to follow up this low‐frequency coding SNP located within the SYPL2 locus and associated with body mass index in order to reveal novel associations with obesity‐related traits.

Methods

The body mass index‐associated SNP (rs62623713 A>G [chr1:109476817/hg19]) and two tagging SNPs within the SYPL2 locus, rs9661614 T>C (chr1:109479215) and rs485660 G>A (chr1:109480810), were genotyped in the obesity (n = 3,017) and the infogene (n = 676) cohorts, which were further combined, leading to a larger cohort of 3,693 individuals. Association testing was performed by general linear models in the obesity cohort and validated by joint analysis in the combined cohort.

Results

rs9661614 and rs485660 were significantly associated with hip circumference (HC) in the obesity cohort, with heterozygotes exhibiting a significantly lower HC. These results were validated by joint analysis for rs9661614 (false discovery rate [FDR]‐corrected P = 7.5 × 10⁻⁴) and, to a lesser extent, for rs485660 (FDR corrected P = 3.9 × 10⁻²). The association with HC remained significant for rs9661614 when tested independently in women (FDR‐corrected P = 1.7 × 10⁻²), but not for rs485660 (FDR‐corrected P = 0.2). Both associations were absent in men.

Conclusions

This study reveals strong evidence for a novel association between rs9661614 (T>C) and HC in women, which likely reflects a preferential association of SYPL2 to a gynoid profile of fat distribution. The study findings support a clinical significance of SYPL2 worth considering when assessing risk factors associated with obesity.

High-throughput allele-specific expression across 250 environmental conditions

Gene-by-environment (GxE) interactions determine common disease risk factors and biomedically relevant complex traits. However, quantifying how the environment modulates genetic effects on human quantitative phenotypes presents unique challenges. Environmental covariates are complex and difficult to measure and control at the organismal level, as found in GWAS and epidemiological studies. An alternative approach focuses on the cellular environment using in vitro treatments as a proxy for the organismal environment. These cellular environments simplify the organism-level environmental exposures to provide a tractable influence on subcellular phenotypes, such as gene expression. Expression quantitative trait loci (eQTL) mapping studies identified GxE interactions in response to drug treatment and pathogen exposure. However, eQTL mapping approaches are infeasible for large-scale analysis of multiple cellular environments. Recently, allele-specific expression (ASE) analysis emerged as a powerful tool to identify GxE interactions in gene expression patterns by exploiting naturally occurring environmental exposures. Here we characterized genetic effects on the transcriptional response to 50 treatments in five cell types. We discovered 1455 genes with ASE (FDR < 10%) and 215 genes with GxE interactions. We demonstrated a major role for GxE interactions in complex traits. Genes with a transcriptional response to environmental perturbations showed sevenfold higher odds of being found in GWAS. Additionally, 105 genes that indicated GxE interactions (49%) were identified by GWAS as associated with complex traits. Examples include GIPR–caffeine interaction and obesity and include LAMP3–selenium interaction and Parkinson disease. Our results demonstrate that comprehensive catalogs of GxE interactions are indispensable to thoroughly annotate genes and bridge epidemiological and genome-wide association studies.

Tissue Specificity and Sex-Specific Regulatory Variation Permit the Evolution of Sex-Biased Gene Expression

Genetic correlations between males and females are often thought to constrain the evolution of sexual dimorphism. However, sexually dimorphic traits and the underlying sexually dimorphic gene expression patterns are often rapidly evolving. We explore this apparent paradox by measuring the genetic correlation in gene expression between males and females (Cmf) across broad evolutionary timescales, using two RNA-sequencing data sets spanning multiple populations and multiple species. We find that unbiased genes have higher Cmf than sex-biased genes, consistent with intersexual genetic correlations constraining the evolution of sexual dimorphism. However, we found that highly sex-biased genes (both male and female biased) also had higher tissue specificity, and unbiased genes had greater expression breadth, suggesting that pleiotropy may constrain the breakdown of intersexual genetic correlations. Finally, we show that genes with high Cmf showed some degree of sex-specific changes in gene expression in males and females. Together, our results suggest that genetic correlations between males and females may be less important in constraining the evolution of sex-biased gene expression than pleiotropy. Sex-specific regulatory variation and tissue specificity may resolve the paradox of widespread sex bias within a largely shared genome.

Epigenomic elements analyses for promoters identify ESRRG as a new susceptibility gene for obesity-related traits

OBJECTIVES

With ENCODE epigenomic data and results from published genome-wide association studies (GWASs), we aimed to find regulatory signatures of obesity genes and discover novel susceptibility genes.

METHODS

Obesity genes were obtained from public GWAS databases and their promoters were annotated based on the regulatory element information. Significantly enriched or depleted epigenomic elements in the promoters of obesity genes were evaluated and all human genes were then prioritized according to the existence of the selected elements to predict new candidate genes. Top-ranked genes were subsequently applied to validate their associations with obesity-related traits in three independent in-house GWAS samples.

RESULTS

We identified RAD21 and EZH2 as over-represented, and STAT2 (signal transducer and activator of transcription 2) and IRF3 (interferon regulatory transcription factor 3) as depleted transcription factors. Histone modification of H3K9me3 and chromatin state segmentation of 'poised promoter' and 'repressed' were over-represented. All genes were prioritized and we selected the top five genes for validation at the population level. Combining results from the three GWAS samples, rs7522101 in ESRRG (estrogen-related receptor-γ) remained significantly associated with body mass index after multiple testing corrections (P=7.25 × 10(-5)). It was also associated with β-cell function (P=1.99 × 10(-3)) and fasting glucose level (P<0.05) in the meta-analyses of glucose and insulin-related traits consortium (MAGIC) data set.Cnoclusions:In summary, we identified epigenomic characteristics for obesity genes and suggested ESRRG as a novel obesity-susceptibility gene.

Integrative approaches for large-scale transcriptome-wide association studies

Many genetic variants influence complex traits by modulating gene expression, thus altering the abundance of one or multiple proteins. Here we introduce a powerful strategy that integrates gene expression measurements with summary association statistics from large-scale genome-wide association studies (GWAS) to identify genes whose cis-regulated expression is associated with complex traits. We leverage expression imputation from genetic data to perform a transcriptome-wide association study (TWAS) to identify significant expression-trait associations. We applied our approaches to expression data from blood and adipose tissue measured in ∼ 3,000 individuals overall. We imputed gene expression into GWAS data from over 900,000 phenotype measurements to identify 69 new genes significantly associated with obesity-related traits (BMI, lipids and height). Many of these genes are associated with relevant phenotypes in the Hybrid Mouse Diversity Panel. Our results showcase the power of integrating genotype, gene expression and phenotype to gain insights into the genetic basis of complex traits.

Effects of Obesity Related Genetic Variations on Visceral and Subcutaneous Fat Distribution in a Chinese Population

Genome-wide association studies (GWAS) have uncovered numerous variants associated with body mass index (BMI), waist circumference, and waist-to-hip ratio. Our study aims to investigate how these variants are linked to fat distribution. We genotyped 56 validated variants of BMI, waist circumference, and waist-to-hip ratio in 2958 subjects from Chinese community-based populations and performed linear regression analyses to determine the association with visceral fat area (VFA) and subcutaneous fat area (SFA) imaged by magnetic resonance imaging (MRI). We found rs671 in ALDH2 exhibited the significant associations with VFA and the VFA-SFA ratio in all subjects (P = 9.64 × 10⁻⁵ and 6.54 × 10⁻⁴). rs17782313 near MC4R for VFA and rs4846567 near LYPLAL1 for SFA were found in females only (P = 2.93 × 10⁻⁴ and 0.0015), whereas rs671 in ALDH2 for VFA and the VFA-SFA ratio was restricted to males (P = 1.75 × 10⁻⁸ and 4.43 × 10⁻⁸). Given the robust association of rs671 with alcohol consumption, we next demonstrated the primary effects of rs671 on VFA and the VFA-SFA ratio were restricted to drinkers (P = 1.45 × 10⁻⁴ and 4.65 × 10⁻³). Our data implied that variants of MC4R and LYPLAL1 modulated body fat distribution with sexual dimorphism and that alcohol consumption may mediate the impact of the ALDH2 locus on visceral fat in a Chinese population.

Association of common variants identified by recent genome-wide association studies with obesity in Chinese children: a case-control study

Background

Large-scale genome-wide association studies have identified multiple genetic variants that are associated with elevated body mass index (BMI) or the risk of obesity in Caucasian or Asian populations. We examined whether these variants are individually associated with obesity in Chinese children, and also assessed their cumulative effects and predictive value for obesity risk in Chinese children.

Methods

We genotyped 40 single nucleotide polymorphisms (SNPs) and conducted association analyses for 32/40 SNPs with an estimated minor allele frequency >1 % in 2 030 unrelated Chinese children, including 607 normal-weight, 718 overweight, and 705 obese individuals from two cross-sectional study groups. Logistic regression and linear regression under the additive model were used to examine associations, and the area under the receiver operating characteristic curve (AUC_ROC) was reported as prediction summary.

Results

We identified obesity association for 6 SNPs near SEC16B, RBJ, CDKAL1, TFAP2B, MAP2K5 and FTO (odds ratios (ORs) ranged from 1.19 to 1.41, nominal two-sided P-values < 0.05). Association (Bonferroni corrected) of rs543874 near SEC16B and rs2241423 near MAP2K5 had presumably stronger effects on obesity in Chinese children than in Caucasian populations. Their risk alleles were also associated with BMI standard deviation score (BMI-SDS) variability. We demonstrated the cumulative effects of the 32 SNPs on obesity risk (per risk allele: OR = 1.06, 95 % CI: 1.03-1.11, P = 4.84 × 10^-4) and BMI-SDS (β = 0.04, 95 % CI: 0.02-0.06, P = 3.69 × 10^-7). The difference in AUC_ROC for a model with covariates (age, age square, sex and study group) and the model including covariates and all 32 SNPs was 2.8 % (P = 0.0002).

Conclusion

While six SNPs were individually associated with obesity in Chinese children, the 32 common variants identified by recent GWA studies had cumulative effects and resulted in a limited increase in the AUC_ROC predictive value for childhood obesity.

Electronic supplementary material

The online version of this article (doi:10.1186/s12881-016-0268-4) contains supplementary material, which is available to authorized users.

Genome-wide association studies in East Asians identify new loci for waist-hip ratio and waist circumference

Sixty genetic loci associated with abdominal obesity, measured by waist circumference (WC) and waist-hip ratio (WHR), have been previously identified, primarily from studies conducted in European-ancestry populations. We conducted a meta-analysis of associations of abdominal obesity with approximately 2.5 million single nucleotide polymorphisms (SNPs) among 53,052 (for WC) and 48,312 (for WHR) individuals of Asian descent, and replicated 33 selected SNPs among 3,762 to 17,110 additional individuals. We identified four novel loci near the EFEMP1, ADAMTSL3 , CNPY2, and GNAS genes that were associated with WC after adjustment for body mass index (BMI); two loci near the NID2 and HLA-DRB5 genes associated with WHR after adjustment for BMI, and three loci near the CEP120, TSC22D2, and SLC22A2 genes associated with WC without adjustment for BMI. Functional enrichment analyses revealed enrichment of corticotropin-releasing hormone signaling, GNRH signaling, and/or CDK5 signaling pathways for those newly-identified loci. Our study provides additional insight on genetic contribution to abdominal obesity.