Investigation of Genetic Variation Underlying Central Obesity amongst South Asians
Investigation of Genetic Variation Underlying Central Obesity amongst South Asians
William R. Scott 0 1 2 3
Weihua Zhang 1 2 3
Marie Loh 1 2 3
Sian-Tsung Tan 0 1 3
Benjamin Lehne 1 2 3
Uzma Afzal 1 2 3
Juan Peralta 1 3
Richa Saxena 1 3
Sarju Ralhan 1 3
Gurpreet S. Wander 1 3
Kiymet Bozaoglu 1 3
Dharambir K. Sanghera 1 3 4
Paul Elliott 1 2 3
James Scott 0 1 3
John C. Chambers 1 2 3
Jaspal S. Kooner 0 1 3
0 National Heart and Lung Institute, Imperial College London, Hammersmith Hospital , Du Cane Road, London , United Kingdom , 3 Ealing Hospital NHS Trust , Southall, Middlesex , United Kingdom , 4 Genomics Computer Centre, South Texas Diabetes and Obesity Institute, University of Texas at the Rio Grande Valley, Brownsville, Texas, United States of America, 5 Broad Institute of Massachusetts Institute of Technology and Harvard, Massachusetts General Hospital , Cambridge, MA , United States of America, 6 Hero DMC Heart Institute , Ludhiana, Punjab , India , 7 Genomics and Systems Biology, Baker IDI Heart and Diabetes Institute , Melbourne , VIC Australia, 8 Department of Pediatrics, Section of Genetics, College of Medicine, University of Oklahoma Health Sciences Center , Oklahoma City, OK , United States of America
1 Funding: The LOLIPOP study is supported by the National Institute for Health Research Comprehensive Biomedical Research Centre Imperial College Healthcare NHS Trust , the British Heart Foundation (SP/04/002, https://
2 Epidemiology and Biostatistics, Imperial College London , Norfolk Place, London , United Kingdom
3 Editor: Nicholette D. Palmer, Wake Forest School of Medicine , UNITED STATES
4 Department of Pharmaceutical Sciences, College of Pharmacy, University of Oklahoma Health Sciences Center, Oklahoma City, OK, United States of America, 10 Oklahoma Center for Neuroscience, University of Oklahoma Health Sciences Center, Oklahoma City, OK, United States of America, 11 MRC-PHE Centre for Environment and Health, Imperial College London , Norfolk Place, London , United Kingdom , 12 Imperial College Healthcare NHS Trust , Du Cane Road, London , United Kingdom
South Asians are 1/4 of the world's population and have increased susceptibility to central obesity and related cardiometabolic disease. Knowledge of genetic variants affecting risk of central obesity is largely based on genome-wide association studies of common SNPs in Europeans. To evaluate the contribution of DNA sequence variation to the higher levels of central obesity (defined as waist hip ratio adjusted for body mass index, WHR) among South Asians compared to Europeans we carried out: i) a genome-wide association analysis of >6M genetic variants in 10,318 South Asians with focused analysis of population-specific SNPs; ii) an exome-wide association analysis of ~250K SNPs in protein-coding regions in 2,637 South Asians; iii) a comparison of risk allele frequencies and effect sizes of 48 known WHR SNPs in 12,240 South Asians compared to Europeans. In genome-wide analyses, we found no novel associations between common genetic variants and WHR in South Asians at P<5x10-8; variants showing equivocal association with WHR (P<1x10-5) did not replicate at P<0.05 in an independent cohort of South Asians (N = 1,922) or in published, predominantly European meta-analysis data. In the targeted analyses of 122,391 popula-
Data Availability Statement: All relevant data are
within the paper and its Supporting Information files.
tion-specific SNPs we also found no associations with WHR in South Asians at P<0.05 after
multiple testing correction. Exome-wide analyses showed no new associations between
genetic variants and WHR in South Asians, either individually at P<1.5x10-6 or grouped by
gene locus at P<2.5x10−6. At known WHR loci, risk allele frequencies were not higher in
South Asians compared to Europeans (P = 0.77), while effect sizes were unexpectedly
smaller in South Asians than Europeans (P<5.0x10-8). Our findings argue against an
collection and analysis, decision to publish, or
preparation of the manuscript.
Competing Interests: The authors have declared
that no competing interests exist.
important contribution for population-specific or cosmopolitan genetic variants underlying
the increased risk of central obesity in South Asians compared to Europeans.
Central obesity is a major risk factor for insulin resistance, type-2 diabetes (T2D), and
cardiovascular disease (CVD) [
]. South Asians, who comprise 1/4 of the world’s population,
are at high risk of developing central obesity compared to Europeans [
], and this is
widely considered to underlie their 2–4 fold excess risk of T2D and CVD [
Central obesity is heritable amongst South Asians [
], consistent with an important genetic
contribution. Genome-wide association studies (GWAS) predominantly carried out in people
of European ancestry have identified genetic loci underlying central obesity
]. Few studies have investigated genetic factors underlying central obesity
amongst South Asians [
To investigate whether genetic variation accounts for the increased risk of central obesity
amongst South Asians compared to Europeans, we carried out South Asian-specific
genomeand exome-wide association analyses and related our results to published European findings.
We selected waist hip ratio adjusted for body mass index (WHR) as a measure of central
obesity. We focused on population-specific SNPs present amongst South Asians but not
Europeans to test whether these influence WHR amongst South Asians. We also targeted
cosmopolitan SNPs associated with WHR in previous GWAS, and SNPs in protein coding
regions which are considered to have a higher probability of functional relevance.
Materials and Methods
We investigated the contribution of population-specific and cosmopolitan DNA sequence
variation to the increased risk of central obesity in South Asians using genome- and exome-wide
association analyses. The study design is described in Fig 1.
Population and phenotype
Samples for the discovery analyses were selected from the London Life Sciences Population
(LOLIPOP) study, an ongoing population-based cohort of 17,606 South Asian and 7,766
Northern European men and women, aged 35–75 years, recruited from the lists of 58 general
practitioners in West London. South Asians were recruited to the study if all 4 grandparents
were born in the Indian Subcontinent (countries of India, Pakistan, Sri Lanka or Bangladesh).
Data on medical history, current prescribed medication, and cardiovascular risk factors were
obtained by a trained research nurse using an interviewer-administered questionnaire. Country
of birth of participants, parents, and grandparents were recorded together with language and
religion, for assignment of ethnic subgroups. Physical measurements included blood pressure
(mean of 3 readings, taken with an Omron 705CP), height, weight, waist and hip
circumference. Blood was collected after an 8 hour fast for plasma glucose, total and HDL cholesterol,
triglycerides, insulin and high sensitivity C-reactive protein. The research was approved by the
West London Research Ethics Committee (reference number: 07/H0712/150), and all
participants gave written informed consent. Replication testing was done amongst 1,922 South Asian
participants from the Sikh Diabetes and Mauritius Family Studies (S1 Appendix; S1 Table)
2 / 17
Fig 1. Study Design.
Waist hip ratio adjusted for body mass index (WHR) was used as a measure of central
obesity. This allows assessment of abdominal fat independent of differences in overall adiposity
Sample selection. We selected 10,318 South Asians for genome-wide association, and
2,637 South Asians for exome-wide association (2,096 individuals were present in both
groups). We also carried out genome-wide association in 2,148 Europeans to enable
comparison between ethnic groups. Sample selection criteria and characteristics of participants are
summarised in S1 Table.
Genome-wide association. Genotyping platforms, calling algorithms, and quality control
measures used for genome-wide association are summarised in S2 Table. Hidden relatedness
or duplicate samples were sought using identity-by-descent methods implemented in PLINK;
samples with evidence for relatedness (pi_hat 0.5) were identified and one sample was
retained. Principal components analysis was carried out in Eigensoft v3.0; eigenvalues
inconsistent with either South Asian or European ancestry were removed.
Imputation of unmeasured genotypes was carried out amongst South Asians using
sequencing data from 321 South Asian participants in the LOLIPOP study [
], enabling the
3 / 17
investigation of population-specific genetic variants. In Europeans, imputation of unmeasured
genotypes was carried out using the cosmopolitan 1000 genomes reference panel (phase 1,
version 3). All imputation was performed using IMPUTE2; markers with low imputation info
score (<0.4), Hardy-Weinberg equilibrium P<1.0x10−6, or minor allele frequency (MAF)
<2% were excluded. We assessed the accuracy of our South Asian-specific imputation by
comparing imputed genotypes with corresponding direct genotyping results from our exome-wide
dataset among overlapping individuals (6,480 variants with MAF >2%). Compared to the
exome-wide dataset, mean sensitivity and specificity of the imputed SNPs were 99.41% (0.03)
and 99.40% (0.03) respectively. Mean concordance between direct and imputed genotypes was
98.31% (0.04). 6,571,328 SNPs were generated for South Asian-specific genome-wide
Exome-wide association. We genotyped 247,870 variants in protein-coding regions
among 2,637 South Asians, using the Illumina HumanExome Beadchip and
OmniExpressExome (S2 Table). Variants were called using Gencall, followed by post-processing with zCall to
improve calling of rare variants. Samples with a SNP call rate <95% were removed, as were
SNPs with call rate <97%, MAF <1%, or Hardy-Weinberg equilibrium P<1.0x10−6. 245,496
SNPs passed quality control, of which 35,368 SNPs had MAF >1%.
South Asian-specific GWAS. Single variants were examined for association with WHR
using linear regression using SNPTEST; covariate adjustments were made for age, sex, BMI,
and 10 principal components to control for residual population stratification. SNP-trait
associations were examined under additive, dominant and recessive inheritance models. Results were
normalised using an inverse normal transformed ranked scale, to enable comparison with
reported SNPs shown to be associated with WHR in previous GWAS. Each genotyping
platform was analysed separately, and results were combined by fixed effects inverse variance
meta-analysis within METAL. Heterogeneity was evaluated using the Cochran's Q statistic. A
significance threshold of P<5x10-8 was used for the discovery GWAS.
We refined our approach by focusing on a smaller number of genetic variants from our
genome-wide analysis which were population-specific (present in South Asians but not present
in Europeans within the 1000 genomes reference panel phase 1, version 3). Amongst the >6
million SNPs investigated in our genome-wide analysis we identified 122,391 South
Asian-specific SNPs with MAF >2% which we examined for association with WHR. For the analyses of
population-specific genetic variants, statistical significance was inferred at P<0.05 after
Bonferroni correction for the number of SNPs tested, thus enhancing study power compared to
South Asian-specific exome-wide association. Exonic variants were examined for associ
ation with WHR using linear regression and an additive genetic model within
RAREMETALWORKER, with adjustments for age, sex, and 10 principal components of ancestry. Results for
each genotyping platform were analysed separately and combined by meta-analysis within
METAL. Common, low-frequency, and rare variants were also grouped by genetic locus for
gene-based association analysis within RAREMETAL (CMC, Madsen-Browning, Variable
Threshold, and SKAT methods). We selected significance thresholds of P<1.5x10-6 and
P<2.5x10−6 for the single variant and gene-based analyses respectively (Bonferroni corrected).
Replication analysis. We carried out replication testing of all genetic variants showing
suggestive association with WHR amongst South Asians in genome-wide (P<1x10-5) and
exome-wide (P<1x10-3) analyses. Genome-wide results were tested for replication in an
independent cohort of South Asians (N = 1,922; S1 and S2 Tables; S1 Appendix) [
], and in
4 / 17
published, predominantly European meta-analysis data from the Giant Consortium (2012–
2014 release) [
]. For our exome-wide results we performed replication testing among
independent samples from our discovery South Asian GWAS analysis (N = 8,222). A replication
significance threshold of P<0.05 was selected.
Population comparisons. We carried out a literature search to identify previously
reported SNPs associated with central obesity in GWAS at P<5x10-8. We found 38 SNPs
associated with WHR in men and women [
], and 11 SNPs associated in women alone [
We selected 48 of the 49 known WHR SNPs for comparison between South Asians and
Europeans. One SNP, rs7801581 in HOXA11, was not available within our South Asian sample.
We analysed the 48 known WHR SNPs for association with WHR in South Asians at P<0.05.
We examined for systematic differences in risk allele frequencies and in effect sizes between
South Asian participants in LOLIPOP and reported European results using the sign test;
assessing whether the observed South Asian results, categorised individually as higher or lower than
the European reference, differed significantly from the binomial probability distribution
(expected value of 0.5). The contribution of known genetic variants to the excess in WHR
among South Asians compared to Europeans in the LOLIPOP study was quantified by
multivariate linear regression. To extend our analysis, we examined the 48 WHR SNPs in an
independent cohort of South Asians (N = 1,922), and used meta-analysis to combine the results
with our findings among South Asian participants in LOLIPOP.
Our South Asian GWAS has 80% power to identify single variants explaining 0.40% of the
variation in WHR among 10,318 South Asians at P<5x10-8 (S3A Table); power to detect variants
with effects equivalent to the average effect size of reported WHR variants is <1% at
genomewide significance. In our exome-array analysis we have 80% power to detect variants explaining
1.0% of the variation in WHR among 2,637 South Asians at P<1.5x10-6 (S3A Table). For our
population comparison of previously reported WHR SNPs from GWAS, we estimate between
10 and 84% power to replicate individual variants at P<0.05 among 10,318 South Asians (S3B
Table; or 15 and 91% among 12,240 South Asians, combined discovery and replication
To enhance study power we carried out targeted analyses of South Asian-specific SNPs. We
simulated the numbers of population-specific SNPs that would be required to explain the
increase in WHR amongst South Asians compared to Europeans as a function of allele
frequency and effect size (Fig 2). Amongst 10,318 South Asians we have >80% power to detect
common, population-specific variants associated with WHR with an effect size of >0.005 per
allele copy (untransformed β value). At this effect size there would need to be in excess of >20
common, population-specific common genetic variants associated with WHR amongst South
Asians (Fig 2).
The clinical characteristics of the 25,372 South Asian and European participants of the
LOLIPOP Study are shown in Table 1. Compared to Europeans, South Asians have higher waist hip
ratio (0.94 (0.08) v. 0.91 (0.08); P<0.001) despite similar body mass index (BMI) (27.4 (4.6) v.
27.5 (5.2); P = 0.23). South Asians also have a higher prevalence of T2D and CVD, higher
fasting glucose, insulin, HOMA-IR, HbA1c, and triglycerides, with lower HDL-cholesterol
compared to Europeans. Waist hip ratio remained higher amongst South Asians compared to
Europeans after adjustment for differences in age, gender, and BMI (Table 2).
5 / 17
Fig 2. Contribution of population-specific genetic variants to increased WHR amongst South Asians compared to Europeans. Results are shown as
the number of SNPs needed to fully explain the difference in WHR between the populations across a range of effect sizes (per allele copy, dashed black
lines). Superimposed are lines showing the power (10% power, orange line; 80% power, red line) of the current study to detect common and low frequency
genetic variants associated with WHR amongst South Asians. Results show that our study sample size is sufficient to identify common variants with effect
size >0.005 per allele copy, and rare / infrequent variants with effect size >0.025 (untransformed β values). At these effect sizes, there would need to be 10s
to 100s of population-specific genetic variants to explain increased WHR amongst South Asians. At effect sizes smaller than those identifiable in the current
study, there would need to be 100s of common variants or 1000s of rare / low frequency variants that are population-specific and associated with WHR, to
account for increased WHR amongst South Asians.
South Asian-specific GWAS
Genome-wide analysis. There were no genetic variants (MAF >2%) associated with
WHR at genome-wide significance under additive, dominant, or recessive inheritance models
(S4 Table and S1 Fig). Comparisons of the distributions for observed and expected association
results revealed little evidence for test statistic inflation (Lambdas additive = 1.025,
dominant = 1.015, recessive = 1.015, S2 Fig). Variants showing equivocal association with WHR
under an additive model (5x10-8<P<1x10-5) did not replicate at P<0.05, either in an
independent cohort of South Asians (N = 1,922) or in published, predominantly European data for the
Giant Consortium (S4 Table).
We examined population-specific genetic variants for
association with WHR in South Asians, adopting 4 MAF thresholds: i. >2% (N = 122,391); ii. >5%
(N = 38,639); iii. >10% (N = 7,349); and iv. >20% (N = 596). None of these population-specific
genetic variants were found to be associated with WHR in South Asians at P<0.05 after
Bonferroni correction for the number of SNPs tested in their respective allele frequency window
(S5 Table). In addition, the observed distributions of association statistics did not deviate from
the expected null distributions arguing against enrichment for association with WHR amongst
the population-specific SNPs tested (Lambda = 1.022, S3 Fig).
6 / 17
Coronary heart disease (%)
Type-2 diabetes (%)
Treated type-2 diabetes (%)
Treated hypertension (%)
Treated dyslipidaemia (%)
Body mass index (kg/m2)
Waist hip ratio
Systolic BP (mmHg)
Diastolic BP (mmHg)
LDL cholesterol (mmol/L)
HDL cholesterol (mmol/L)
Results are presented as mean (standard deviation), or % for categorical variables; comparisons using student t- and chi-squared tests.
South Asian-specific exome-wide association
None of the examined variants (MAF >1%) were found to be associated with WHR in South
Asians at P<1.5x10-6 (S6 Table and S4 Fig), and we found no evidence of deviation of the
observed distribution of the test statistics from that expected under the null hypothesis
(Lambda = 0.996, S5 Fig). In replication testing of equivocal variants (1.5x10-6<P<1x10-3),
only the rs17778003 SNP in ZFAT showed evidence of association with WHR at P<0.05, in
South Asians (P = 7.2x10-3) but not in Europeans (P = 0.58) (S6 Table). We found no genes
associated with WHR at Bonferroni genome-wide significance threshold of P<2.5x10−6 using
gene-based association analyses (S7 Table).
37 WHR SNPs–previously reported associations with WHR in men and women (P<5x10-8); 11 additional
WHR SNPs–previously reported associations with WHR in women alone (P<5x10-8).
Comparison of known genetic variants
We examined 37 known WHR SNPs for association with variation in WHR amongst South
Asian men and women participating in the LOLIPOP study. Only 4 achieved nominal
significance (rs6556301 near FGFR4, P = 6.3x10-03; rs984222 in TBX15, P = 6.6x10-3; rs1011731 in
DNM3, P = 1.8x10-2; rs12608504 near JUND, P = 2.7x10-02; Table 3). For the 11 SNPs known
to be associated with WHR in women alone, 2 had a nominally significant effect in South
Asian women (rs4684854 near PPARG, P = 1.3x10-2; rs1534696 in SNX10, P = 2.7x10-3;
We compared overall effect allele frequencies and effect sizes for these 48 SNPs between
South Asian participants in the LOLIPOP Study and reported European results. Effect allele
frequencies were similar between the two groups (sign test P = 0.67, Fig 3A). Effect sizes for WHR
were consistently lower in South Asians compared to Europeans (sign test P<1.7x10-6, Fig 3B).
Inclusion of results from a further 1,922 independent South Asians did not increase the number
of SNPs achieving nominal significance (S8 Table), and did not alter comparisons of effect allele
frequencies or effect sizes between populations (sign test P = 0.77 and P<5.0x10-8 respectively,
S6A and S6B Fig).
Together the 48 established WHR SNPs explained 0.7% of variation in WHR, equivalent to
2.5% of the reported ~27% heritability [
], among South Asians (LOLIPOP Study
participants, men and women combined). When included in multivariate regression analyses
examining the differences in WHR between the two populations, these variants did not account for
any of the excess risk of central obesity observed in South Asians compared to Europeans
Discussion and Conclusion
Central obesity, characterised by accumulation of excess body fat in an abdominal distribution,
is a leading risk factor for T2D, CVD, and premature mortality [
]. Central obesity
is more prevalent amongst South Asians compared to Europeans [
], but the
mechanisms underlying this are not well understood. Family studies show that central obesity is
heritable in South Asians [
]. Current knowledge of genetic loci influencing central obesity
risk is however largely based on the study of common variants using GWAS among Europeans
We carried out a GWAS of WHR in South Asians, using imputation of population-specific
and cosmopolitan genetic variants identified through next-generation sequencing in 321 South
]. We found that none of the >6 million evaluated variants were associated with
variation in WHR among 10,318 South Asians at genome-wide significance. We then carried out
a targeted analysis of population-specific SNPs to test whether these might account for
increased WHR amongst South Asians. Despite focus on a smaller number of genetic variants,
enhancing study power compared to genome-wide association, we again found no genetic
variants associated with WHR in South Asians. In parallel we carried out an exome-array analysis
among 2,637 South Asians to provide more comprehensive characterisation of genetic
variation in protein coding regions, which are considered to have a higher probability of functional
relevance. We found a single variant, rs17778003 in ZFAT, showing weak association with
WHR in South Asians, but no novel associations were identified at P<1.5x10-6. Further we
found no evidence of genes showing enrichment for common, low-frequency or rare
proteincoding variants underlying WHR in South Asians.
We modelled the number of population-specific SNPs required to explain the increase in
WHR amongst South Asians compared to Europeans as a function of allele frequency and
effect size, as well as study power to detect these variations. Our simulations show we are well
8 / 17
n D ,2853 ,034 ,3269 ,3663 ,4942 ,0281 ,8102 ,3057 ,9812 ,3532 ,119 ,28009 ,2275 ,8044 ,07890 ,41990 ,95770 ,18470 ,96189 ,28630 ,62380 ,43590 ,44770 ,21559 ,90690 ,86770 ,90990 ,25580 ,69160 ,09418 ,91290 ,34780 ,91490 ,05280 ,84787 ,91490 ,95290 ,02789 iteunnd
TE 11 95 10 11 10 9 9 10 10 11 93 1 8 8 2 2 2 2 1 2 2 2 2 1 2 2 2 2 2 1 2 2 2 2 1 2 2 1 o
P 50 5 05 5 05 5 5 5 5 4 5 05 05 5
n ILPOO ,3018 ,3018 ,3018 ,3018 ,3018 ,3018 ,3018 ,3018 ,3018 ,3018 ,3018 ,3018 ,0318 ,0381 ,0050 ,0050 ,504 ,505 505 ,1005 9504 ,9505 ,01005 AN ,010 ,950 ,100 ,950 ,100 ,950 ,1000 ,1000 ,1000 ,1000 ,950 ,010 ,100 ,950
L 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 9 9 9 0 ,
R T 9 1 5 A O
f48on trseeaG 3SPRO FVEAG -51TXBW L23FEN 41RBG LLLYPA -3PDNM I-2TSPR L86Y 13XCHO -3FZKRN I-SSCHN 4EPBC TSAADM 92DCCC 2PBM 1LEKR EPABC TEPM SCM L2TBN 2EYA 6SADM 11XAHO LLARCC 2TSDC JDNU 1L3FK I1ESM ARBOG LXPDN 7FXR 1AABC 5FDG 4FFRG TPSAA 13FAM ARCM
2 6 7 7 1 2 38 14 16 6 3
n ,74 ,60 ,566 ,895 ,435 ,81 ,66 ,56 ,50 ,226 ,641 R
D 69 73 1 1 1 96 11 11 10 1 1
1 1 1 1 1 H
0 0 0 0 0 0 0 0 0
powered to identify variants with modest phenotypic effects in our sample of 10,318 South
Asians (Fig 2). At these effect sizes there would need to be upwards of 20 population-specific
genetic variants associated with WHR amongst South Asians. In contrast we find none. We
Fig 3. Known WHR SNPs. (A) Comparison of risk allele frequencies in South Asians (N = 10,318) and Europeans
(reported). (B) Comparison of effect sizes (β WHR) in South Asians and Europeans. Green–men and women
combined (37 SNPs in South Asians (N = 10,318) and Europeans (reported)); orange–women alone (11 SNPs in South
Asians (N = 1,463) and Europeans (reported)).
11 / 17
therefore conclude that a small number of common, population-specific genetic variants with
modest effect size do not account for increased WHR amongst South Asians. Indeed, if
increased WHR amongst South Asians is the result population-specific genetic variants, this is
a polygenic disorder comprising very many genetic variants (>20) that are uncommon (MAF
<2%) and / or have small effect size (<0.005 increase in WHR per allele copy).
Differences in risk allele frequencies and effect sizes at known adiposity loci, including FTO
and MC4R, have been observed between South Asians and Europeans [
]. We compared
reported SNPs at 48 established WHR loci in South Asians and in Europeans to investigate
whether these variants underlie the excess risk of central obesity observed in South Asians.
When compared with published European results, we observed limited evidence for replication
at known WHR SNPs among South Asians; only 34 of 48 WHR loci showed directionally
consistent effects on WHR in South Asians. When all 48 known WHR variants were examined
together we observed no systematic differences in risk allele frequencies, and consistently
smaller effect sizes, in South Asians compared to Europeans. Combined these variants
accounted for <1% of phenotypic variation in WHR among South Asians, and did not
contribute to the excess of WHR observed in South Asians compared to Europeans. The reasons for
the lack of replication of known WHR SNPs in South Asians are not known but may include:
(i) genetic loci with European-specific effects; (ii) stronger relationships between tag and causal
SNPs in Europeans because of differences in haplotype structure; (iii) inflated effect sizes in the
European discovery sample due to winner's curse; (iv) greater phenotypic heterogeneity in
South Asians, who have smaller stature but greater central adiposity than Europeans, with
reduced power to detect SNP-trait associations. Nevertheless, our findings robustly
demonstrate that known WHR SNPs do not account for the increased risk of central obesity amongst
South Asians compared to Europeans.
Mechanisms underlying central obesity risk among South Asians remain unclear. We have
not excluded the possibility of a genetic contribution through multiple common variants with
small effects, rare variants with larger phenotypic effects, or structural variations that are
inadequately captured using current analytic platforms. Environmental influences have an
important role in central adiposity [
], and may contribute to increased risk of
central obesity in South Asians. For example differences in lifestyle, such as lower levels of
physical activity [
], and higher levels of total calorie, refined starch and saturated fat
], have been reported in South Asians compared to Europeans. Similarly,
the prevalence of low birth weight, which is associated with accelerated childhood weight gain
] and future development of central obesity and metabolic disturbance
], is higher in South Asians than Europeans [
explanation is that population-specific epigenomic modifications [
], which regulate
gene expression and phenotypic variation without change in DNA sequence [
to central obesity predisposition amongst South Asians. These modifications in DNA
methylation, histone modification, and chromatin remodelling, can be transmitted through the
germline and modified by environmental exposures [
], providing a compelling
putative mechanism for unexplained phenotypic variation. Future initiatives combining new,
more targeted analysis strategies and larger sample sizes will be required to elucidate the
relative roles of these genomic and environmental factors in the excess of central obesity among
S1 Appendix. South Asian replication cohort description.
12 / 17
S1 Fig. South Asian GWAS–Association of SNPs with WHR by chromosome position. (A)
Additive inheritance model. (B) Dominant inheritance model. (C) Recessive inheritance
S2 Fig. South Asian GWAS–Relationship of observed association statistics with those
expected under the null distribution. (A) Additive inheritance model (Lambda = 1.025). (B)
Dominant inheritance model (Lambda = 1.015). (C) Recessive inheritance model
(Lambda = 1.015).
S3 Fig. South-Asian GWAS targeted analysis of population-specific SNPs. Relationship of
observed association statistics with those expected under the null distribution (MAF>2%,
Lambda = 1.022).
S4 Fig. South Asian exome-array analysis–Association of SNPs with WHR by chromosome
S5 Fig. South Asian exome-array analysis–Relationship of observed association statistics
with those expected under the null distribution (Lambda = 0.996).
S6 Fig. Known WHR SNPs. (A) Comparison of risk allele frequencies in South Asians,
extended analysis (N = 12,240) and Europeans (reported). (B) Comparison of effect sizes (β
WHR) in South Asians and Europeans, extended analysis. Green–men and women combined
(37 SNPs in South Asians (N = 12,240) and Europeans (reported)); orange–women alone (11
SNPs in South Asians (N = 2,363) and Europeans (reported)).
S1 Table. Characteristics of participants in the genotyping cohorts.
S2 Table. Genotyping cohort platforms and quality metrics.
S3 Table. Study power to identify single variants associated with variation in WHR
amongst South Asians, and replicate known WHR SNPs amongst South Asians.
S4 Table. South Asian GWAS–Top ranking single markers at P<1x10-5 under an additive,
dominant and recessive inheritance models inheritance model.
S5 Table. South Asian GWAS targeted analysis–Top ranking population-specific markers
under an additive inheritance model.
S6 Table. South Asian exome-array analysis–Top ranking single markers at P<1x10-3, with
replication among South Asians (GWAS) and Europeans (GIANT consortium
13 / 17
S7 Table. South Asian exome-array analysis–Top ranking gene-based association results at
analysis) and Europeans.
S8 Table. Comparison of 48 known WHR SNPs in South Asians (N = 12,240; extended
The LOLIPOP study is supported by the National Institute for Health Research (NIHR)
Comprehensive Biomedical Research Centre Imperial College Healthcare NHS Trust, the British
Heart Foundation (SP/04/002), the Medical Research Council (G0601966,G0700931), the
Wellcome Trust (084723/Z/08/Z), and the NIHR (RP-PG-0407-10371). The work was carried
out in part at the NIHR/Wellcome Trust Imperial Clinical Research Facility. The Sikh Diabetes
Study is supported by National Institute of Health grants KO1TW006087, funded by the
Fogarty International Center, R01DK082766, funded by National Institute of Diabetes and
Digestive and Kidney Diseases, and a seed grant from University of Oklahoma Health Sciences
Center, Oklahoma City, USA. The Mauritius Family Study is supported by the Mauritius
Ministry of Health and Quality of Life, Australian Government National Health and Medical
Research Council NHMRC project grant numbers 1020285 and 1037916, the Victorian
Government’s OIS Program, and partly funded by US National Institutes of Health Grant
DK25446. We thank the participants and research staff who made the study possible.
Conceived and designed the experiments: WRS JCC JSK. Performed the experiments: WRS
WZ ML ST UA JP RS SR GSW JCC JSK. Analyzed the data: WRS WZ ML BL JP RS SR GSW
KB. Contributed reagents/materials/analysis tools: KB DKS PE JS JCC JSK. Wrote the paper:
WRS WZ ST KB DKS PE JS JCC JSK.
14 / 17
15 / 17
16 / 17
1. Despres JP , Lemieux I ( 2006 ) Abdominal obesity and metabolic syndrome . Nature 444 : 881 - 887 . PMID: 17167477
2. Langenberg C , Sharp SJ , Schulze MB , Rolandsson O , Overvad K , Forouhi NG , et al. (2012) Long-term risk of incident type 2 diabetes and measures of overall and regional obesity: the EPIC-InterAct casecohort study . PLoS medicine 9: e1001230. doi: 10.1371/journal.pmed.1001230 PMID: 22679397
3. de Koning L , Merchant AT , Pogue J , Anand SS ( 2007 ) Waist circumference and waist-to-hip ratio as predictors of cardiovascular events: meta-regression analysis of prospective studies . European heart journal 28 : 850 - 856 . PMID: 17403720
4. Yusuf S , Hawken S , Ounpuu S , Bautista L , Franzosi MG , Commerford P , et al. ( 2005 ) Obesity and the risk of myocardial infarction in 27,000 participants from 52 countries: a case-control study . Lancet 366 : 1640 - 1649 . PMID: 16271645
5. Misra A , Khurana L ( 2011 ) Obesity-related non-communicable diseases: South Asians vs White Caucasians . International Journal of Obesity 35 : 167 - 187 . doi: 10 .1038/ijo. 2010 .135 PMID: 20644557
6. Kohli S , Sniderman AD , Tchernof A , Lear SA ( 2010 ) Ethnic-specific differences in abdominal subcutaneous adipose tissue compartments . Obesity 18 : 2177 - 2183 . doi: 10 .1038/oby. 2010 .94 PMID: 20448537
7. Tillin T , Forouhi NG , McKeigue PM , Chaturvedi N ( 2012 ) Southall And Brent REvisited: Cohort profile of SABRE, a UK population-based comparison of cardiovascular disease and diabetes in people of European, Indian Asian and African Caribbean origins . International journal of epidemiology 41 : 33 - 42 . doi: 10 .1093/ije/dyq175 PMID: 21044979
8. Cameron AJ , Sicree RA , Zimmet PZ , Alberti KG , Tonkin AM , Balkau B , et al. ( 2010 ) Cut-points for waist circumference in Europids and South Asians . Obesity 18 : 2039 - 2046 . doi: 10 .1038/oby. 2009 .455 PMID: 20019679
9. McKeigue PM , Pierpoint T , Ferrie JE , Marmot MG ( 1992 ) Relationship of glucose intolerance and hyperinsulinaemia to body fat pattern in south Asians and Europeans . Diabetologia 35 : 785 - 791 . PMID: 1511807
10. McKeigue PM , Shah B , Marmot MG ( 1991 ) Relation of central obesity and insulin resistance with high diabetes prevalence and cardiovascular risk in South Asians . Lancet 337 : 382 - 386 . PMID: 1671422
11. Zaninotto P , Mindell J , Hirani V ( 2007 ) Prevalence of cardiovascular risk factors among ethnic groups: results from the Health Surveys for England . Atherosclerosis 195 : e48 - 57 . PMID: 17395185
12. Wild SH , Fischbacher C , Brock A , Griffiths C , Bhopal R ( 2007 ) Mortality from all causes and circulatory disease by country of birth in England and Wales 2001-2003 . Journal of public health 29 : 191 - 198 . PMID: 17456532
13. Barnett AH , Dixon AN , Bellary S , Hanif MW , O'Hare JP , Raymond NT , et al. ( 2006 ) Type 2 diabetes and cardiovascular risk in the UK south Asian community . Diabetologia 49 : 2234 - 2246 . PMID: 16847701
14. HSE ( 2004 ) Health Survery for England: Health of ethnic minorities .
15. Davey G , Ramachandran A , Snehalatha C , Hitman GA , McKeigue PM ( 2000 ) Familial aggregation of central obesity in Southern Indians . International journal of obesity and related metabolic disorders: journal of the International Association for the Study of Obesity 24 : 1523 - 1527 .
16. Zabaneh D , Chambers JC , Elliott P , Scott J , Balding DJ , Kooner JS ( 2009 ) Heritability and genetic correlations of insulin resistance and component phenotypes in Asian Indian families using a multivariate analysis . Diabetologia 52 : 2585 - 2589 . doi: 10 .1007/s00125-009 -1504-7 PMID: 19763535
17. Lindgren CM , Heid IM , Randall JC , Lamina C , Steinthorsdottir V , Qi L , et al. ( 2009 ) Genome-wide association scan meta-analysis identifies three Loci influencing adiposity and fat distribution . PLoS genetics 5: e1000508. doi: 10.1371/journal.pgen.1000508 PMID: 19557161
18. Heard-Costa NL , Zillikens MC , Monda KL , Johansson A , Harris TB , Fu M , et al. ( 2009 ) NRXN3 is a novel locus for waist circumference: a genome-wide association study from the CHARGE Consortium . PLoS genetics 5: e1000539. doi: 10.1371/journal.pgen.1000539 PMID: 19557197
19. Heid IM , Jackson AU , Randall JC , Winkler TW , Qi L , Steinthorsdottir V , et al. ( 2010 ) Meta-analysis identifies 13 new loci associated with waist-hip ratio and reveals sexual dimorphism in the genetic basis of fat distribution . Nature genetics 42 : 949 - 960 . doi: 10 .1038/ng.685 PMID: 20935629
20. Randall JC , Winkler TW , Kutalik Z , Berndt SI , Jackson AU , Monda KL , et al. ( 2013 ) Sex-stratified genome-wide association studies including 270,000 individuals show sexual dimorphism in genetic loci for anthropometric traits . PLoS genetics 9: e1003500. doi: 10.1371/journal.pgen.1003500 PMID: 23754948
21. Fox CS , Liu Y , White CC , Feitosa M , Smith AV , Heard-Costa N , et al. ( 2012 ) Genome-wide association for abdominal subcutaneous and visceral adipose reveals a novel locus for visceral fat in women . PLoS genetics 8: e1002695. doi: 10.1371/journal.pgen.1002695 PMID: 22589738
22. Shungin D , Winkler TW , Croteau-Chonka DC , Ferreira T , Locke AE , Mägi R , et al. ( 2015 ) New genetic loci link adipose and insulin biology to body fat distribution . Nature 518 : 187 - 196 . doi: 10 .1038/ nature14132 PMID: 25673412
23. Chambers JC , Elliott P , Zabaneh D , Zhang W , Li Y , Froguel P , et al. ( 2008 ) Common genetic variation near MC4R is associated with waist circumference and insulin resistance . Nature genetics 40 : 716 - 718 . doi: 10 .1038/ng.156 PMID: 18454146
24. Zabaneh D , Balding DJ ( 2010 ) A genome-wide association study of the metabolic syndrome in Indian Asian men . PLoS One 5 : e11961 . doi: 10.1371/journal.pone.0011961 PMID: 20694148
25. Sanghera DK , Bhatti JS , Bhatti GK , Ralhan SK , Wander GS , Singh JR , et al. ( 2006 ) The Khatri Sikh Diabetes Study (SDS): study design, methodology, sample collection, and initial results . Human biology 78 : 43 - 63 . PMID: 16900881
26. Soderberg S , Zimmet P , Tuomilehto J , de Courten M , Dowse GK , Chitson P , et al. ( 2004 ) High incidence of type 2 diabetes and increasing conversion rates from impaired fasting glucose and impaired glucose tolerance to diabetes in Mauritius . Journal of internal medicine 256 : 37 - 47 . PMID: 15189364
27. Soderberg S , Zimmet P , Tuomilehto J , de Courten M , Dowse GK , Chitson P , et al. ( 2005 ) Increasing prevalence of Type 2 diabetes mellitus in all ethnic groups in Mauritius . Diabetic medicine: a journal of the British Diabetic Association 22 : 61 - 68 .
28. Jowett JB , Diego VP , Kotea N , Kowlessur S , Chitson P , Dyer TD , et al. ( 2009 ) Genetic influences on type 2 diabetes and metabolic syndrome related quantitative traits in Mauritius. Twin research and human genetics: the official journal of the International Society for Twin Studies 12 : 44 - 52 .
29. Chambers JC , Abbott J , Zhang W , Turro E , Scott WR , Tan ST , et al. ( 2014 ) The South Asian genome . PLoS One 9 : e102645 . doi: 10.1371/journal.pone.0102645 PMID: 25115870
30. Pischon T , Boeing H , Hoffmann K , Bergmann M , Schulze MB , Overvad K , et al. ( 2008 ) General and abdominal adiposity and risk of death in Europe . The New England journal of medicine 359 : 2105 - 2120 . doi: 10 .1056/NEJMoa0801891 PMID: 19005195
31. Boggs DA , Rosenberg L , Cozier YC , Wise LA , Coogan PF , Ruiz-Narvaez EA , et al. ( 2011 ) General and abdominal obesity and risk of death among black women . The New England journal of medicine 365 : 901 - 908 . doi: 10 .1056/NEJMoa1104119 PMID: 21899451
32. Li H , Kilpelainen TO , Liu C , Zhu J , Liu Y , Hu C , et al. ( 2012 ) Association of genetic variation in FTO with risk of obesity and type 2 diabetes with data from 96,551 East and South Asians . Diabetologia 55 : 981 - 995 . doi: 10 .1007/s00125-011 -2370-7 PMID: 22109280
33. Schleinitz D , Bottcher Y , Bluher M , Kovacs P ( 2014 ) The genetics of fat distribution . Diabetologia 57 : 1276 - 1286 . doi: 10 .1007/s00125-014-3214-z PMID: 24632736
34. Murakami K , McCaffrey TA , Livingstone MB ( 2013 ) Associations of dietary glycaemic index and glycaemic load with food and nutrient intake and general and central obesity in British adults . The British journal of nutrition 110 : 2047 - 2057 . doi: 10 .1017/S0007114513001414 PMID: 23656860
35. Bradlee ML , Singer MR , Qureshi MM , Moore LL ( 2010 ) Food group intake and central obesity among children and adolescents in the Third National Health and Nutrition Examination Survey (NHANES III) . Public health nutrition 13 : 797 - 805 . doi: 10 .1017/S1368980009991546 PMID: 19772691
36. Wang Y , Beydoun MA ( 2009 ) Meat consumption is associated with obesity and central obesity among US adults . International Journal of Obesity 33 : 621 - 628 . doi: 10 .1038/ijo. 2009 .45 PMID: 19308071
37. Brandhagen M , Forslund HB , Lissner L , Winkvist A , Lindroos AK , Carlsson LMS , et al. ( 2012 ) Alcohol and macronutrient intake patterns are related to general and central adiposity . European journal of clinical nutrition 66 : 305 - 313 . doi: 10 .1038/ejcn. 2011 .189 PMID: 22085868
38. Kondoh T , Takase H , Yamaguchi TF , Ochiai R , Katashima M , Katsuragi Y , et al. ( 2014 ) Association of dietary factors with abdominal subcutaneous and visceral adiposity in Japanese men . Obesity research & clinical practice 8 : e16 - 25 .
39. Hayes L , White M , Unwin N , Bhopal R , Fischbacher C , Harland J , et al. ( 2002 ) Patterns of physical activity and relationship with risk markers for cardiovascular disease and diabetes in Indian, Pakistani, Bangladeshi and European adults in a UK population . Journal of public health medicine 24 : 170 - 178 . PMID: 12831085
40. Fischbacher CM , Hunt S , Alexander L ( 2004 ) How physically active are South Asians in the United Kingdom? A literature review . Journal of public health 26 : 250 - 258 . PMID: 15454592
41. ( 1994 ) UK Prospective Diabetes Study. XII: Differences between Asian, Afro-Caribbean and white Caucasian type 2 diabetic patients at diagnosis of diabetes . UK Prospective Diabetes Study Group. Diabetic medicine: a journal of the British Diabetic Association 11 : 670 - 677 .
42. Dhawan J , Bray CL ( 1997 ) Asian Indians, coronary artery disease, and physical exercise . Heart 78 : 550 - 554 . PMID: 9470869
43. Misra A , Khurana L , Isharwal S , Bhardwaj S ( 2009 ) South Asian diets and insulin resistance . The British journal of nutrition 101 : 465 - 473 . doi: 10 .1017/S0007114508073649 PMID: 18842159
44. Merchant AT , Anand SS , Kelemen LE , Vuksan V , Jacobs R , Davis B , et al. ( 2007 ) Carbohydrate intake and HDL in a multiethnic population . The American Journal of Clinical Nutrition 85 : 225 - 230 . PMID: 17209200
45. Donin AS , Nightingale CM , Owen CG , Rudnicka AR , McNamara MC , Prynne CJ , et al. ( 2010 ) Nutritional composition of the diets of South Asian, black African-Caribbean and white European children in the United Kingdom: the Child Heart and Health Study in England (CHASE) . The British journal of nutrition 104 : 276 - 285 . doi: 10 .1017/S000711451000070X PMID: 20230652
46. Donin AS , Nightingale CM , Owen CG , Rudnicka AR , McNamara MC , Prynne CJ , et al. ( 2010 ) Ethnic differences in blood lipids and dietary intake between UK children of black African, black Caribbean, South Asian, and white European origin: the Child Heart and Health Study in England (CHASE) . The American Journal of Clinical Nutrition 92 : 776 - 783 . doi: 10 .3945/ajcn. 2010 .29533 PMID: 20739425
47. Meas T ( 2010 ) Fetal origins of insulin resistance and the metabolic syndrome: a key role for adipose tissue? Diabetes & metabolism 36 : 11 - 20 .
48. Amador-Licona N , Martinez-Cordero C , Guizar-Mendoza JM , Malacara JM , Hernandez J , Alcala JF ( 2007 ) Catch-up growth in infants born small for gestational age-a longitudinal study . Journal of pediatric endocrinology & metabolism: JPEM 20 : 379 - 386 .
49. Mericq V , Ong KK , Bazaes R , Pena V , Avila A , Salazar T , et al. ( 2005 ) Longitudinal changes in insulin sensitivity and secretion from birth to age three years in small- and appropriate-for-gestational-age children . Diabetologia 48 : 2609 - 2614 . PMID: 16283238
50. Oken E , Gillman MW ( 2003 ) Fetal origins of obesity . Obesity research 11 : 496 - 506 . PMID: 12690076
51. Rogers I ( 2003 ) The influence of birthweight and intrauterine environment on adiposity and fat distribution in later life . International journal of obesity and related metabolic disorders: journal of the International Association for the Study of Obesity 27 : 755 - 777 .
52. Lindsay RS , Dabelea D , Roumain J , Hanson RL , Bennett PH , Knowler WC ( 2000 ) Type 2 diabetes and low birth weight: the role of paternal inheritance in the association of low birth weight and diabetes . Diabetes 49 : 445 - 449 . PMID: 10868967
53. Martinson ML , McLanahan S , Brooks-Gunn J ( 2015 ) Variation in child body mass index patterns by race/ethnicity and maternal nativity status in the United States and England . Maternal and child health journal 19 : 373 - 380 . doi: 10 .1007/s10995-014 -1519-7 PMID: 24894727
54. Ramadhani MK , Grobbee DE , Bots ML , Castro Cabezas M , Vos LE , Oren A , et al. ( 2006 ) Lower birth weight predicts metabolic syndrome in young adults: the Atherosclerosis Risk in Young Adults (ARYA)- study . Atherosclerosis 184 : 21 - 27 . PMID: 16326169
55. Iniguez G , Ong K , Bazaes R , Avila A , Salazar T , Dunger D , et al. ( 2006 ) Longitudinal changes in insulin-like growth factor-I, insulin sensitivity, and secretion from birth to age three years in small-for-gestational-age children . The Journal of clinical endocrinology and metabolism 91 : 4645 - 4649 . PMID: 16912131
56. WHO ( 2004 ) World Health Organisation. Low birthweight.Country, regional and global estimates .
57. Wells JC , Sharp G , Steer PJ , Leon DA ( 2013 ) Paternal and maternal influences on differences in birth weight between Europeans and Indians born in the UK . PLoS One 8 : e61116 . doi: 10.1371/journal. pone.0061116 PMID: 23667432
58. Veena SR , Kumaran K , Swarnagowri MN , Jayakumar MN , Leary SD , Stein CE , et al. ( 2004 ) Intergenerational effects on size at birth in South India . Paediatric and perinatal epidemiology 18 : 361 - 370 . PMID: 15367323
59. Leon DA , Moser KA ( 2012 ) Low birth weight persists in South Asian babies born in England and Wales regardless of maternal country of birth. Slow pace of acculturation, physiological constraint or both? Analysis of routine data . Journal of epidemiology and community health 66 : 544 - 551 . doi: 10 .1136/ jech. 2010 .112516 PMID: 21118952
60. Veena SR , Geetha S , Leary SD , Saperia J , Fisher DJ , Kumaran K , et al. ( 2007 ) Relationships of maternal and paternal birthweights to features of the metabolic syndrome in adult offspring: an inter-generational study in South India . Diabetologia 50 : 43 - 54 . PMID: 17143606
61. Harding S , Rosato MG , Cruickshank JK ( 2004 ) Lack of change in birthweights of infants by generational status among Indian, Pakistani, Bangladeshi, Black Caribbean, and Black African mothers in a British cohort study . International journal of epidemiology 33 : 1279 - 1285 . PMID: 15155695
62. Fraser HB , Lam LL , Neumann SM , Kobor MS ( 2012 ) Population-specificity of human DNA methylation . Genome biology 13: R8 . doi: 10 .1186/gb-2012 -13-2-r8 PMID: 22322129
63. Moen EL , Zhang X , Mu W , Delaney SM , Wing C , McQuade J , et al. ( 2013 ) Genome-wide variation of cytosine modifications between European and African populations and the implications for complex traits . Genetics 194 : 987 - 996 . doi: 10 .1534/genetics.113.151381 PMID: 23792949
64. Barfield RT , Almli LM , Kilaru V , Smith AK , Mercer KB , Duncan R , et al. ( 2014 ) Accounting for population stratification in DNA methylation studies . Genetic epidemiology 38 : 231 - 241 . doi: 10 .1002/gepi.21789 PMID: 24478250
65. Heyn H , Moran S , Hernando-Herraez I , Sayols S , Gomez A , Sandoval J , et al. ( 2013 ) DNA methylation contributes to natural human variation . Genome research 23 : 1363 - 1372 . doi: 10 .1101/gr.154187.112 PMID: 23908385
66. Petronis A ( 2010 ) Epigenetics as a unifying principle in the aetiology of complex traits and diseases . Nature 465 : 721 - 727 . doi: 10 .1038/nature09230 PMID: 20535201
67. Feinberg AP ( 2007 ) Phenotypic plasticity and the epigenetics of human disease . Nature 447 : 433 - 440 . PMID: 17522677
68. Feil R , Fraga MF ( 2011 ) Epigenetics and the environment: emerging patterns and implications . Nature reviews Genetics 13 : 97 - 109 .
69. Johnstone SE , Baylin SB ( 2010 ) Stress and the epigenetic landscape: a link to the pathobiology of human diseases? Nature reviews Genetics 11 : 806 - 812 . doi: 10 .1038/nrg2881 PMID: 20921961
70. Jirtle RL , Skinner MK ( 2007 ) Environmental epigenomics and disease susceptibility . Nature reviews Genetics 8 : 253 - 262 . PMID: 17363974
71. Portela A , Esteller M ( 2010 ) Epigenetic modifications and human disease . Nature biotechnology 28 : 1057 - 1068 . doi: 10 .1038/nbt.1685 PMID: 20944598