Identification of novel high-impact recessively inherited type 2 diabetes risk variants in the Greenlandic population

Genome-wide association study and functional characterization identifies candidate genes for insulin-stimulated glucose uptake

Distinct tissue-specific mechanisms mediate insulin action in fasting and postprandial states. Previous genetic studies have largely focused on insulin resistance in the fasting state, where hepatic insulin action dominates. Here we studied genetic variants influencing insulin levels measured 2 h after a glucose challenge in >55,000 participants from three ancestry groups. We identified ten new loci (P < 5 × 10-8) not previously associated with postchallenge insulin resistance, eight of which were shown to share their genetic architecture with type 2 diabetes in colocalization analyses. We investigated candidate genes at a subset of associated loci in cultured cells and identified nine candidate genes newly implicated in the expression or trafficking of GLUT4, the key glucose transporter in postprandial glucose uptake in muscle and fat. By focusing on postprandial insulin resistance, we highlighted the mechanisms of action at type 2 diabetes loci that are not adequately captured by studies of fasting glycemic traits.

A novel splice-affecting HNF1A variant with large population impact on diabetes in Greenland

Background

The genetic disease architecture of Inuit includes a large number of common high-impact variants. Identification of such variants contributes to our understanding of the genetic aetiology of diseases and improves global equity in genomic personalised medicine. We aimed to identify and characterise novel variants in genes associated with Maturity Onset Diabetes of the Young (MODY) in the Greenlandic population.

Methods

Using combined data from Greenlandic population cohorts of 4497 individuals, including 448 whole genome sequenced individuals, we screened 14 known MODY genes for previously identified and novel variants. We functionally characterised an identified novel variant and assessed its association with diabetes prevalence and cardiometabolic traits and population impact.

Findings

We identified a novel variant in the known MODY gene HNF1A with an allele frequency of 1.9% in the Greenlandic Inuit and absent elsewhere. Functional assays indicate that it prevents normal splicing of the gene. The variant caused lower 30-min insulin (β = -232 pmol/L, β_SD = -0.695, P = 4.43 × 10^-4) and higher 30-min glucose (β = 1.20 mmol/L, β_SD = 0.441, P = 0.0271) during an oral glucose tolerance test. Furthermore, the variant was associated with type 2 diabetes (OR 4.35, P = 7.24 × 10^-6) and HbA1c (β = 0.113 HbA1c%, β_SD = 0.205, P = 7.84 × 10^-3). The variant explained 2.5% of diabetes variance in Greenland.

Interpretation

The reported variant has the largest population impact of any previously reported variant within a MODY gene. Together with the recessive TBC1D4 variant, we show that close to 1 in 5 cases of diabetes (18%) in Greenland are associated with high-impact genetic variants compared to 1-3% in large populations.

Funding

Novo Nordisk Foundation, Independent Research Fund Denmark, and Karen Elise Jensen's Foundation.

Recessive Genome-Wide Meta-analysis Illuminates Genetic Architecture of Type 2 Diabetes

Most genome-wide association studies (GWAS) of complex traits are performed using models with additive allelic effects. Hundreds of loci associated with type 2 diabetes have been identified using this approach. Additive models, however, can miss loci with recessive effects, thereby leaving potentially important genes undiscovered. We conducted the largest GWAS meta-analysis using a recessive model for type 2 diabetes. Our discovery sample included 33,139 case subjects and 279,507 control subjects from 7 European-ancestry cohorts, including the UK Biobank. We identified 51 loci associated with type 2 diabetes, including five variants undetected by prior additive analyses. Two of the five variants had minor allele frequency of <5% and were each associated with more than a doubled risk in homozygous carriers. Using two additional cohorts, FinnGen and a Danish cohort, we replicated three of the variants, including one of the low-frequency variants, rs115018790, which had an odds ratio in homozygous carriers of 2.56 (95% CI 2.05-3.19; P = 1 × 10-16) and a stronger effect in men than in women (for interaction, P = 7 × 10-7). The signal was associated with multiple diabetes-related traits, with homozygous carriers showing a 10% decrease in LDL cholesterol and a 20% increase in triglycerides; colocalization analysis linked this signal to reduced expression of the nearby PELO gene. These results demonstrate that recessive models, when compared with GWAS using the additive approach, can identify novel loci, including large-effect variants with pathophysiological consequences relevant to type 2 diabetes.

Identification of Recessively Inherited Genetic Variants Potentially Linked to Pancreatic Cancer Risk

Although 21 pancreatic cancer susceptibility loci have been identified in individuals of European ancestry through genome-wide association studies (GWASs), much of the heritability of pancreatic cancer risk remains unidentified. A recessive genetic model could be a powerful tool for identifying additional risk variants. To discover recessively inherited pancreatic cancer risk loci, we performed a re-analysis of the largest pancreatic cancer GWAS, the Pancreatic Cancer Cohort Consortium (PanScan) and the Pancreatic Cancer Case-Control Consortium (PanC4), including 8,769 cases and 7,055 controls of European ancestry. Six single nucleotide polymorphisms (SNPs) showed associations with pancreatic cancer risk according to a recessive model of inheritance. We replicated these variants in 3,212 cases and 3,470 controls collected from the PANcreatic Disease ReseArch (PANDoRA) consortium. The results of the meta-analyses confirmed that rs4626538 (7q32.2), rs7008921 (8p23.2) and rs147904962 (17q21.31) showed specific recessive effects (p<10-5) compared with the additive effects (p>10-3), although none of the six SNPs reached the conventional threshold for genome-wide significance (p < 5×10-8). Additional bioinformatic analysis explored the functional annotations of the SNPs and indicated a possible relationship between rs36018702 and expression of the BCL2L11 and BUB1 genes, which are known to be involved in pancreatic biology. Our findings, while not conclusive, indicate the importance of considering non-additive genetic models when performing GWAS analysis. The SNPs associated with pancreatic cancer in this study could be used for further meta-analysis for recessive association of SNPs and pancreatic cancer risk and might be a useful addiction to improve the performance of polygenic risk scores.

The Opportunities and Challenges of Integrating Population Histories Into Genetic Studies for Diverse Populations: A Motivating Example From Native Hawaiians

There is a well-recognized need to include diverse populations in genetic studies, but several obstacles continue to be prohibitive, including (but are not limited to) the difficulty of recruiting individuals from diverse populations in large numbers and the lack of representation in available genomic references. These obstacles notwithstanding, studying multiple diverse populations would provide informative, population-specific insights. Using Native Hawaiians as an example of an understudied population with a unique evolutionary history, I will argue that by developing key genomic resources and integrating evolutionary thinking into genetic epidemiology, we will have the opportunity to efficiently advance our knowledge of the genetic risk factors, ameliorate health disparity, and improve healthcare in this underserved population.

Clinical significance of LARGE1 in progression of liver cancer and the underlying mechanism

Liver cancer is a malignant disease and causes thousands of death each year. The prognosis is dismal for patients with metastasis and recurrence. It is urgent to disclose the cause and mechanism underlying liver cancer. LARGE1 encodes a glycosyltransferase and was reported to promote progression in cancer. But its role in liver cancer is unknown. In this study, LARGE1 displayed upregulated expression in liver cancer cells. When LARGE1 was knocked down in SMMC-7721 and Huh-7 cells, the ability of cell proliferation and colony formation were decreased significantly. Cell migration and invasion were suppressed. The number of cells in G1 phase increased but decreased in S phase. Cell apoptosis was not affected. Tumor growth in vivo was also inhibited. Tumor volume was decreased from 1270 mm3 to 721 mm3 (p < 0.05) and tumor weight from 0.95 g to 0.63 g (p < 0.05). Furthermore, the expression of β-catenin, TCF and Cyclin D1 was reduced when LARGE1 was knocked down but increased in LARGE1-overexpressed cells. LGK-974, a specific inhibitor in canonical Wnt signaling, inhibited cell proliferation even when LARGE1 was over-expressed. In tumor tissues, LARGE1 was increased by 4.8 folds compared to paratumoral tissues. And higher LARGE1 expression caused shorter survival. Clinicopathological analysis demonstrated that LARGE1 was associated with TNM stage (Ⅰ/Ⅱ vs III/IV, p = 0.005). Therefore, LARGE1 promotes progression and regulates Wnt/β-catenin signaling pathway in liver cancer.

The impact of non-additive genetic associations on age-related complex diseases

Genome-wide association studies (GWAS) are not fully comprehensive, as current strategies typically test only the additive model, exclude the X chromosome, and use only one reference panel for genotype imputation. We implement an extensive GWAS strategy, GUIDANCE, which improves genotype imputation by using multiple reference panels and includes the analysis of the X chromosome and non-additive models to test for association. We apply this methodology to 62,281 subjects across 22 age-related diseases and identify 94 genome-wide associated loci, including 26 previously unreported. Moreover, we observe that 27.7% of the 94 loci are missed if we use standard imputation strategies with a single reference panel, such as HRC, and only test the additive model. Among the new findings, we identify three novel low-frequency recessive variants with odds ratios larger than 4, which need at least a three-fold larger sample size to be detected under the additive model. This study highlights the benefits of applying innovative strategies to better uncover the genetic architecture of complex diseases.

Ethnic differences in adiposity and diabetes risk - insights from genetic studies

Type 2 diabetes is more common in non-Europeans and starts at a younger age and at lower BMI cut-offs. This review discusses the insights from genetic studies about pathophysiological mechanisms which determine risk of disease with a focus on the role of adiposity and body fat distribution in ethnic disparity in risk of type 2 diabetes. During the past decade, genome-wide association studies (GWAS) have identified more than 400 genetic variants associated with the risk of type 2 diabetes. The Eurocentric nature of these genetic studies has made them less effective in identifying mechanisms that make non-Europeans more susceptible to higher risk of disease. One possible mechanism suggested by epidemiological studies is the role of ethnic difference in body fat distribution. Using genetic variants associated with an ability to store extra fat in a safe place, which is subcutaneous adipose tissue, we discuss how different ethnic groups could be genetically less susceptible to type 2 diabetes by developing a more favourable fat distribution.

The derived allele of a novel intergenic variant at chromosome 11 associates with lower body mass index and a favorable metabolic phenotype in Greenlanders

The genetic architecture of the small and isolated Greenlandic population is advantageous for identification of novel genetic variants associated with cardio-metabolic traits. We aimed to identify genetic loci associated with body mass index (BMI), to expand the knowledge of the genetic and biological mechanisms underlying obesity. Stage 1 BMI-association analyses were performed in 4,626 Greenlanders. Stage 2 replication and meta-analysis were performed in additional cohorts comprising 1,058 Yup'ik Alaska Native people, and 1,529 Greenlanders. Obesity-related traits were assessed in the stage 1 study population. We identified a common variant on chromosome 11, rs4936356, where the derived G-allele had a frequency of 24% in the stage 1 study population. The derived allele was genome-wide significantly associated with lower BMI (beta (SE), -0.14 SD (0.03), p = 3.2x10-8), corresponding to 0.64 kg/m2 lower BMI per G allele in the stage 1 study population. We observed a similar effect in the Yup'ik cohort (-0.09 SD, p = 0.038), and a non-significant effect in the same direction in the independent Greenlandic stage 2 cohort (-0.03 SD, p = 0.514). The association remained genome-wide significant in meta-analysis of the Arctic cohorts (-0.10 SD (0.02), p = 4.7x10-8). Moreover, the variant was associated with a leaner body type (weight, -1.68 (0.37) kg; waist circumference, -1.52 (0.33) cm; hip circumference, -0.85 (0.24) cm; lean mass, -0.84 (0.19) kg; fat mass and percent, -1.66 (0.33) kg and -1.39 (0.27) %; visceral adipose tissue, -0.30 (0.07) cm; subcutaneous adipose tissue, -0.16 (0.05) cm, all p<0.0002), lower insulin resistance (HOMA-IR, -0.12 (0.04), p = 0.00021), and favorable lipid levels (triglyceride, -0.05 (0.02) mmol/l, p = 0.025; HDL-cholesterol, 0.04 (0.01) mmol/l, p = 0.0015). In conclusion, we identified a novel variant, where the derived G-allele possibly associated with lower BMI in Arctic populations, and as a consequence also leaner body type, lower insulin resistance, and a favorable lipid profile.

ALOX5, LPA, MMP9 and TPO gene polymorphisms increase atherothrombosis susceptibility in middle-aged Mexicans

Atherothrombosis is the cornerstone of cardiovascular diseases and the primary cause of death worldwide. Genetic contribution to disturbances in lipid metabolism, coagulation, inflammation and oxidative stress increase the susceptibility to its development and progression. Given its multifactorial nature, the multiloci studies have been proposed as potential predictors of susceptibility. A cross-sectional study was conducted to explore the contribution of nine genes involved in oxidative stress, inflammatory and thrombotic processes in 204 subjects with atherothrombosis matched by age and gender with a healthy group (n = 204). To evaluate the possibility of spurious associations owing to the Mexican population genetic heterogeneity as well as its ancestral origins, 300 unrelated mestizo individuals and 329 Native Americans were also included. ALOX5, LPA, MMP9 and TPO gene polymorphisms, as well as their multiallelic combinations, were twice to four times more frequent in those individuals with clinical manifestations of atherothrombosis than in the healthy group. Once adjusting for population stratification was done, these differences remained. Our results add further evidence on the contribution of ALOX5, LPA, MMP9 and TPO polymorphisms to atherothrombosis development in the middle-aged group, emphasizing the multiethnic studies in search of gene risk polymorphisms.

Point mutations in the PDX1 transactivation domain impair human β-cell development and function

OBJECTIVE

Hundreds of missense mutations in the coding region of PDX1 exist; however, if these mutations predispose to diabetes mellitus is unknown.

METHODS

In this study, we screened a large cohort of subjects with increased risk for diabetes and identified two subjects with impaired glucose tolerance carrying common, heterozygous, missense mutations in the PDX1 coding region leading to single amino acid exchanges (P33T, C18R) in its transactivation domain. We generated iPSCs from patients with heterozygous PDX1^P33T/+, PDX1^C18R/+ mutations and engineered isogenic cell lines carrying homozygous PDX1^P33T/P33T, PDX1^C18R/C18R mutations and a heterozygous PDX1 loss-of-function mutation (PDX1^+/-).

RESULTS

Using an in vitro β-cell differentiation protocol, we demonstrated that both, heterozygous PDX1^P33T/+, PDX1^C18R/+ and homozygous PDX1^P33T/P33T, PDX1^C18R/C18R mutations impair β-cell differentiation and function. Furthermore, PDX1^+/- and PDX1^P33T/P33T mutations reduced differentiation efficiency of pancreatic progenitors (PPs), due to downregulation of PDX1-bound genes, including transcription factors MNX1 and PDX1 as well as insulin resistance gene CES1. Additionally, both PDX1^P33T/+ and PDX1^P33T/P33T mutations in PPs reduced the expression of PDX1-bound genes including the long-noncoding RNA, MEG3 and the imprinted gene NNAT, both involved in insulin synthesis and secretion.

CONCLUSIONS

Our results reveal mechanistic details of how common coding mutations in PDX1 impair human pancreatic endocrine lineage formation and β-cell function and contribute to the predisposition for diabetes.

Identification of 28 new susceptibility loci for type 2 diabetes in the Japanese population

To understand the genetics of type 2 diabetes in people of Japanese ancestry, we conducted A meta-analysis of four genome-wide association studies (GWAS; 36,614 cases and 155,150 controls of Japanese ancestry). We identified 88 type 2 diabetes-associated loci (P < 5.0 × 10^-8) with 115 independent signals (P < 5.0 × 10^-6), of which 28 loci with 30 signals were novel. Twenty-eight missense variants were in linkage disequilibrium (r² > 0.6) with the lead variants. Among the 28 missense variants, three previously unreported variants had distinct minor allele frequency (MAF) spectra between people of Japanese and European ancestry (MAF_JPN > 0.05 versus MAF_EUR < 0.01), including missense variants in genes related to pancreatic acinar cells (GP2) and insulin secretion (GLP1R). Transethnic comparisons of the molecular pathways identified from the GWAS results highlight both ethnically shared and heterogeneous effects of a series of pathways on type 2 diabetes (for example, monogenic diabetes and beta cells).