Abstract

The prevalence of Type 2 diabetes is rising rapidly in both developed and developing countries. Asia is developing as the epicentre of the escalating pandemic, reflecting rapid transitions in demography, migration, diet, and lifestyle patterns. The effective management of Type 2 diabetes in Asia may be complicated by differences in prevalence, risk factor profiles, genetic risk allele frequencies, and gene-environment interactions between different Asian countries, and between Asian and other continental populations. To reduce the worldwide burden of T2D, it will be important to understand the architecture of T2D susceptibility both within and between populations. This review will provide an overview of known genetic and nongenetic risk factors for T2D, placing the results from Asian studies in the context of broader global research. Given recent evidence from large-scale genetic studies of T2D, we place special emphasis on emerging knowledge about the genetic architecture of T2D and the potential contribution of genetic effects to population differences in risk.

1. Introduction

Type 2 diabetes (T2D) is one of the top five noncommunicable diseases globally, comprising a major, growing cause of morbidity and premature death. In 2012, the International Diabetes Federation (IDF) estimated that 371 million people worldwide were living with diabetes, of which about half live in South Asia, the Western Pacific, and Eastern Mediterranean regions [1]. Asia is now the epicenter of an escalating diabetes epidemic, chiefly due to population growth and ageing in India and China. Projections suggest that by 2030, more than 60% of worldwide diabetes cases will come from Asia [2, 3], with the vast majority of these being Type 2 diabetes (T2D) [4]. T2D has an enormous economic, psychosocial, and physical impact on individuals, their families, and communities, both directly and indirectly. The direct economic burden of T2D includes both recorded expenditure by health services and unrecorded costs borne by patients and their families. Indirect costs such as loss of productivity and disability are also substantial and may match or surpass direct costs. The proportion of worldwide disability-adjusted life years (DALYs) due to T2D has soared in recent decades, rising from 43% in 1990 to 54% in 2010 [5]. Temporary and permanent disability, excess morbidity, and premature death are the consequences of T2D vascular complications, including cardiovascular disease, retinopathy (blindness), nephropathy (kidney failure), and neuropathy (nerve problems) which can lead to amputation. Intangible costs due to psychosocial effects on quality of life, diminished contribution to family tasks, and reduced income of care-giving family members are also likely substantial but difficult to assess. The enormous, growing global burden of T2D—particularly in Asia—is now viewed as a crisis by the World Health Organisation (WHO) and the United Nations (UN) [6]. There is a major worldwide push to decrease the prevalence and impact of T2D by identifying risk factors, both genetic and nongenetic. Explaining the distribution and variation of T2D susceptibility across Asia will be vital for reducing the global burden of disease, due to the demographic, cultural, and genetic heterogeneity of Asian populations, and T2D risk factor profiles between these populations [710].

2. Epidemiology

2.1. Burden of the Disease

The vast majority of T2D (about 80%) occurs in low- and middle-income countries (LMICs), with India and China providing the largest absolute contributions. The prevalence of T2D is also rising most swiftly in LMICs [6], particularly in Asian countries experiencing rapid economic growth (Figure 2). However, there are disparities in T2D prevalence among Asian populations; Asians from the Indian subcontinent (India, Pakistan and Bangladesh) have the highest prevalence (15.9% to 24.9%), with intermediate prevalence in Malays (11.4% to 16.9%) and reduced prevalence in Chinese (6.4% to 13.8%) [1113]. These risk profile differences may reflect population differences in T2D risk due to ethnicity-specific diet and lifestyle, body composition, genetic effects, or gene-environment interactions, as discussed further in the sections below.

2.2. Pathophysiology

The pathogenesis of Type 2 diabetes (T2D) involves deficient insulin secretion by pancreatic β-cells, and diminished insulin effectiveness in target tissues (insulin resistance) T2D aetiology differs from that of Type 1 diabetes (T1D), in which there is absolute insulin deficiency due to the destruction of insulin-producing β-cells [14]. T2D represents 90% of all diabetes cases worldwide [4]. Impaired insulin secretion and insulin action led to an accumulation of glucose in the blood (hyperglycaemia), with adverse effects on health. Clinical features of hyperglycaemia and T2D include excessive excretion of urine (polyuria), thirst (polydipsia), constant hunger, weight loss, vision change, and fatigue [15]. These symptoms may occur suddenly but are often less marked, and T2D patients may be unaware of their illness for several years until further complications develop.

2.2.1. Insulin Resistance

Glucose homeostasis depends upon a highly regulated feedback system comprising both insulin-secreting β-cells and insulin-sensitive target tissues. The function of either component—while accounting for the associated homeostatic response of the other—can be evaluated using Homeostasis Model Assessment (HOMA) [16]. Studies assessing insulin resistance using HOMA (HOMA-IR) report continental differences in the relative contribution of insulin deficiency versus insulin resistance to T2D. Compared to healthy European-ancestry participants matched for age and body mass index (BMI), Asian Indian individuals exhibit higher insulin resistance [17] and a greater contribution of insulin resistance—relative to insulin secretion—to T2D pathogenesis [18]. One study evaluating insulin response to a fixed glucose load also showed that Japanese-Americans displayed an insulin response more similar to native Japanese than European-Americans, in spite of sharing a highly Westernized lifestyle with their European-American counterparts [8].

There is also variation in the predisposition to insulin resistance between Asian populations [19]. For several decades, it has been recognised that the highest propensity is present in Asian Indians, in whom insulin resistance contributes substantially to T2D pathogenesis [20], potentially reflecting ancestry-related predisposition to abdominal obesity [21, 22]. A recent population based study of 4,136 Chinese, Malays, and Asian Indians living in Singapore supported these findings, reporting substantially higher insulin resistance in Asian Indians, intermediate levels in Malays, and the lowest levels in Chinese () [19]. Differences between Malays and Chinese were removed after adjusting for body mass index (BMI); the remaining additional resistance in Indians appeared to be mediated by a tendency to higher BMI and BMI-adjusted waist circumference, together with other unexplained factors [19].

Dickinson and colleagues studied postprandial hyperglycemia and insulin sensitivity after a 75 gram carbohydrate challenge in 60 lean, healthy individuals from five ethnic groups with similar age, BMI, waist circumference, and birth weight distributions. Prior to carbohydrate consumption, fasting insulin was significantly higher in South Indians and South East Asians, compared to European Caucasian, Arabic, and Chinese individuals () [23]. Following the challenge, hyperglycemia was significantly higher in South East Asian and Chinese participants compared with European Caucasians, while Indians and South East Asians showed a 2-3-fold higher insulin response than Europeans [23]. A small Singapore-based study of 30 individuals also showed significantly reduced insulin sensitivity in South Indians compared with Chinese or European individuals matched for age, BMI, and physical activity [24].

2.2.2. Insulin Secretion

Impaired insulin secretion is associated with β-cell dysfunction that results in a reduced insulin-secretion response to rises in blood glucose after eating [25]. The insulin secretion response to various foods can be quantified using the insulin index; more complex relationship between insulin secretion and insulin sensitivity can be measured using the disposition index (DI), which is assessed by an intravenous glucose tolerance-test [26]. A recent family-based study found that a high-fat, low-carbohydrate dietary pattern contributed to obesity, insulin resistance, and reduced β-cell function [27]. This finding might be explained by increased free fatty acids (FFAs) reducing the expression of β-cell—specific transcription factors and impairing the β-cells’ ability to respond to glucose with appropriate insulin secretion [28].

Similar to insulin resistance, insulin secretion also shows evidence of racial differences, being reduced in Asians compared with Europeans. The insulin index of Asians is reduced almost 70% in the progression from impaired glucose tolerance (IGT) to T2D, whereas in Europeans the corresponding reduction is only 50% [29, 30]. A population based-cohort study of insulin resistance and β-cell function during pregnancy also found a significantly lower β-cell secretory response to pregnancy-induced insulin resistance in South Asian and East Asian women, compared to European participants with a similar level of insulin resistance [31].

2.2.3. Complications

T2D complications can be life-threatening and include cardiovascular disease, nephropathy (kidney disease), retinopathy (blindness), and neuropathy (nerve impairment). Observational studies in European American and African American population report that cardiovascular disease risk in individuals with T2D is more than double the rate in the general population [32] and 50% of people with T2D die from cardiovascular disease, primarily heart disease and stroke [33].

There is evidence for population differences in the rate of T2D complications, between Asian populations and broader continental groups. A cross-sectional study of 5,707 Chinese, Indians, and Malays showed that the population attributable risk of ischaemic heart disease related to T2D was the highest in Indians (40.9%), intermediate in Malays (27.9%), and the lowest in Chinese (11%) [34]. A cohort study found that the progression of kidney dysfunction in T2D was faster in Indo-Asian (Indian, Pakistani, and Bangladeshi) subjects - with an estimated 2-3-fold increase in the mean rate of rise of serum creatinine over a constant follow-up period—compared to European-ancestry subjects [35]. The prevalence of diabetic end stage renal disease (ESRD) has also been reported as significantly higher in Asian T2D subjects (52.6%) compared to Caucasians (36.2%) [36].

Another microvascular complication of T2D, diabetic retinopathy, represents about five percent of all cases of global blindness [37]. Visual impairment occurs as a result of long-term, accumulated damage to small blood vessels in the retina. A recent cross-sectional study conducted by The Diabetic Retinopathy in Various Ethnic groups in UK (DRIVE UK) found that South Asian T2D populations have significantly higher prevalence of diabetic retinopathy (42.3% versus 38%) and sight threatening diabetic retinopathy (10.3% versus 5.5%) compared to white Europeans [38].

Combined with reduced blood flow, neuropathy (nerve damage) in the feet increases the risk of foot ulcers, infections, poor wound healing, and poor distal circulation, eventually increasing the risk of limb amputation [39]. Due to the elevated risk of these life-threatening complications, mortality risk among people with diabetes is at least double that of individuals without diabetes [40].

2.3. Conventional Risk Factors

A range of lifestyle and clinical factors contribute to risk of insulin resistance and T2D, including elevated body mass index (BMI), high waist-to-hip ratio (WHR), physical inactivity, and diet (Figure 1).

2.3.1. Body Mass Index (BMI) and Obesity

According to the World Health Organisation (WHO), body mass index (BMI) is a simple index of weight-for-height that can be widely used to classify overweight and obesity in adults [41]. It is defined as a person’s weight in kilograms divided by the square of their height in meters (kg/m2). Individuals with BMI greater than or equal to 30 kg/m2 are classified as obese for international standardised comparison. Obesity elevates serum fatty acid concentrations, reducing glucose uptake and increasing fatty acid uptake by the liver, skeletal muscle, and pancreatic β-cells. Reduced glucose uptake elevates serum glucose, stimulating further insulin secretion; it is the lack of response to this secreted insulin that induces insulin resistance [42]. Continually high insulin secretion in turn produces metabolic stress in pancreatic β-cell mitochondria, inducing the release of reactive oxygen species that damage mitochondria. Over time, mitochondria lose their ability to maintain cellular processes and β-cells undergo apoptosis, irreversibly reducing insulin secretion potential [43].

Associations between BMI, percentage of body fat, and body fat distribution differ across populations, influencing the thresholds at which T2D risk increases. Asian T2D patients have lower average BMI compared to European patients [44], which might reflect higher percentage body fat in Asians (3–5% higher) than Europeans for a given BMI [45, 46]. Similarly, for a fixed body fat percentage, Asians have a 3 to 4 unit lower BMI than Europeans [45]. The body fat percentage is also different between Asian groups; for fixed BMI, it tends to be the highest in Indians, followed by Malays and Chinese [47]. One study also showed that among Asians, Indians have the highest prevalence of obesity (35.8% (95% CI: 32.4–39.3)), followed by Malays (32.0% (95% CI: 30.6–33.4)) and Chinese (19.7 (95% CI: 17.9–21.6)) [13]. However, due to differences in body composition, recent studies have shown that waist circumference (WC) measurement or waist-to-hip ratio (WHR) is a better predictor of T2D in Asian populations than simple BMI or body fat percentage [48, 49], since these latter measures are insensitive to differences in body fat distribution.

2.3.2. Abdominal Obesity (High Waist-to-Hip Ratio/High Waist Circumference)

High waist-to-hip ratio (WHR) and waist circumference (WC), or abdominal obesity, is a major cause of insulin resistance since subcutaneous abdominal adipocytes drain their lipolytic products (free fatty acids) directly into the portal vein [50]. These free fatty acids are thought to decrease hepatic clearance of insulin and worsen systemic hyperinsulinemia [51], a precursor to T2D. Additional factors such as reduced secretion of adiponectin by adipose tissue may also contribute to the insulin-resistant state in individuals with abdominal obesity [52]. Adiponectin is an adipose tissue-specific protein that controls a number of metabolic processes, including insulin sensitivity and fatty acid oxidation [53].

The prevalence of abdominal obesity differs between ancestral groups and seems particularly marked in certain ethnic populations such as Native Americans, African-Americans, Asians, and Pacific Islanders [5456]. The Multi-Ethnicity Study of Atherosclerosis found that for a given waist circumference, Chinese have the highest diabetes incidence, followed by Hispanic, African, and European-ancestry individuals [57], a finding that may be explained by higher levels of visceral adipose tissue (VAT) in Chinese compared with Europeans, at a fixed waist circumference [58]. The same study also found that South Asians have substantially higher visceral adipose tissue compared to Europeans for given waist circumference [58]. This might explain increased lipid and insulin levels observed in South Asians compared with Europeans at the same WC and/or WHR [59]. Such differences are apparent not only between Asian and other continental populations, but also among Asian populations. Among three major Asian groups, the prevalence of abdominal obesity seems to be significantly higher among Indians (61.8% (95% CI: 58.3–65.2)) compared with Malays (45.3% (95% CI: 43.8–46.8) or Chinese (40.4% (95% CI: 38.0–42.7)) [13].

2.3.3. Diet and Physical Activity

The increasing global prevalence of T2D parallels escalating obesity rates resulting from reduced physical activity, increased intake of total calories, saturated fat (especially in fast food), and sugar-sweetened beverages (SSBs) in many societies. Asian populations are undergoing a nutrition transition in conjunction with the increasing adoption of Westernized lifestyles. In India and China, for example, caloric intake from animal fat has almost doubled in recent decades [60, 61]. High consumption of red and processed meat, SSBs, and refined grains with associated low consumption of cruciferous and yellow vegetables is strongly associated with increased in T2D [62]. At the same time, physical activity has reduced in Asian populations due to rapid urbanization and modernization [63, 64], further increasing T2D risk.

2.3.4. Metabolic Features

Metabolic features including elevated blood pressure, hyperglycaemia, and hyperlipidaemia increase T2D risk by several-fold [65]. A recent multi-ethnic population-based survey indicated population differences in the prevalence of metabolic syndrome features, irrespective of T2D status. Indians seem to have higher levels of triglycerides and hyperglycaemia and lower HDL cholesterol, compared with Malay and Chinese [66]. These findings parallel those of a case-control study in which Indians from UK and Indians from India had higher total insulin and triglycerides and lower HDL cholesterol compared to European individuals, irrespective of shared environmental influences [22].

2.3.5. Other Factors

Other factors that have been associated with T2D risk include short sleep duration [67, 68], increasing age, which may reflect reduced exercise and muscle mass [14], history of gestational diabetes, polycystic ovary syndrome, severe mental illness, and having a family history of the disease [54]. A recent randomized, crossover study found that sleep deprivation impairs peripheral metabolic pathways, thereby reducing insulin sensitivity [69]. The loss of skeletal muscle mass with age, or sarcopenia, is also related to insulin resistance, with sarcopenia thought to cause insulin resistance and thereby increase risk of diabetes [70]. In turn, insulin resistance results in further loss of muscle strength [71]. Finally, patients with severe mental illness such as schizophrenia or bipolar disorder have 3-fold higher risk of developing T2D compared to the general population; this may result from underlying lifestyle factors, adverse effects of pharmacotherapy, and possible common genetic and/or environmentally linked pathophysiologic processes [72].

2.4. Genetic Susceptibility

In addition to conventional risk factors, family, twin and genetic studies show that T2D susceptibility has a substantial genetic component [73]. Full siblings of T2D probands have a 30–60% increased risk of disease, compared with the general population [74, 75] and children with one affected parent have a 40% lifetime risk of developing T2D, which rises to almost 70% if both parents are affected [76]. Twin studies also show higher T2D concordance in monozygotic (60–70%) compared with dizygotic twins (20–30%) [7779].

The proportion of trait variance due to additive genetic effects is termed “heritability” and can be formally estimated from twin studies. Twin study heritability estimates are on the order of 30–70% for T2D and about 60% for impaired glucose tolerance (IGT) [80, 81]. Twin studies also demonstrate a substantial genetic component for quantitative phenotypes related to glucose homeostasis, with heritability estimates of 75–85% for in vivo insulin secretion, ~50% for peripheral insulin sensitivity, and ~50% for glucose metabolism [82].

Population differences in T2D pathophysiologic and risk factor profiles have been discussed in previous sections. It has been suggested that such differences may partly reflect population differences in the frequency of particular genetic risk factors and/or population-specific interactions between genetic and environmental factors [83].

2.5. Methods of Gene Identification for Common Complex Disease

Observed patterns of T2D inheritance, combined with the results of recent large-scale genetic studies, suggest that the genetic component of T2D is complex, involving multiple genetic variants of individually small effect (polygenic model) [84]. There have been three main approaches employed to identify genetic risk variants for such common complex disorders: linkage studies, candidate gene association studies (CGAS), and, more recently, genome-wide association studies (GWAS).

2.5.1. Linkage Studies

Familial linkage studies seek to identify broad genomic regions harboring disease risk variants by tracking disease and genetic marker segregation through multiple generations of families. Familial linkage studies are challenging for disorders with advanced age at onset, as parents may no longer be alive. Further challenges include difficulty in collecting accurate genealogical information and genetic (locus) heterogeneity, meaning that a particular risk locus contributes to disease in only a subset of families [85]. More broadly, this approach is limited by low power for common variants of small effect [86] and its inability to precisely localise underlying risk variants [87]. Earlier linkage studies found four (4) genetic loci linked with T2D; CAPN10 [88], ENPP1 [89], HNF4A [90], and ACDC (ADIPOQ) [91]. However, only the HNF4A locus has been confirmed by recent large-scale genome-wide association studies (GWAS) [92].

HNF4A, together with the related locus, HNF1A and also GCK also account for up to 80% of rare monogenic forms of diabetes. These diabetes cases present as familial, young onset, noninsulin dependent diabetes mellitus (maturity onset diabetes of the young or MODY) and are inherited in a Mendelian dominant pattern [75]. Unlike common polygenic forms of T2D, these monogenic forms require only one pathogenic genetic variant to produce disease.

2.5.2. Candidate Gene Association Studies (CGAS)

In contrast to linkage studies, the candidate gene approach searches for associations between common genetic variants and disease, restricting the search region to prespecified genes of interest. Candidate genes are typically selected based on a priori knowledge or hypotheses reflecting the gene’s presumed biological role in disease [93]. The most common study design is the case control study; for a particular genetic variant, this involves comparing the frequency of genetic alleles between individuals with and without disease, aiming to identify alleles that are associated with disease status [87]. Although a mainstay of the initial era of disease gene mapping, the candidate gene approach has been limited by small sample sizes, restriction of hypotheses to known biology, and an inability to replicate many results [94]. While candidate gene studies have reported numerous variants as beeing associated with T2D [75], just three loci, PPARG [95], KCNJ11 [96], and TCF7L2 [97], have been robustly confirmed by recent GWAS [74, 98, 99]. We note that the TCF7L2 association study was informed by prior genome-wide linkage study showing suggestive linkage between T2D and the 10q genomic region harboring TCF7L2 [97].

2.5.3. Genome-Wide Association Studies (GWAS)

Within the last five years, genome-wide association studies (GWAS) have emerged as the method of choice for identifying common genetic variants associated with complex disease. GWAS were facilitated by completion of the Human Genome Project in 2003, the International HapMap Project in 2005 that catalogued millions of common single nucleotide polymorphisms (SNPs), and the parallel development of high throughput genotyping arrays. Single nucleotide polymorphisms (SNPs) are DNA sequence variants in which a single nucleotide differs between individuals. SNPs have a low historical mutation rate, are amenable to high throughput genotyping, and are distributed abundantly across all 22 autosomes and the sex chromosomes. Typically several million variants are screened genome-wide; appropriate adjustment of the prespecified significance level is thus necessary to avoid an increase in false positive results. Based on patterns of human genomic correlation, Bonferroni correction for one million independent tests is the accepted approach, with variants required to reach a pointwise value < 5 × 10−8 (or 0.05/1,000,000) to be declared “genome-wide significant” [100]. Due to this stringent significance level, very large sample sizes are necessary to identify associations of modest effect, which is often achieved via international collaboration and the formation of consortia. The existence of such collaborations also facilitates rapid replication of findings in independent samples, a requirement for publication.

It has been eight years since the first notable GWAS finding in 2005, identifying a common allele of large effect associated with age-related macular degeneration. The year 2007 was coined the “Year of GWAS”, due to the explosion of GWAS publications in that year. From 2005 to September 2013, there have been more than 1,600 GWAS published reports for a range of human diseases and traits, with an online catalogue established by the National Human Genome Research Institute at the US National Institutes of Health to collate major findings (http://www.genome.gov/gwastudies/). Although complicated and costly, GWAS have successfully identified thousands of genetic loci associated with common complex diseases under the common disease common variant (CDCV) hypothesis.

2.6. Genome-Wide Association Studies of T2D

The first T2D GWAS was published in 2007 [99], and as of September 2013, there were more than 40 publications on T2D and its complications listed in the online catalogue of published genome-wide association studies (http://www.genome.gov/gwastudies/). At the time of writing, the catalogue describes over 100 individual SNPs showing genome-wide significant association ( < 5 × 10−8) with T2D and related metabolic traits across diverse populations (Table 1) and over 60 SNPs showing suggestive association ( < 1 × 10−5) (Table 2). This section will provide a review of T2D GWAS findings to date.

The first T2D GWAS was conducted in European-ancestry participants 2007 by Sladek and colleagues [99], using a discovery sample of 600 cases and 600 controls and a separate European replication sample of nearly 3,000 cases and 3,000 controls. This small study of early onset T2D reported T2D-associated variants in three novel susceptibility genes: TCF7L2 and HHEX/IDE which are associated with β-cell function [101] and SLC30A8, encoding a zinc transporter highly expressed in pancreatic islets [102]. Several months later, four additional European studies [74, 98, 103, 104] confirmed association of variants at these loci and identified additional associated variants in IGF2BP2, associated with β-cell dysfunction [105], and CDKN2A/CDKN2B and CDKAL1, which are both associated with β-cell development [105, 106]. During this time, variants in the FTO (fat mass and obesity associated) gene were also identified with important effects on obesity and hence, indirectly, T2D [107, 108]. Interestingly, as the effect of FTO variants on T2D is only via obesity, the FTO locus was not identified in T2D GWAS using cases and controls matched for BMI. Two of these early publications also showcased the output of international consortia: The UK-based Wellcome Trust Case Control Consortium (WTCCC) and the USA-based Diabetes Genetics Initiative (DGI), highlighting the benefits of large-scale collaboration in the GWAS era.

Since these initial GWAS had modest power to detect variants with modest effects on disease risk, follow-up studies employed meta-analysis to increase sample size and hence power to detect additional loci of similar or smaller effect. The first T2D GWAS meta-analysis was published in 2008 and was also a European study [109], representing collaboration between three different consortia; the WTCCC, DGI, and the Finland—United States Investigation of NIDDM Genetics (FUSION) which combined to form the Diabetes Genetics Replication and Meta-Analysis (DIAGRAM) consortium. This study utilised an enlarged discovery sample of 4,549 cases and 5,579 controls with replication in further 24,194 cases and 55,598 controls, all of European-ancestry. This study identified associated variants at six additional novel loci: JAZF1, CDC123, TSPAN8 and THADA which are associated with β-cell dysfunction [110, 111], ADAMTS9 which is associated with insulin action [111], and NOTCH2, associated with glucose-stimulated glucagon secretion by pancreatic islet cells [112].

These initial T2D GWAS were all restricted to populations of European-ancestry. The first two large-scale T2D GWAS conducted in Asian populations were reported in 2008, each employing a multi-stage design in East Asian groups. Both studies reported association of variants in KCNQ1, encoding the alpha subunit of a voltage-gated potassium channel expressed in the pancreas [113, 114]. These studies clearly demonstrated the utility of extending T2D GWAS to non-European populations; association of the KCNQ1 variants with T2D was not detected in previous European GWAS, due to markedly lower frequency of the risk allele in Europeans (5% versus 40%), resulting in dramatically reduced power. A European meta-analysis subsequently confirmed association of the KCNQ1 variants with T2D in Europeans but at significance levels far below thresholds usually inspiring replication or follow-up studies ().

A European study published in 2009 used multiple samples of French, Danish, and Finnish ancestry to identify association of variants in the insulin receptor substrate 1 gene (IRS1), showing that the risk allele is also associated with insulin resistance and hyperinsulinaemia in large population-based cohorts [115]. This contrasted with the apparent biology of previous associations, which principally related to impaired pancreatic β-cell function.

This first wave of T2D GWAS was succeeded by a second wave beginning in 2010, in which existing and new datasets were combined into expanded meta-analyses. The most notable was a large European study reported by Voight and colleagues, involving ~42,000 T2D cases and 100,000 controls split between discovery and replication stages and identifying twelve new associated loci. These included X-chromosomal association and an HNFA1A locus overlapping with the locus implicated in Mendelian monogenic (single gene) forms of diabetes [116]. Other studies reported at this time included three East Asian studies, one African American, and one European study, which together identified nine (9) additional loci [117121]. Three of these genes have unknown function (RBMS1, PTPRD, and SRR) [117, 118], while RPS12, LIMK2, and AUH are associated with diabetic nephropathy [121], C2CD4A is associated with β-cell dysfunction [119, 122], SPRY2 is associated with obesity and insulin resistance [120, 123], and SASH1 is associated with insulin growth factors [121].

A subsequent 2011 meta-analysis included three Southeast Asian populations: Chinese (3955 subjects), Indian (2146 subjects), and Malay (2034 subjects) and it further emphasized the value of surveying diverse ethnic groups [124]. This study was the first to include individuals from the Malay population, the largest group in Southeast Asia, with a population size of more than 300 million [124]. This study alone contributed an additional 16 novel loci, in spite of its relatively modest sample size; this partly reflected higher minor allele frequencies in Southeast Asian populations at some associated loci (e.g., rs3792615, number 18 in Table 2).

The first T2D GWAS in South Asian populations was also published in 2011, including individuals from India, Pakistan, Sri Lanka and Bangladesh. Using a relatively modest sample size (5,561 cases and 14,458 controls in the discovery step) five additional novel T2D loci were discovered [92]: HNF4A, involved in monogenic forms of diabetes and associated with β-cell development [125]; GRB14 which is associated with obesity and insulin resistance [126]; and another three loci with less clear functions; AP3S2, ST6GAL1 and VPS26A [92].

Another large meta-analysis in East Asian groups were performed in 2012 and identified 10 further novel loci [127] with mostly unknown function except for GLIS3, associated with β-cell development [128]. It is interesting that the MAEA variant discovered in this study is unique to East Asian and African populations, being monomorphic in Europeans and South Asians (Table 1; number 24) [129]. Several other variants identified in this study have substantially higher risk allele frequency (RAF) in East Asians than Europeans, for example, ZFAND3 (Table 1; number 28; 34% versus 12%), FITM2-R3HDML (Table 1; number 111; 41% versus 18%), and RPS3P7-MAF (Table 2; number 56, 18% versus 1%), enhancing their detection in East Asian samples of relatively modest size.

Reflecting the success of initial T2D GWAS and the fast pace of technology, in 2012 Voight [130] and colleagues designed the “Metabochip,” a custom genotyping array enriched for variants shown to be associated with cardiometabolic traits via GWAS. These traits include T2D, coronary heart disease, myocardial infarction, body mass index, glucose and insulin level, lipid levels, and blood pressure. This new platform offers a powerful and cost-effective approach to both the discovery and follow-up of variants associated with these related traits, due to comprehensive coverage of previously associated loci (~120,000 SNPs) [130]. Morris and colleagues used the Metabochip to discover and characterize T2D-associated variants via meta-analysis of 34,840 cases and 114,981 controls of European descent, finding ten novel loci [131] not reported in previous European studies, all of which reached genome-wide significance. Another study using the Metabochip combined newly available samples with samples from previous discovery meta-analyses, using genotype data for 66,000 follow-up SNPs. This study identified 41 novel glycaemic associations, 33 of which were also associated with T2D [132]. This study implicated new loci in the aetiology of T2D and increased the overlap between loci associated with both glycaemia and T2D. These studies highlight the Metabochip as a promising tool for identifying novel and robustly associated loci, facilitating future research into underlying biology.

Taken together, the results of T2D GWAS signify tremendous progress in our quest to understand the genetic causes of T2D. Alternatively, they also highlight the genetic complexity of this disease. Genetic variants showing replicable association with T2D uniformly exert only a modest effect on disease risk, with per-allele odds ratios typically in the range of 1.1–1.3 (Table 1). The combined effect of all variants reported to date explains only about 10% of observed familial clustering [116]. Furthermore, the functional significance of various loci remains unclear. While some appear to be associated with β-cell function and insulin resistance, the biological role of many of them remains unknown. This suggests that the findings to date represent the first stage of a long journey to understanding T2D genetic risk.

2.7. Polygenic Models of T2D

The distribution of odds ratios observed for validated T2D-associated SNPs suggests that numerous, associated loci exist with even smaller effects than those identified to date. One would not expect such loci to have reached genome-wide significance in previous GWAS due to insufficient power. The existence of such additional small effect loci is consistent with the pattern of additional associated variants being discovered as sample sizes have increased; it is also consistent with validated SNP associations explaining only a small proportion of the T2D heritability estimated from twin studies, known as the “missing heritability” problem.

Two methods have recently been developed for assessing the contribution of common SNPs not reaching genome-wide significance to the heritability of a trait. These are polygenic scoring [133] and mixed linear modelling [134]. Both methods test the combined effects of thousands (or hundreds of thousands) of SNPs upon a trait of interest. A recent study by Stahl and colleagues used polygenic analyses and linear mixed modelling to show that thousands of SNPs contribute to T2D risk, estimating that about 50% of observed variance in T2D risk could be attributed to the combined effects of all SNPs genome-wide [135]. These investigators suggested that at least 70% of T2D heritability can be attributed to common SNPs represented on GWAS arrays [135], with most having very small individual effects upon disease risk.

2.8. Population Differences in T2D Risk Alleles

The frequency of T2D risk alleles often varies between populations, producing population differences in the attributable risk due to a particular genetic risk factor or combination of risk factors. The discovery of KCNQ1 emphasized the impact of such frequency differences upon genetic discoveries [136]. Association of KCNQ1 variants was found in East Asian populations [113, 114] using a considerably smaller sample size than that required to detect the association in with European populations [116], reflecting higher allele frequency (33% in East Asian versus 8% in Europeans) and hence statistical power in Asian groups. In addition, variants at the TCF7L2 locus showed the inverse; a high risk allele frequency in Europeans (30%) compared to a low frequency in East Asians (3%) enhanced the detection in European studies [137]. Similarly, the SYK variant demonstrates a RAF of 26% in East Asians [138]and only 2% in Europeans and is monomorphic in Africans (Table 2; number 36). Further, a number of T2D risk variants are monomorphic (not variable) in some populations, preventing the detection of an association in these groups. The recently reported SCGG variant is unique to Indian Punjabi Sikh, being monomorphic in both European and African populations (Table 1; number 92). Other instances include the THADA variant, which was discovered in European populations but is monomorphic in East Asians (Table 1; number 4), RND3-FABP5P10 that was discovered in African Americans but is monomorphic in Europeans (Table 1; number 6), and G6PC2, discovered in Europeans but monomorphic in Africans (Table 1; number 10). For a set of SNPs showing association with T2D across multiple populations, Table 3 shows risk allele frequencies and odds ratios for different populations in which associations have been reported. For these 14 SNPs, risk allele frequencies commonly differ across populations; however, allelic effects upon disease seem markedly consistent in both direction and magnitude, given overlapping confidence intervals for allelic odds ratios. Taken together, these results suggest that population differences can have important effects on power to detect common genetic associations for T2D in samples of diverse ancestry but may have less impact upon disease risk within individuals carrying the identified risk alleles. Nevertheless, at the population level, the attributable risk of such genetic variants will increase with allele frequency, thus potentially influencing population disease burden.

Significantly, a recent study assessing thousands of genetic associations showed that T2D has the most extreme population differentiation of risk allele frequencies among a broad range of complex diseases [139]. T2D risk allele frequencies demonstrated clear gradient matching paths of early human migration, suggesting potential differences in evolutionary adaptation to food availability, dietary patterns, or agricultural practices. This is consistent with “thrifty genotype” hypothesis [139, 140], which proposes that susceptibility to obesity and T2D in some populations reflects historical positive selection for genotypes promoting efficiency of metabolism, and energy and fat storage, thus providing an advantage in times of nutrient shortage [141]. This might explain the extraordinarily high prevalence of diabetes now seen among certain populations [34, 142, 143], potentially reflecting historical feast and famine cycles [62], increasing the frequency of thrifty genotypes and genetic predisposition to obesity and diabetes. While being unproven, this may partly explain higher susceptibility to abdominal obesity at lower BMI and reduced muscle mass with increased insulin resistance in Asian compared with Caucasian populations [7]. Nevertheless, pronounced population differentiation of T2D risk allele frequencies provides a strong rationale for further comprehensive genetic studies of T2D in diverse populations, expanding on the comprehensive studies in European samples.

To date, a range of non-European T2D GWAS have been conducted, including studies in Japanese [114, 119, 138, 144], Chinese [117, 145, 146], African-American [121, 147, 148], Southeast Asian [124], Hispanic [149], Mexican-American [150], South Asian [92], Indo-European [151], and Indian Punjabi Sikh [152]. These studies have led to new discoveries, including novel loci and loci that seem specific to certain populations [119, 127, 151, 152]. While many loci appear to contribute broadly to T2D risk, some loci have currently been confirmed in European populations only, including WFS1, NOTCH2, THADA, ADAMTS9, TSPAN8/LGR5, INS-IGF2, ADCY5, GCK, MTRNR1B, HMGA2, HNF1A, ZBED3, KLF14, ZFAND6, PRC1, TLES/CHCHD9, and RBMS1 [109, 116, 153155]. Other loci currently show association specifically in East Asian populations, including PTPRD, SRR, CDC123/CAMK1D, PSMD6, MAEA, ZFAND3, KCNK16, GCC1/PAX4, GLIS3, and PEPD [117, 119, 120, 127]. On the other hand, TMEM163 [151] and SGCG [152] appear unique to South Asian and Indian Punjabi Sikh, respectively. Some of these discoveries may reflect the impact of population allele frequency differences, as previously discussed. In such cases, larger studies may eventually show that some loci contribute to disease across a broader range of populations.

Seemingly population-specific genetic associations for T2D may also reflect differences in the patterns of genomic correlation, or linkage disequilibrium (LD), between associated marker loci and the underlying unobserved functional variants. Populations with different demographic histories will often display different patterns of LD reflecting population differences in evolutionary recombination [156]. Older populations such as those in Africa have lower LD and can be helpful for finely localizing a risk variant following an initial association finding. This is because the genomic distance between disease-associated markers and true risk variants is likely to be smaller in such populations [157].

Thus, the apparent population-specificity of some known T2D risk alleles may reflect population differences in risk allele frequencies or LD between tagging and causal variants, rather than actual population-specificity of the underlying functional risk loci. We note that population-specific estimates of disease variance explained by all known T2D loci—although not widely reported—do seem largely similar between European and Asian populations. In their large 2012 study, Morris and colleagues [131] estimated that known common variants explain about 5.7% of T2D disease variance in Europeans. In 2013, Tabassum and colleagues [151] estimated that known loci combined with one novel Indian-specific locus explained 7.65% of T2D risk variance in South Asian Indians. The slightly higher estimate in Indians may potentially be explained by the inclusion of additional variants discovered between publications of the two studies, together with the inclusion of the Indian-specific locus discovered in the Tabassum study. Thus, available evidence thus does not strongly suggest that differences in the cumulative variance explained by known common T2D risk alleles can explain the markedly higher T2D prevalence observed in South Asians.

Interestingly, however, very recent studies show that population differences in linkage disequilibrium (LD) and the presence of multiple independent variants within a locus can markedly influence estimates of variance explained by known risk variants [158, 159]. Detailed fine mapping of T2D susceptibility loci in diverse populations, combined with the identification of underlying functional variants, may thus reveal population differences in the contribution of known loci to disease. Future research may also show the extent to which population differences in T2D risk can be explained by rare alleles, gene-environment interactions, or epigenetic effects.

2.9. Gene-Environment Interactions in T2D

In addition to the effects of specific genetic and environmental risk factors, gene-environment interactions are likely important mediators of population differences in T2D risk and contributors to the “missing heritability” problem. Indeed, given the relative stability of DNA sequence within populations over decades, recent increases in T2D prevalence must largely reflect environmental changes. Accordingly, the single largest contributor to T2D risk is obesity, and the global T2D epidemic chronologically parallels the global obesity epidemic.

A paucity of studies has examined gene-environment interactions in the context of T2D in general, let alone in Asian populations. A study by Qi and colleagues [160] found that a high genetic risk score formed from 10 T2D-associated SNPs was further increased by the presence of a “Westernized” dietary pattern characterised by increased red and processed meat intake and reduced dietary fibre [160]. The Westernized diet was not associated with increased risk among those with a low genetic risk score. Several studies have also found evidence for interactions between T2D-associated variants in TCF7L2 and the quality and quantity of ingested carbohydrates in the context of T2D risk [161163]. These studies support a possible contribution of gene-environment interactions to T2D risk, together with a potential model where interactions between recent lifestyle transitions and genetic risk factors may be contributing to the rapidly increasing prevalence of T2D in Asian populations. However, these preliminary findings require validation. Future analyses in well-designed, well-powered studies will help to clarify the role of gene-environment interactions in population differences in T2D risk.

2.10. Epigenetics

Similar to the “thrifty genotype” hypothesis, the “thrifty phenotype” hypothesis considers the adaptive consequences of the environment in utero. The hypothesis relates to the metabolic consequences of fetal malnutrition, proposing that adaptation to a low-calorie intrauterine environment induces permanent changes in chromatin structure and gene expression that influence insulin secretion and resistance, promoting more efficient energy utilisation and thus fetal survival [164]. According to the hypothesis, such epigenetic changes may predispose to insulin dysregulation, obesity, and T2D in later life. In support, epidemiologic studies have shown associations between small birth size, a marker for fetal malnutrition, and adult-onset T2D [165, 166]. A study by van Hoek and colleagues [167] in the Dutch Famine Birth Cohort detected an interaction between an IGF2BP2 polymorphism and prenatal famine upon glucose level in the offspring. Interactions between other T2D risk variant alleles and birthweight have also been associated with increased T2D risk [168, 169].

3. Conclusions

We have discussed differences in prevalence, risk factor profiles, and genetic risk allele frequencies between different Asian countries and between Asian and other continental populations. Given these differences, continued T2D genetic studies in diverse populations are likely to contribute crucially to the broadening terrain of shared and unique population genetic effects for this disorder. Future studies will ideally include large, population-specific characterisation of risk variants, studies of gene-environment interaction, and epigenetic studies. Well-powered, well-designed studies performed in diverse Asian populations should enhance the benefits of genetic discoveries and their ultimate clinical translation for these large susceptible groups.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

Noraidatulakma Abdullah is sponsored by Ministry of Higher Education of Malaysia and Universiti Kebangsaan Malaysia. Elizabeth Holliday is funded by the NHMRC fellowship scheme.