The expression of genes in top obesity-associated loci is enriched in insula and substantia nigra brain regions involved in addiction and reward.

BACKGROUND: Genome-wide association studies (GWAS) have identified more than 250 loci associated with body mass index (BMI) and obesity. However, post-GWAS functional genomic investigations have been inadequate for understanding how these genetic loci physiologically impact disease development. METHODS: We performed a PCR-free expression assay targeting genes located nearby the GWAS-identified SNPs associated with BMI/obesity in a large panel of human tissues. Furthermore, we analyzed several genetic risk scores (GRS) summing GWAS-identified alleles associated with increased BMI in 4236 individuals. RESULTS: We found that the expression of BMI/obesity susceptibility genes was strongly enriched in the brain, especially in the insula (p = 4.7 × 10-9) and substantia nigra (p = 6.8 × 10-7), which are two brain regions involved in addiction and reward. Inversely, we found that top obesity/BMI-associated loci, including FTO, showed the strongest gene expression enrichment in the two brain regions. CONCLUSIONS: Our data suggest for the first time that the susceptibility genes for common obesity may have an effect on eating addiction and reward behaviors through their high expression in substantia nigra and insula, i.e., a different pattern from monogenic obesity genes that act in the hypothalamus and cause hyperphagia. Further epidemiological studies with relevant food behavior phenotypes are necessary to confirm these findings.

Loss-of-function mutations in MRAP2 are pathogenic in hyperphagic obesity with hyperglycemia and hypertension.

The G-protein-coupled receptor accessory protein MRAP2 is implicated in energy control in rodents, notably via the melanocortin-4 receptor1. Although some MRAP2 mutations have been described in people with obesity1-3, their functional consequences on adiposity remain elusive. Using large-scale sequencing of MRAP2 in 9,418 people, we identified 23 rare heterozygous variants associated with increased obesity risk in both adults and children. Functional assessment of each variant shows that loss-of-function MRAP2 variants are pathogenic for monogenic hyperphagic obesity, hyperglycemia and hypertension. This contrasts with other monogenic forms of obesity characterized by excessive hunger, including melanocortin-4 receptor deficiency, that present with low blood pressure and normal glucose tolerance4. The pleiotropic metabolic effect of loss-of-function mutations in MRAP2 might be due to the failure of different MRAP2-regulated G-protein-coupled receptors in various tissues including pancreatic islets.

Expression and functional assessment of candidate type 2 diabetes susceptibility genes identify four new genes contributing to human insulin secretion.

OBJECTIVES: Genome-wide association studies (GWAS) have identified >100 loci independently contributing to type 2 diabetes (T2D) risk. However, translational implications for precision medicine and for the development of novel treatments have been disappointing, due to poor knowledge of how these loci impact T2D pathophysiology. Here, we aimed to measure the expression of genes located nearby T2D associated signals and to assess their effect on insulin secretion from pancreatic beta cells. METHODS: The expression of 104 candidate T2D susceptibility genes was measured in a human multi-tissue panel, through PCR-free expression assay. The effects of the knockdown of beta-cell enriched genes were next investigated on insulin secretion from the human EndoC-ßH1 beta-cell line. Finally, we performed RNA-sequencing (RNA-seq) so as to assess the pathways affected by the knockdown of the new genes impacting insulin secretion from EndoC-ßH1, and we analyzed the expression of the new genes in mouse models with altered pancreatic beta-cell function. RESULTS: We found that the candidate T2D susceptibility genes' expression is significantly enriched in pancreatic beta cells obtained by laser capture microdissection or sorted by flow cytometry and in EndoC-ßH1 cells, but not in insulin sensitive tissues. Furthermore, the knockdown of seven T2D-susceptibility genes (CDKN2A, GCK, HNF4A, KCNK16, SLC30A8, TBC1D4, and TCF19) with already known expression and/or function in beta cells changed insulin secretion, supporting our functional approach. We showed first evidence for a role in insulin secretion of four candidate T2D-susceptibility genes (PRC1, SRR, ZFAND3, and ZFAND6) with no previous knowledge of presence and function in beta cells. RNA-seq in EndoC-ßH1 cells with decreased expression of PRC1, SRR, ZFAND6, or ZFAND3 identified specific gene networks related to T2D pathophysiology. Finally, a positive correlation between the expression of Ins2 and the expression of Prc1, Srr, Zfand6, and Zfand3 was found in mouse pancreatic islets with altered beta-cell function. CONCLUSIONS: This study showed the ability of post-GWAS functional studies to identify new genes and pathways involved in human pancreatic beta-cell function and in T2D pathophysiology.

Beneficial effect of a high number of copies of salivary amylase AMY1 gene on obesity risk in Mexican children.

AIMS/HYPOTHESIS: Childhood obesity is a major public health problem in Mexico, affecting one in every three children. Genome-wide association studies identified genetic variants associated with childhood obesity, but a large missing heritability remains to be elucidated. We have recently shown a strong association between a highly polymorphic copy number variant encompassing the salivary amylase gene (AMY1 also known as AMY1A) and obesity in European and Asian adults. In the present study, we aimed to evaluate the association between AMY1 copy number and obesity in Mexican children. METHODS: We evaluated the number of AMY1 copies in 597 Mexican children (293 obese children and 304 normal weight controls) through highly sensitive digital PCR. The effect of AMY1 copy number on obesity status was assessed using a logistic regression model adjusted for age and sex. RESULTS: We identified a marked effect of AMY1 copy number on reduced risk of obesity (OR per estimated copy 0.84, with the number of copies ranging from one to 16 in this population; p = 4.25 × 10(-6)). The global association between AMY1 copy number and reduced risk of obesity seemed to be mostly driven by the contribution of the highest AMY1 copy number. Strikingly, all children with >10 AMY1 copies were normal weight controls. CONCLUSIONS/INTERPRETATION: Salivary amylase initiates the digestion of dietary starch, which is highly consumed in Mexico. Our current study suggests putative benefits of high number of AMY1 copies (and related production of salivary amylase) on energy metabolism in Mexican children.

Low copy number of the salivary amylase gene predisposes to obesity.

Common multi-allelic copy number variants (CNVs) appear enriched for phenotypic associations compared to their biallelic counterparts. Here we investigated the influence of gene dosage effects on adiposity through a CNV association study of gene expression levels in adipose tissue. We identified significant association of a multi-allelic CNV encompassing the salivary amylase gene (AMY1) with body mass index (BMI) and obesity, and we replicated this finding in 6,200 subjects. Increased AMY1 copy number was positively associated with both amylase gene expression (P = 2.31 × 10(-14)) and serum enzyme levels (P < 2.20 × 10(-16)), whereas reduced AMY1 copy number was associated with increased BMI (change in BMI per estimated copy = -0.15 (0.02) kg/m(2); P = 6.93 × 10(-10)) and obesity risk (odds ratio (OR) per estimated copy = 1.19, 95% confidence interval (CI) = 1.13-1.26; P = 1.46 × 10(-10)). The OR value of 1.19 per copy of AMY1 translates into about an eightfold difference in risk of obesity between subjects in the top (copy number > 9) and bottom (copy number < 4) 10% of the copy number distribution. Our study provides a first genetic link between carbohydrate metabolism and BMI and demonstrates the power of integrated genomic approaches beyond genome-wide association studies.

Alternative human liver transcripts of TCF7L2 bind to the gluconeogenesis regulator HNF4α at the protein level.

AIMS/HYPOTHESIS: Gene polymorphisms of TCF7L2 are associated with increased risk of type 2 diabetes and transcription factor 7-like 2 (TCF7L2) plays a role in hepatic glucose metabolism. We therefore addressed the impact of TCF7L2 isoforms on hepatocyte nuclear factor 4α (HNF4α) and the regulation of gluconeogenesis genes. METHODS: Liver TCF7L2 transcripts were analysed by quantitative PCR in 33 non-diabetic and 31 type 2 diabetic obese individuals genotyped for TCF7L2 rs7903146. To analyse transcriptional regulation by TCF7L2, small interfering RNA transfection, luciferase reporter and co-immunoprecipitation assays were performed in human hepatoma HepG2 cells. RESULTS: In livers of diabetic compared with normoglycaemic individuals, five C-terminal TCF7L2 transcripts showed increased expression. The type 2 diabetes risk allele of rs7903146 positively correlated with TCF7L2 expression in livers from normoglycaemic individuals only. In HepG2 cells, transcript and TCF7L2 protein levels were increased upon incubation in high glucose and insulin. Of the exon 13 transcripts, six were increased in a glucose dose-responsive manner. TCF7L2 transcriptionally regulated 29 genes related to glucose metabolism, including glucose-6-phosphatase. In cultured HepG2 cells, TCF7L2 did not regulate HNF4Α and FOXO1 transcription, but did affect HNF4α protein expression. The TCF7L2 isoforms T6 and T8 (without exon 13 and with exon 15/14, respectively) specifically interacted with HNF4α. CONCLUSIONS/INTERPRETATION: The different levels of expression of alternative C-terminal TCF7L2 transcripts in HepG2 cells, in livers of normoglycaemic individuals carrying the rs7901346 type 2 diabetes risk allele and in livers of diabetic individuals suggest that these transcripts play a role in the pathophysiology of type 2 diabetes. We also report for the first time a protein interaction in HepG2 cells between HNF4α and the T6 and T8 isoforms of TCF7L2, which suggests a distinct role for these specific alternative transcripts.

Whole-exome sequencing and high throughput genotyping identified KCNJ11 as the thirteenth MODY gene.

BACKGROUND: Maturity-onset of the young (MODY) is a clinically heterogeneous form of diabetes characterized by an autosomal-dominant mode of inheritance, an onset before the age of 25 years, and a primary defect in the pancreatic beta-cell function. Approximately 30% of MODY families remain genetically unexplained (MODY-X). Here, we aimed to use whole-exome sequencing (WES) in a four-generation MODY-X family to identify a new susceptibility gene for MODY. METHODOLOGY: WES (Agilent-SureSelect capture/Illumina-GAIIx sequencing) was performed in three affected and one non-affected relatives in the MODY-X family. We then performed a high-throughput multiplex genotyping (Illumina-GoldenGate assay) of the putative causal mutations in the whole family and in 406 controls. A linkage analysis was also carried out. PRINCIPAL FINDINGS: By focusing on variants of interest (i.e. gains of stop codon, frameshift, non-synonymous and splice-site variants not reported in dbSNP130) present in the three affected relatives and not present in the control, we found 69 mutations. However, as WES was not uniform between samples, a total of 324 mutations had to be assessed in the whole family and in controls. Only one mutation (p.Glu227Lys in KCNJ11) co-segregated with diabetes in the family (with a LOD-score of 3.68). No KCNJ11 mutation was found in 25 other MODY-X unrelated subjects. CONCLUSIONS/SIGNIFICANCE: Beyond neonatal diabetes mellitus (NDM), KCNJ11 is also a MODY gene ('MODY13'), confirming the wide spectrum of diabetes related phenotypes due to mutations in NDM genes (i.e. KCNJ11, ABCC8 and INS). Therefore, the molecular diagnosis of MODY should include KCNJ11 as affected carriers can be ideally treated with oral sulfonylureas.

Low-frequency variants in HMGA1 are not associated with type 2 diabetes risk.

It has recently been suggested that the low-frequency c.136-14_136-13insC variant in high-mobility group A1 (HMGA1) may strongly contribute to insulin resistance and type 2 diabetes risk. In our study, we attempted to confirm that HMGA1 is a novel type 2 diabetes locus in French Caucasians. The gene was sequenced in 368 type 2 diabetic case subjects with a family history of type 2 diabetes and 372 normoglycemic control subjects without a family history of type 2 diabetes. None of the 41 genetic variations identified were associated with type 2 diabetes. The lack of association between the c.136-14_136-13insC variant and type 2 diabetes was confirmed in an independent French group of 4,538 case subjects and 4,015 control subjects and in a large meta-analysis of 16,605 case subjects and 46,179 control subjects. Finally, this variant had no effects on metabolic traits and was not involved in variations of HMGA1 and insulin receptor (INSR) expressions. The c.136-14_136-13insC variant was not associated with type 2 diabetes in individuals of European descent. Our study emphasizes the need to analyze a large number of subjects to reliably assess the association of low-frequency variants with the disease.