Maternal Fiber Intake During Pregnancy and Development of Attention-Deficit/Hyperactivity Disorder Symptoms Across Childhood: The Norwegian Mother, Father, and Child Cohort Study

BACKGROUND

Epidemiological studies suggest that maternal diet quality during pregnancy may influence the risk of neurodevelopmental disorders in offspring. Here, we investigated associations between maternal intake of dietary fiber and attention-deficit/hyperactivity disorder (ADHD) symptoms in early childhood.

METHODS

We used longitudinal data of up to 21,852 mother-father-child trios (49.2% female offspring) from MoBa (the Norwegian Mother, Father, and Child Cohort Study). The relationships between maternal fiber intake during pregnancy and offspring ADHD symptoms at ages 3, 5, and 8 years were examined using 1) multivariate regression (overall levels of ADHD symptoms), 2) latent class analysis (subclasses of ADHD symptoms by sex at each age), and 3) latent growth curves (longitudinal change in offspring ADHD symptoms). Covariates were ADHD polygenic scores in child and parents, total energy intake and energy-adjusted sugar intake, parental ages at birth of the child, and sociodemographic factors.

RESULTS

Higher maternal prenatal fiber intake was associated with lower offspring ADHD symptom scores at all ages (Bage3 = -0.14 [95% CI, -0.18 to -0.10]; Bage5 = -0.14 [95% CI, -0.19 to -0.09]; Bage8 = -0.14 [95% CI, -0.20 to -0.09]). Of the derived low/middle/high subclasses of ADHD symptoms, fiber was associated with lower risk of belonging to the middle subclass for boys and girls and to the high subclass for girls only (middle: odds ratioboys 0.91 [95% CI, 0.86 to 0.97]/odds ratiogirls 0.86 [95% CI, 0.81 to 0.91]; high: odds ratiogirls 0.82 [95% CI, 0.72 to 0.94]). Maternal fiber intake and rate of change in child ADHD symptoms across ages were not associated.

CONCLUSIONS

Low prenatal maternal fiber intake may increase symptom levels of ADHD in offspring during childhood, independently of genetic predisposition to ADHD, unhealthy dietary exposures, and sociodemographic factors.

Bioinformatics pipeline for the systematic mining genomic and proteomic variation linked to rare diseases: The example of monogenic diabetes

Monogenic diabetes is characterized as a group of diseases caused by rare variants in single genes. Like for other rare diseases, multiple genes have been linked to monogenic diabetes with different measures of pathogenicity, but the information on the genes and variants is not unified among different resources, making it challenging to process them informatically. We have developed an automated pipeline for collecting and harmonizing data on genetic variants linked to monogenic diabetes. Furthermore, we have translated variant genetic sequences into protein sequences accounting for all protein isoforms and their variants. This allows researchers to consolidate information on variant genes and proteins linked to monogenic diabetes and facilitates their study using proteomics or structural biology. Our open and flexible implementation using Jupyter notebooks enables tailoring and modifying the pipeline and its application to other rare diseases.

Intrauterine Growth and Offspring Neurodevelopmental Traits: A Mendelian Randomization Analysis of the Norwegian Mother, Father and Child Cohort Study (MoBa)

Importance

Conventional epidemiological analyses have suggested that lower birth weight is associated with later neurodevelopmental difficulties; however, it is unclear whether this association is causal.

Objective

To investigate the relationship between intrauterine growth and offspring neurodevelopmental difficulties.

Design, Setting, and Participants

MoBa is a population-based pregnancy cohort that recruited pregnant women from June 1999 to December 2008 included approximately 114 500 children, 95 200 mothers, and 75 200 fathers. Observational associations between birth weight and neurodevelopmental difficulties were assessed with a conventional epidemiological approach. Mendelian randomization analyses were performed to investigate the potential causal association between maternal allele scores for birth weight and offspring neurodevelopmental difficulties conditional on offspring allele scores.

Exposures

Birth weight and maternal allele scores for birth weight (derived from genetic variants robustly associated with birth weight) were the exposures in the observational and mendelian randomization analyses, respectively.

Main Outcomes and Measures

Clinically relevant maternal ratings of offspring neurodevelopmental difficulties at 6 months, 18 months, 3 years, 5 years, and 8 years of age assessing language and motor difficulties, inattention and hyperactivity-impulsivity, social communication difficulties, and repetitive behaviors.

Results

The conventional epidemiological sample included up to 46 970 offspring, whereas the mendelian randomization sample included up to 44 134 offspring (median offspring birth year, 2005 [range, 1999-2009]; mean [SD] maternal age at birth, 30.1 [4.5] years; mean [SD] paternal age at birth, 32.5 [5.1] years). The conventional epidemiological analyses found evidence that birth weight was negatively associated with several domains at multiple offspring ages (outcome of autism-related trait scores: Social Communication Questionnaire [SCQ]-full at 3 years, β = -0.046 [95% CI, -0.057 to -0.034]; SCQ-Restricted and Repetitive Behaviors subscale at 3 years, β = -0.049 [95% CI, -0.060 to -0.038]; attention-deficit/hyperactivity disorder [ADHD] trait scores: Child Behavior Checklist [CBCL]-ADHD subscale at 18 months, β = -0.035 [95% CI, -0.045 to -0.024]; CBCL-ADHD at 3 years, β = -0.032 [95% CI, -0.043 to -0.021]; CBCL-ADHD at 5 years, β = -0.050 [95% CI, -0.064 to -0.037]; Rating Scale for Disruptive Behavior Disorders [RS-DBD]-ADHD at 8 years, β = -0.036 [95% CI, -0.049 to -0.023]; RS-DBD-Inattention at 8 years, β = -0.037 [95% CI, -0.050 to -0.024]; RS-DBD-Hyperactive-Impulsive Behavior at 8 years, β = -0.027 [95% CI, -0.040 to -0.014]; Conners Parent Rating Scale-Revised [Short Form] at 5 years, β = -0.041 [95% CI, -0.054 to -0.028]; motor scores: Ages and Stages Questionnaire-Motor Difficulty [ASQ-MOTOR] at 18 months, β = -0.025 [95% CI, -0.035 to -0.015]; ASQ-MOTOR at 3 years, β = -0.029 [95% CI, -0.040 to -0.018]; and Child Development Inventory-Gross and Fine Motor Skills at 5 years, β = -0.028 [95% CI, -0.042 to -0.015]). Mendelian randomization analyses did not find any evidence for an association between maternal allele scores for birth weight and offspring neurodevelopmental difficulties.

Conclusions and Relevance

This study found that the maternal intrauterine environment, as proxied by maternal birth weight genetic variants, is unlikely to be a major determinant of offspring neurodevelopmental outcomes.

Congenital hyperinsulinism

This clinical review will give doctors who work with children and neonates an introduction to the diagnosis and treatment of congenital hyperinsulinism, the most common cause of persistent neonatal hypoglycaemia. The condition is a rare monogenic disorder characterised by elevated insulin secretion and is a result of mutations in genes that regulate insulin secretion from pancreatic beta cells. The anabolic effect of insulin induces systemic glucose uptake and inhibits gluconeogenesis, glycogenolysis, ketogenesis and lipolysis. Low levels of glucose and ketone bodies in the blood are harmful to the central nervous system and can lead to brain damage or death. Early diagnosis and treatment of congenital hyperinsulinism are therefore crucial for a good prognosis.

Structural properties of the HNF-1A transactivation domain

Hepatocyte nuclear factor 1α (HNF-1A) is a transcription factor with important gene regulatory roles in pancreatic β-cells. HNF1A gene variants are associated with a monogenic form of diabetes (HNF1A-MODY) or an increased risk for type 2 diabetes. While several pancreatic target genes of HNF-1A have been described, a lack of knowledge regarding the structure-function relationships in HNF-1A prohibits a detailed understanding of HNF-1A-mediated gene transcription, which is important for precision medicine and improved patient care. Therefore, we aimed to characterize the understudied transactivation domain (TAD) of HNF-1A in vitro. We present a bioinformatic approach to dissect the TAD sequence, analyzing protein structure, sequence composition, sequence conservation, and the existence of protein interaction motifs. Moreover, we developed the first protocol for the recombinant expression and purification of the HNF-1A TAD. Small-angle X-ray scattering and synchrotron radiation circular dichroism suggested a disordered conformation for the TAD. Furthermore, we present functional data on HNF-1A undergoing liquid-liquid phase separation, which is in line with in silico predictions and may be of biological relevance for gene transcriptional processes in pancreatic β-cells.

Extending protein interaction networks using proteoforms and small molecules

MOTIVATION

Biological network analysis for high-throughput biomedical data interpretation relies heavily on topological characteristics. Networks are commonly composed of nodes representing genes or proteins that are connected by edges when interacting. In this study, we use the rich information available in the Reactome pathway database to build biological networks accounting for small molecules and proteoforms modeled using protein isoforms and post-translational modifications to study the topological changes induced by this refinement of the network representation.

RESULTS

We find that improving the interactome modeling increases the number of nodes and interactions, but that isoform and post-translational modification annotation is still limited compared to what can be expected biologically. We also note that small molecule information can distort the topology of the network due to the high connectedness of these molecules, which does not necessarily represent the reality of biology. However, by restricting the connections of small molecules to the context of biochemical reactions, we find that these improve the overall connectedness of the network and reduce the prevalence of isolated components and nodes. Overall, changing the representation of the network alters the prevalence of articulation points and bridges globally but also within and across pathways. Hence, some molecules can gain or lose in biological importance depending on the level of detail of the representation of the biological system, which might in turn impact network-based studies of diseases or druggability.

AVAILABILITY AND IMPLEMENTATION

Networks are constructed based on data publicly available in the Reactome Pathway knowledgebase: reactome.org.

Second international consensus report on gaps and opportunities for the clinical translation of precision diabetes medicine

Precision medicine is part of the logical evolution of contemporary evidence-based medicine that seeks to reduce errors and optimize outcomes when making medical decisions and health recommendations. Diabetes affects hundreds of millions of people worldwide, many of whom will develop life-threatening complications and die prematurely. Precision medicine can potentially address this enormous problem by accounting for heterogeneity in the etiology, clinical presentation and pathogenesis of common forms of diabetes and risks of complications. This second international consensus report on precision diabetes medicine summarizes the findings from a systematic evidence review across the key pillars of precision medicine (prevention, diagnosis, treatment, prognosis) in four recognized forms of diabetes (monogenic, gestational, type 1, type 2). These reviews address key questions about the translation of precision medicine research into practice. Although not complete, owing to the vast literature on this topic, they revealed opportunities for the immediate or near-term clinical implementation of precision diabetes medicine; furthermore, we expose important gaps in knowledge, focusing on the need to obtain new clinically relevant evidence. Gaps include the need for common standards for clinical readiness, including consideration of cost-effectiveness, health equity, predictive accuracy, liability and accessibility. Key milestones are outlined for the broad clinical implementation of precision diabetes medicine.

The Composite Autonomic Symptom Score 31 Questionnaire: A Sensitive Test to Detect Risk for Autonomic Neuropathy

Background and Aims

Autonomic neuropathy is a common but often neglected complication of diabetes, prediabetes, and even in individuals with an elevated risk of diabetes. The Composite Autonomic Symptom Score (COMPASS) 31 is a validated and easy-to-use questionnaire regarding autonomic symptoms. We aimed to use a digitally, Norwegian version of the COMPASS 31 in people with different durations of diabetes and healthy controls to consider feasibility and to investigate if scores could discriminate between positive and negative outcomes for established tests for diabetic neuropathy, including cardiovascular autonomic neuropathy (CAN) and a novel method of examining the gastrointestinal visceral sensitivity.

Method

We included 21 participants with longstanding type 2 diabetes, 15 with early type 2 diabetes, and 30 matched controls. The mean age for all groups was 69 years. Participants were phenotyped by cardiovascular autonomic reflex tests, electrical skin conductance, sural nerve electrophysiology, and the monofilament test. As a proxy for gastrointestinal visceral and autonomic nerve function, evoked potentials were measured following rapid rectal balloon distention.

Results

Participants with longstanding diabetes scored a median (IQR) of 14.9 (10.8-28.7) points, early diabetes of 7.3 (1.6-15.2), and matched controls of 8.6 (4.1-21.6), p = 0.04. Women and men scored 14.4 (5.5-28.7) and 7.8 (3.6-14.6) points, respectively, p = 0.01. Participants with definite or borderline CAN scored 14.3 (10.4-31.9) points, compared to participants with no CAN, 8.3 (3.2-21.5), p = 0.04. Lowering the diagnostic cut-off from 16 to 10 points increased the sensitivity from 0.33 to 0.83, with a decreased specificity from 0.68 to 0.55.

Conclusion

We successfully used COMPASS 31 in Norwegian. Thus, following the guidelines, we suggest clinical implementation for the assessment of autonomic neuropathy. Participants with longstanding diabetes had an increased likelihood of symptoms and signs of autonomic neuropathy. For screening purposes, the sensitivity was improved by lowering the cut-off to 10 points, with a lower score nearly excluding the diagnosis.

Developmental milestones in early childhood and genetic liability to neurodevelopmental disorders

BACKGROUND

Timing of developmental milestones, such as age at first walking, is associated with later diagnoses of neurodevelopmental disorders. However, its relationship to genetic risk for neurodevelopmental disorders in the general population is unknown. Here, we investigate associations between attainment of early-life language and motor development milestones and genetic liability to autism, attention deficit hyperactivity disorder (ADHD), and schizophrenia.

METHODS

We use data from a genotyped sub-set (N = 25699) of children in the Norwegian Mother, Father and Child Cohort Study (MoBa). We calculate polygenic scores (PGS) for autism, ADHD, and schizophrenia and predict maternal reports of children's age at first walking, first words, and first sentences, motor delays (18 months), and language delays and a generalised measure of concerns about development (3 years). We use linear and probit regression models in a multi-group framework to test for sex differences.

RESULTS

We found that ADHD PGS were associated with earlier walking age (β = -0.033, padj < 0.001) in both males and females. Additionally, autism PGS were associated with later walking (β = 0.039, padj = 0.006) in females only. No robust associations were observed for schizophrenia PGS or between any neurodevelopmental PGS and measures of language developmental milestone attainment.

CONCLUSIONS

Genetic liabilities for neurodevelopmental disorders show some specific associations with the age at which children first walk unsupported. Associations are small but robust and, in the case of autism PGS, differentiated by sex. These findings suggest that early-life motor developmental milestone attainment is associated with genetic liability to ADHD and autism in the general population.

Early manifestations of genetic risk for neurodevelopmental disorders

BACKGROUND

Attention deficit/hyperactivity disorder (ADHD), autism spectrum disorder (autism) and schizophrenia are highly heritable neurodevelopmental disorders, affecting the lives of many individuals. It is important to increase our understanding of how the polygenic risk for neurodevelopmental disorders manifests during childhood in boys and girls.

METHODS

Polygenic risk scores (PRS) for ADHD, autism and schizophrenia were calculated in a subsample of 15 205 children from the Norwegian Mother, Father and Child Cohort Study (MoBa). Mother-reported traits of repetitive behavior, social communication, language and motor difficulties, hyperactivity and inattention were measured in children at 6 and 18 months, 3, 5 and 8 years. Linear regression models in a multigroup framework were used to investigate associations between the three PRS and dimensional trait measures in MoBa, using sex as a grouping variable.

RESULTS

Before the age of 2, the ADHD PRS was robustly associated with hyperactivity and inattention, with increasing strength up to 8 years, and with language difficulties at age 5 and 8. The autism PRS was robustly associated with language difficulties at 18 months, motor difficulties at 36 months, and hyperactivity and inattention at 8 years. We did not identify robust associations for the schizophrenia PRS. In general, the PRS associations were similar in boys and girls. The association between ADHD PRS and hyperactivity at 18 months was, however, stronger in boys.

CONCLUSIONS

Polygenic risk for autism and ADHD in the general population manifests early in childhood and broadly across behavioral measures of neurodevelopmental traits.

Within-sibship genome-wide association analyses decrease bias in estimates of direct genetic effects

Estimates from genome-wide association studies (GWAS) of unrelated individuals capture effects of inherited variation (direct effects), demography (population stratification, assortative mating) and relatives (indirect genetic effects). Family-based GWAS designs can control for demographic and indirect genetic effects, but large-scale family datasets have been lacking. We combined data from 178,086 siblings from 19 cohorts to generate population (between-family) and within-sibship (within-family) GWAS estimates for 25 phenotypes. Within-sibship GWAS estimates were smaller than population estimates for height, educational attainment, age at first birth, number of children, cognitive ability, depressive symptoms and smoking. Some differences were observed in downstream SNP heritability, genetic correlations and Mendelian randomization analyses. For example, the within-sibship genetic correlation between educational attainment and body mass index attenuated towards zero. In contrast, analyses of most molecular phenotypes (for example, low-density lipoprotein-cholesterol) were generally consistent. We also found within-sibship evidence of polygenic adaptation on taller height. Here, we illustrate the importance of family-based GWAS data for phenotypes influenced by demographic and indirect genetic effects.

Structural and biophysical characterization of transcription factor HNF-1A as a tool to study MODY3 diabetes variants

Hepatocyte nuclear factor 1A (HNF-1A) is a transcription factor expressed in several embryonic and adult tissues, modulating the expression of numerous target genes. Pathogenic variants in the HNF1A gene are known to cause maturity-onset diabetes of the young 3 (MODY3 or HNF1A MODY), a disease characterized by dominant inheritance, age of onset before 25 to 35 years of age, and pancreatic β-cell dysfunction. A precise diagnosis can alter management of this disease, as insulin can be exchanged with sulfonylurea tablets and genetic counseling differs from polygenic forms of diabetes. Therefore, more knowledge on the mechanisms of HNF-1A function and the level of pathogenicity of the numerous HNF1A variants is required for precise diagnostics. Here, we structurally and biophysically characterized an HNF-1A protein containing both the DNA-binding domain and the dimerization domain, and determined the folding and DNA-binding capacity of two established MODY3 HNF-1A variant proteins (P112L, R263C) and one variant of unknown significance (N266S). All three variants showed reduced functionality compared to the WT protein. Furthermore, while the R263C and N266S variants displayed reduced binding to an HNF-1A target promoter, we found the P112L variant was unstable in vitro and in cells. Our results support and mechanistically explain disease causality for these investigated variants and present a novel approach for the dissection of structurally unstable and DNA-binding defective variants. This study indicates that structural and biochemical investigation of HNF-1A is a valuable tool in reliable variant classification needed for precision diabetes diagnostics and management.

Fucosylated AGP glycopeptides as biomarkers of HNF1A-Maturity onset diabetes of the young

AIMS

We previously demonstrated that antennary fucosylated N-glycans on plasma proteins are regulated by HNF1A and can identify cases of Maturity-Onset Diabetes of the Young caused by HNF1A variants (HNF1A-MODY). Based on literature data, we further postulated that N-glycans with best diagnostic value mostly originate from alpha-1-acid glycoprotein (AGP). In this study we analyzed fucosylation of AGP in subjects with HNF1A-MODY and other types of diabetes aiming to evaluate its diagnostic potential.

METHODS

A recently developed LC-MS method for AGP N-glycopeptide analysis was utilized in two independent cohorts: a) 466 subjects with different diabetes subtypes to test the fucosylation differences, b) 98 selected individuals to test the discriminative potential for pathogenic HNF1A variants.

RESULTS

Our results showed significant reduction in AGP fucosylation associated to HNF1A-MODY when compared to other diabetes subtypes. Additionally, ROC curve analysis confirmed significant discriminatory potential of individual fucosylated AGP glycopeptides, where the best performing glycopeptide had an AUC of 0.94 (95% CI 0.90-0.99).

CONCLUSIONS

A glycopeptide based diagnostic tool would be beneficial for patient stratification by providing information about the functionality of HNF1A. It could assist the interpretation of DNA sequencing results and be a useful addition to the differential diagnostic process.

Subtype-Specific Surface Proteins on Adipose Tissue Macrophages and Their Association to Obesity-Induced Insulin Resistance

A chronic low-grade inflammation, originating in the adipose tissue, is considered a driver of obesity-associated insulin resistance. Macrophage composition in white adipose tissue is believed to contribute to the pathogenesis of metabolic diseases, but a detailed characterization of pro- and anti-inflammatory adipose tissue macrophages (ATMs) in human obesity and how they are distributed in visceral- and subcutaneous adipose depots is lacking. In this study, we performed a surface proteome screening of pro- and anti-inflammatory ATMs in both subcutaneous- (SAT) and visceral adipose tissue (VAT) and evaluated their relationship with systemic insulin resistance. From the proteomics screen we found novel surface proteins specific to M1-like- and M2-like macrophages, and we identified depot-specific immunophenotypes in SAT and VAT. Furthermore, we found that insulin resistance, assessed by HOMA-IR, was positively associated with a relative increase in pro-inflammatory M1-like macrophages in both SAT and VAT.

Maternal and offspring genetic risk score analyses of fetal alcohol exposure and attention-deficit hyperactivity disorder risk in offspring

Background

Studies investigating the effects of prenatal alcohol exposure on childhood attention‐deficit hyperactivity disorder (ADHD) symptoms using conventional observational designs have reported inconsistent findings, which may be affected by unmeasured confounding and maternal and fetal ability to metabolize alcohol. We used genetic variants from the alcohol metabolizing genes, alcohol dehydrogenase (ADH) and aldehyde dehydrogenase (ALDH), as proxies for fetal alcohol exposure to investigate their association with risk of offspring ADHD symptoms around age 7–8 years.

Methods

We used data from 3 longitudinal pregnancy cohorts: Avon Longitudinal Study of Parents and Children (ALSPAC), Generation R study (GenR), and the Norwegian Mother, Father and Child Cohort study (MoBa). Genetic risk scores (GRS) for alcohol use and metabolism using 36 single nucleotide polymorphisms (SNPs) from ADH and ALDH genes were calculated for mothers (N_ALSPAC = 8196; N_MOBA = 13,614), fathers (N_MOBA = 13,935), and offspring (N_ALSPAC=8,237; N_MOBA=14,112; N_GENR=2,661). Associations between maternal GRS and offspring risk of ADHD symptoms were tested in the full sample to avoid collider bias. Offspring GRS analyses were stratified by maternal drinking status.

Results

The pooled estimate in maternal GRS analyses adjusted for offspring GRS in ALSPAC and MoBa was OR = 0.99, 95%CI 0.97–1.02. The pooled estimate in offspring GRS analyses stratified by maternal drinking status across all the cohorts was as follows: OR_DRINKING = 0.98, 95% CI 0.94–1.02; OR_{NO DRINKING} = 0.99, 95% CI 0.97–1.02. These findings remained similar after accounting for maternal genotype data in ALSPAC and maternal and paternal genotype data in MoBa.

Conclusions

We did not find evidence for a causal effect of fetal alcohol exposure on risk of ADHD symptoms in offspring. The results may be affected by limited power to detect small effects and outcome assessment.

Genetic Liability for Schizophrenia and Childhood Psychopathology in the General Population

Genetic liability for schizophrenia is associated with psychopathology in early life. It is not clear if these associations are time dependent during childhood, nor if they are specific across different forms of psychopathology. Using genotype and questionnaire data on children (N = 15 105) from the Norwegian Mother, Father and Child Cohort Study, we used schizophrenia polygenic risk scores to test developmental stability in associations with measures of emotional and behavioral problems between 18 months and 5 years, and domain specificity in associations with symptoms of depression, anxiety, conduct problems, oppositionality, inattention, and hyperactivity at 8 years. We then sought to identify symptom profiles-across development and domains-associated with schizophrenia polygenic liability. We found evidence for developmental stability in associations between schizophrenia polygenic risk scores and emotional and behavioral problems, with the latter being mediated specifically via the rate of change in symptoms (β slope = 0.032; 95% CI: 0.007-0.057). At age 8, associations were better explained by a model of symptom-specific polygenic effects rather than effects mediated via a general psychopathology factor or by domain-specific factors. Overall, individuals with higher schizophrenia polygenic risk scores were more likely (OR = 1.310 [95% CIs: 1.122-1.528]) to have a profile of increasing behavioral and emotional symptoms in early childhood, followed by elevated symptoms of conduct disorder, oppositionality, hyperactivity, and inattention by age 8. Schizophrenia-associated alleles are linked to specific patterns of early-life psychopathology. The associations are small, but findings of this nature can help us better understand the developmental emergence of schizophrenia.

Unsupervised Clustering of Missense Variants in HNF1A Using Multidimensional Functional Data Aids Clinical Interpretation

Exome sequencing in diabetes presents a diagnostic challenge because depending on frequency, functional impact, and genomic and environmental contexts, HNF1A variants can cause maturity-onset diabetes of the young (MODY), increase type 2 diabetes risk, or be benign. A correct diagnosis matters as it informs on treatment, progression, and family risk. We describe a multi-dimensional functional dataset of 73 HNF1A missense variants identified in exomes of 12,940 individuals. Our aim was to develop an analytical framework for stratifying variants along the HNF1A phenotypic continuum to facilitate diagnostic interpretation. HNF1A variant function was determined by four different molecular assays. Structure of the multi-dimensional dataset was explored using principal component analysis, k-means, and hierarchical clustering. Weights for tissue-specific isoform expression and functional domain were integrated. Functionally annotated variant subgroups were used to re-evaluate genetic diagnoses in national MODY diagnostic registries. HNF1A variants demonstrated a range of behaviors across the assays. The structure of the multi-parametric data was shaped primarily by transactivation. Using unsupervised learning methods, we obtained high-resolution functional clusters of the variants that separated known causal MODY variants from benign and type 2 diabetes risk variants and led to reclassification of 4% and 9% of HNF1A variants identified in the UK and Norway MODY diagnostic registries, respectively. Our proof-of-principle analyses facilitated informative stratification of HNF1A variants along the continuum, allowing improved evaluation of clinical significance, management, and precision medicine in diabetes clinics. Transcriptional activity appears a superior readout supporting pursuit of transactivation-centric experimental designs for high-throughput functional screens.

PathwayMatcher: proteoform-centric network construction enables fine-granularity multiomics pathway mapping

BACKGROUND

Mapping biomedical data to functional knowledge is an essential task in bioinformatics and can be achieved by querying identifiers (e.g., gene sets) in pathway knowledge bases. However, the isoform and posttranslational modification states of proteins are lost when converting input and pathways into gene-centric lists.

FINDINGS

Based on the Reactome knowledge base, we built a network of protein-protein interactions accounting for the documented isoform and modification statuses of proteins. We then implemented a command line application called PathwayMatcher (github.com/PathwayAnalysisPlatform/PathwayMatcher) to query this network. PathwayMatcher supports multiple types of omics data as input and outputs the possibly affected biochemical reactions, subnetworks, and pathways.

CONCLUSIONS

PathwayMatcher enables refining the network representation of pathways by including proteoforms defined as protein isoforms with posttranslational modifications. The specificity of pathway analyses is hence adapted to different levels of granularity, and it becomes possible to distinguish interactions between different forms of the same protein.

Sonographic pancreas echogenicity in cystic fibrosis compared to exocrine pancreatic function and pancreas fat content at Dixon-MRI

OBJECTIVE

Fatty infiltration of the pancreas is a dominating feature in cystic fibrosis (CF). We evaluate the association between pancreatic fat content assessed by Dixon magnetic resonance imaging (MRI), pancreatic echogenicity at ultrasonography (US) and exocrine function in CF patients and healthy controls (HC).

MATERIAL AND METHODS

Transabdominal US, pancreatic Dixon-MRI and diffusion-weighted imaging (DWI) were performed in 21 CF patients and 15 HCs. Exocrine function was assessed by endoscopic secretin test and fecal elastase.

RESULTS

CF patients were grouped according to exocrine pancreatic function as subjects with normal (CFS: n = 11) or reduced (CFI: n = 10) function. Among CFI 90% (9/10) had visual hyperechogenicity. CFI also had increased echo-level values (p<0.05 vs others). All CFI (10/10) had markedly increased pancreatic fat content estimated by MRI compared to sufficient groups, p<0.001). Among CFS patients and HC, 27% (3/11) and 33% (5/15), respectively, had hyperechoic pancreas. However, all these had low pancreatic fat-content at MRI compared to CFI. In CFI, pancreatic fat content was correlated to ADC (r = -0.93, p<0.001).

CONCLUSION

Pancreas insufficient CF patients exhibit severe pancreatic fatty-infiltration at MRI and hyperechoic pancreas at US. Pancreas hyperechogenicity in pancreatic sufficient subjects does not co-exist with fatty infiltration at MRI. MRI evaluates pancreatic fatty infiltration more accurately than US and fat infiltration estimated by MRI outperforms sonographic hyper-echogenicity as a marker for exocrine pancreatic failure in CF.

Erratum: Sequence data and association statistics from 12,940 type 2 diabetes cases and controls

This corrects the article DOI: 10.1038/sdata.2017.179.

Sequence data and association statistics from 12,940 type 2 diabetes cases and controls

To investigate the genetic basis of type 2 diabetes (T2D) to high resolution, the GoT2D and T2D-GENES consortia catalogued variation from whole-genome sequencing of 2,657 European individuals and exome sequencing of 12,940 individuals of multiple ancestries. Over 27M SNPs, indels, and structural variants were identified, including 99% of low-frequency (minor allele frequency [MAF] 0.1-5%) non-coding variants in the whole-genome sequenced individuals and 99.7% of low-frequency coding variants in the whole-exome sequenced individuals. Each variant was tested for association with T2D in the sequenced individuals, and, to increase power, most were tested in larger numbers of individuals (>80% of low-frequency coding variants in ~82 K Europeans via the exome chip, and ~90% of low-frequency non-coding variants in ~44 K Europeans via genotype imputation). The variants, genotypes, and association statistics from these analyses provide the largest reference to date of human genetic information relevant to T2D, for use in activities such as T2D-focused genotype imputation, functional characterization of variants or genes, and other novel analyses to detect associations between sequence variation and T2D.

The HNF1A mutant Ala180Val: Clinical challenges in determining causality of a rare HNF1A variant in familial diabetes

AIMS

Heterozygous mutations in hepatocyte nuclear factor-1A (HNF1A) cause maturity-onset diabetes of the young type 3 (MODY3). Our aim was to compare two families with suspected dominantly inherited diabetes and a new HNF1A variant of unknown clinical significance.

METHODS

The HNF1A gene was sequenced in two independently recruited families from the Norwegian MODY Registry. Both familes were phenotyped clinically and biochemically. Microsatellite markers around and within the HNF1A locus were used for haplotyping. Chromosomal linkage analysis was performed in one family, and whole-exome sequencing was undertaken in two affected family members from each family. Transactivation activity, DNA binding and nuclear localization of wild type and mutant HNF-1A were assessed.

RESULTS

The novel HNF1A variant c.539C>T (p.Ala180Val) was found in both families. The variant fully co-segregated with diabetes in one family. In the other family, two subjects with diabetes mellitus and one with normal glucose levels were homozygous variant carriers. Chromosomal linkage of diabetes to the HNF1A locus or to other genomic regions could not be established. The protein functional studies did not reveal significant differences between wild type and variant HNF-1A. In each family, whole-exome sequencing failed to identify any other variant that could explain the disease.

CONCLUSIONS

The HNF1A variant p.Ala180Val does not seem to cause MODY3, although it may confer risk for type 2 diabetes mellitus. Our data demonstrate challenges in causality evaluation of rare variants detected in known diabetes genes.

Nuclear import of glucokinase in pancreatic beta-cells is mediated by a nuclear localization signal and modulated by SUMOylation

The localization of glucokinase in pancreatic beta-cell nuclei is a controversial issue. Although previous reports suggest such a localization, the mechanism for its import has so far not been identified. Using immunofluorescence, subcellular fractionation and mass spectrometry, we present evidence in support of glucokinase localization in beta-cell nuclei of human and mouse pancreatic sections, as well as in human and mouse isolated islets, and murine MIN6 cells. We have identified a conserved, seven-residue nuclear localization signal (³⁰LKKVMRR³⁶) in the human enzyme. Substituting the residues KK^31,32 and RR^35,36 with AA led to a loss of its nuclear localization in transfected cells. Furthermore, our data indicates that SUMOylation of glucokinase modulates its nuclear import, while high glucose concentrations do not significantly alter the enzyme nuclear/cytosolic ratio. Thus, for the first time, we provide data in support of a nuclear import of glucokinase mediated by a redundant mechanism, involving a nuclear localization signal, and which is modulated by its SUMOylation. These findings add new knowledge to the functional role of glucokinase in the pancreatic beta-cell.

Functional Investigations of HNF1A Identify Rare Variants as Risk Factors for Type 2 Diabetes in the General Population

Variants in HNF1A encoding hepatocyte nuclear factor 1α (HNF-1A) are associated with maturity-onset diabetes of the young form 3 (MODY 3) and type 2 diabetes. We investigated whether functional classification of HNF1A rare coding variants can inform models of diabetes risk prediction in the general population by analyzing the effect of 27 HNF1A variants identified in well-phenotyped populations (n = 4,115). Bioinformatics tools classified 11 variants as likely pathogenic and showed no association with diabetes risk (combined minor allele frequency [MAF] 0.22%; odds ratio [OR] 2.02; 95% CI 0.73-5.60; P = 0.18). However, a different set of 11 variants that reduced HNF-1A transcriptional activity to <60% of normal (wild-type) activity was strongly associated with diabetes in the general population (combined MAF 0.22%; OR 5.04; 95% CI 1.99-12.80; P = 0.0007). Our functional investigations indicate that 0.44% of the population carry HNF1A variants that result in a substantially increased risk for developing diabetes. These results suggest that functional characterization of variants within MODY genes may overcome the limitations of bioinformatics tools for the purposes of presymptomatic diabetes risk prediction in the general population.

Causes of death in childhood-onset Type 1 diabetes: long-term follow-up

AIMS

To assess the causes of death and cause-specific standardized mortality ratios in two nationwide, population-based cohorts diagnosed with Type 1 diabetes during the periods 1973-1982 and 1989-2012, and to evaluate changes in causes of death during the follow-up period.

METHODS

People with Type 1 diabetes who were aged < 15 years at diagnosis were identified in the Norwegian Childhood Diabetes Registry and followed from diagnosis until death, emigration or September 2013 (n = 7871). We assessed causes of death by linking data to the nationwide Cause of Death Registry and through a review committee that evaluated medical records, autopsy reports and death certificates.

RESULTS

During a mean (range) follow-up of 16.8 (0-40.7) years, 241 individuals (3.1%) died, representing 132 143 person-years. The leading cause of death before the age of 30 years was acute complications (41/119, 34.5%). After the age of 30 years cardiovascular disease was predominant (41/122, 33.6%), although death attributable to acute complications was still important in this age group (22/122, 18.0%). A total of 5% of deaths were caused by 'dead-in-bed' syndrome. The standardized mortality ratio was elevated for cardiovascular disease [11.9 (95% CI 8.6-16.4)] and violent death [1.7 (95% CI 1.3-2.1)] in both sexes combined, but was elevated for suicide only in women [2.5 (95% CI 1.2-5.3)]. The risk of death from acute complications was approximately half in women compared with men [hazard ratio 0.43 (95% CI 0.25-0.76)], and did not change with more recent year of diagnosis [hazard ratio 1.02 (0.98-1.05)].

CONCLUSIONS

There was no change in mortality attributable to acute complications during the study period. To reduce premature mortality in people with childhood-onset diabetes focus should be on prevention of acute complications. Male gender implied increased risk.

Length of Variable Numbers of Tandem Repeats in the Carboxyl Ester Lipase (CEL) Gene May Confer Susceptibility to Alcoholic Liver Cirrhosis but Not Alcoholic Chronic Pancreatitis

BACKGROUND

Carboxyl-ester lipase (CEL) contributes to fatty acid ethyl ester metabolism, which is implicated in alcoholic pancreatitis. The CEL gene harbours a variable number of tandem repeats (VNTR) region in exon 11. Variation in this VNTR has been linked to monogenic pancreatic disease, while conflicting results were reported for chronic pancreatitis (CP). Here, we aimed to investigate a potential association of CEL VNTR lengths with alcoholic CP.

METHODS

Overall, 395 alcoholic CP patients, 218 patients with alcoholic liver cirrhosis (ALC) serving as controls with a comparable amount of alcohol consumed, and 327 healthy controls from Germany and the United Kingdom (UK) were analysed by determination of fragment lengths by capillary electrophoresis. Allele frequencies and genotypes of different VNTR categories were compared between the groups.

RESULTS

Twelve repeats were overrepresented in UK ACP patients (P = 0.04) compared to controls, whereas twelve repeats were enriched in German ALC compared to alcoholic CP patients (P = 0.03). Frequencies of CEL VNTR lengths of 14 and 15 repeats differed between German ALC patients and healthy controls (P = 0.03 and 0.008, respectively). However, in the genotype and pooled analysis of VNTR lengths no statistical significant association was depicted. Additionally, the 16-16 genotype as well as 16 repeats were more frequent in UK ALC than in alcoholic CP patients (P = 0.034 and 0.02, respectively). In all other calculations, including pooled German and UK data, allele frequencies and genotype distributions did not differ significantly between patients and controls or between alcoholic CP and ALC.

CONCLUSIONS

We did not obtain evidence that CEL VNTR lengths are associated with alcoholic CP. However, our results suggest that CEL VNTR lengths might associate with ALC, a finding that needs to be clarified in larger cohorts.

Structure-function studies of HNF1A (MODY3) gene mutations in South Indian patients with monogenic diabetes

Maturity-onset diabetes of the young (MODY) is a genetically heterogeneous monogenic form of diabetes characterized by onset of diabetes below 25 years of age, autosomal dominant mode of inheritance and primary defect in insulin secretion. Mutations in the gene (HNF1A) encoding transcription factor hepatocyte nuclear factor 1A (HNF-1A) results in one of the most common forms of MODY (MODY3). HNF-1A is mainly enriched in pancreatic β-cells and hepatocytes and important for organ development and normal pancreatic function. We here report on the functional interrogation of eight missense HNF1A mutations associated with MODY3 in South Indian subjects, and the contributing effect of common variant (S487N) within HNF1A. Of the eight mutations, three mutations (p.R171G, p.G245R and p.R263H), in particular, affected HNF-1A function in transfected HeLa cells by reducing both transcriptional activity and nuclear localization, possibly due to disruption of the integrity of the three dimensional structure. The common variant p.S487N contributed further to the loss-of-function of p.R271Q (p.R271Q+p.S487N double mutant), in vitro, on both activity and localization. Our data on the first functional study of HNF1A mutations in South India subjects confers that the defect of the HNF-1A mutant proteins are responsible for MODY3 diabetes in these patients.

Clinical reappraisal of SHORT syndrome with PIK3R1 mutations: toward recommendation for molecular testing and management

SHORT syndrome has historically been defined by its acronym: short stature (S), hyperextensibility of joints and/or inguinal hernia (H), ocular depression (O), Rieger abnormality (R) and teething delay (T). More recently several research groups have identified PIK3R1 mutations as responsible for SHORT syndrome. Knowledge of the molecular etiology of SHORT syndrome has permitted a reassessment of the clinical phenotype. The detailed phenotypes of 32 individuals with SHORT syndrome and PIK3R1 mutation, including eight newly ascertained individuals, were studied to fully define the syndrome and the indications for PIK3R1 testing. The major features described in the SHORT acronym were not universally seen and only half (52%) had four or more of the classic features. The commonly observed clinical features of SHORT syndrome seen in the cohort included intrauterine growth restriction (IUGR) <10th percentile, postnatal growth restriction, lipoatrophy and the characteristic facial gestalt. Anterior chamber defects and insulin resistance or diabetes were also observed but were not as prevalent. The less specific, or minor features of SHORT syndrome include teething delay, thin wrinkled skin, speech delay, sensorineural deafness, hyperextensibility of joints and inguinal hernia. Given the high risk of diabetes mellitus, regular monitoring of glucose metabolism is warranted. An echocardiogram, ophthalmological and hearing assessments are also recommended.

Assessing the Causal Relationship of Maternal Height on Birth Size and Gestational Age at Birth: A Mendelian Randomization Analysis

BACKGROUND

Observational epidemiological studies indicate that maternal height is associated with gestational age at birth and fetal growth measures (i.e., shorter mothers deliver infants at earlier gestational ages with lower birth weight and birth length). Different mechanisms have been postulated to explain these associations. This study aimed to investigate the casual relationships behind the strong association of maternal height with fetal growth measures (i.e., birth length and birth weight) and gestational age by a Mendelian randomization approach.

METHODS AND FINDINGS

We conducted a Mendelian randomization analysis using phenotype and genome-wide single nucleotide polymorphism (SNP) data of 3,485 mother/infant pairs from birth cohorts collected from three Nordic countries (Finland, Denmark, and Norway). We constructed a genetic score based on 697 SNPs known to be associated with adult height to index maternal height. To avoid confounding due to genetic sharing between mother and infant, we inferred parental transmission of the height-associated SNPs and utilized the haplotype genetic score derived from nontransmitted alleles as a valid genetic instrument for maternal height. In observational analysis, maternal height was significantly associated with birth length (p = 6.31 × 10-9), birth weight (p = 2.19 × 10-15), and gestational age (p = 1.51 × 10-7). Our parental-specific haplotype score association analysis revealed that birth length and birth weight were significantly associated with the maternal transmitted haplotype score as well as the paternal transmitted haplotype score. Their association with the maternal nontransmitted haplotype score was far less significant, indicating a major fetal genetic influence on these fetal growth measures. In contrast, gestational age was significantly associated with the nontransmitted haplotype score (p = 0.0424) and demonstrated a significant (p = 0.0234) causal effect of every 1 cm increase in maternal height resulting in ~0.4 more gestational d. Limitations of this study include potential influences in causal inference by biological pleiotropy, assortative mating, and the nonrandom sampling of study subjects.

CONCLUSIONS

Our results demonstrate that the observed association between maternal height and fetal growth measures (i.e., birth length and birth weight) is mainly defined by fetal genetics. In contrast, the association between maternal height and gestational age is more likely to be causal. In addition, our approach that utilizes the genetic score derived from the nontransmitted maternal haplotype as a genetic instrument is a novel extension to the Mendelian randomization methodology in casual inference between parental phenotype (or exposure) and outcomes in offspring.

An Egyptian case of congenital hyperinsulinism of infancy due to a novel mutation in KCNJ11 encoding Kir6.2 and response to octreotide

Congenital hyperinsulinism of infancy (CHI) is a rare heterogeneous disease mostly attributable to mutations in the genes encoding the KATP channel subunits found in pancreatic β-cells. Here, we report a child presenting at day 1 with persistent hyperinsulinemic hypoglycemia and who underwent open laparotomy and subtotal pancreatectomy with resection of tail and body of pancreas at 30 days of age. Normoglycemia was restored by Octreotide that was discontinued when the child was 7-month old. However, 3 months later Octreotide was re-administered as hypoglycemic attacks recurred. On follow-up, the child has adequate glycemic control and is thriving well with no neurodevelopmental morbidity. Genetic analysis revealed the novel mutation c.407G > A [p.R136H] in KCNJ11 encoding Kir6.2, confirming the diffuse form of CHI. This is to our knowledge the first reported Egyptian case of CHI due to a mutation in KCNJ11.

An Egyptian family with H syndrome due to a novel mutation in SLC29A3 illustrating overlapping features with pigmented hypertrichotic dermatosis with insulin-dependent diabetes and Faisalabad histiocytosis

The SLC29A3 gene, encoding hENT3, a member of the equilibrative nucleoside transporter family, has recently been found mutated in Faisalabad histiocytosis, pigmented hypertrichotic dermatosis with insulin-dependent diabetes, familial sinus histiocytosis with massive lymphadenopathy (SHML), and H syndromes. We here report clinical and genetic findings of an Egyptian family with H syndrome. We describe two siblings, a 19-yr old girl and a 15-yr old boy, of consanguineous parents. From 5 yr of age, the girl developed bilateral flexion deformity of interphalengeal joints and insulin-dependent diabetes mellitus. At age 7 yr, prominent hyperpigmented patches appeared on the skin at lower limbs, genitalia, and trunk. On clinical examination, she had hepatosplenomegaly, generalized lymphadenopathy, left ventricular hypertrophy, sensorineural hearing loss, hypogonadism, short stature, and characteristic dysmorphic features. Her brother had fixed flexion contractures of the feet, profound sensorineural hearing loss, characteristic dysmorphic features, but no diabetes. DNA sequence analysis revealed a homozygous mutation (c.300+1G>C) in SLC29A3 in both siblings. The phenotype and genotype of the siblings were compatible with that of the H syndrome, although the features were overlapping with those found in pigmented hypertrichotic dermatosis with insulin-dependent diabetes, familial SHML, and Faisalabad histiocytosis, indicating that these four syndromes may be regarded as one disease with varying phenotypic features. A new, common name for these conditions is warranted.

Exocrine pancreatic function in hepatocyte nuclear factor 1β-maturity-onset diabetes of the young (HNF1B-MODY) is only moderately reduced: compensatory hypersecretion from a hypoplastic pancreas

OBJECTIVES

To examine the exocrine pancreatic function in carriers of the hepatocyte nuclear factor 1β gene (HNF1B) mutation by direct testing.

METHODS

Patients with HNF1B mutations and control subjects were assessed using rapid endoscopic secretin tests and secretin-stimulated magnetic resonance imaging. Seven patients and 25 controls underwent endoscopy, while eight patients and 20 controls had magnetic resonance imaging. Ductal function was assessed according to peak bicarbonate concentrations and acinar function was assessed according to peak digestive enzyme activities in secretin-stimulated duodenal juice. The association of pancreatic exocrine function and diabetes status with pancreatic gland volume was examined.

RESULTS

The mean increase in secretin-stimulated duodenal fluid was smaller in patients than controls (4.0 vs 6.4 ml/min; P = 0.003). We found lower ductal function in patients than controls (median peak bicarbonate concentration: 73 vs 116 mEq/L; P < 0.001) and lower acinar function (median peak lipase activity: 6.4 vs 33.5 kU/ml; P = 0.01; median peak elastase activity: 0.056 vs 0.130 U/ml; P = 0.01). Pancreatic fluid volume outputs correlated significantly with pancreatic gland volumes (r² = 0.71, P = 0.008) in patients. The total fluid output to pancreatic gland volume ratios were higher in patients than controls (4.5 vs 1.3 ml/cm³; P = 0.03), suggesting compensatory hypersecretion in the remaining gland.

CONCLUSION

Carriers of the HNF1B mutation have lower exocrine pancreatic function involving both ductal and acinar cells. Compensatory hypersecretion suggests that the small pancreas of HNF1B mutation carriers is attributable to hypoplasia, not atrophy.

SHORT syndrome with partial lipodystrophy due to impaired phosphatidylinositol 3 kinase signaling

The phosphatidylinositol 3 kinase (PI3K) pathway regulates fundamental cellular processes such as metabolism, proliferation, and survival. A central component in this pathway is the p85α regulatory subunit, encoded by PIK3R1. Using whole-exome sequencing, we identified a heterozygous PIK3R1 mutation (c.1945C>T [p.Arg649Trp]) in two unrelated families affected by partial lipodystrophy, low body mass index, short stature, progeroid face, and Rieger anomaly (SHORT syndrome). This mutation led to impaired interaction between p85α and IRS-1 and reduced AKT-mediated insulin signaling in fibroblasts from affected subjects and in reconstituted Pik3r1-knockout preadipocytes. Normal PI3K activity is critical for adipose differentiation and insulin signaling; the mutated PIK3R1 therefore provides a unique link among lipodystrophy, growth, and insulin signaling.

Prevalence of monogenic diabetes in the population-based Norwegian Childhood Diabetes Registry

AIMS/HYPOTHESIS

Monogenic diabetes (MD) might be misdiagnosed as type 1 diabetes. The prevalence of MD among children with apparent type 1 diabetes has not been established. Our aim was to estimate the prevalence of common forms of MD in childhood diabetes.

METHODS

We investigated 2,756 children aged 0-14 years with newly diagnosed diabetes who had been recruited to the nationwide population-based Norwegian Childhood Diabetes Registry (NCDR), from July 2002 to March 2012. Completeness of ascertainment was 91%. Children diagnosed with diabetes who were under12 months of age were screened for mutations in KCNJ11, ABCC8 and INS. Children without GAD and protein tyrosine phosphatase-like protein antibodies were screened in two ways. Those who had a parent with diabetes were screened for mutations in HNF1A, HNF4A, INS and MT-TL1. Children with HbA1c <7.5% (<58 mmol/mol) and no insulin requirement were screened for mutations in GCK. Finally, we searched the Norwegian MODY Registry for children with genetically verified MD.

RESULTS

We identified 15 children harbouring a mutation in HNF1A, nine with one in GCK, four with one in KCNJ11, one child with a mutation in INS and none with a mutation in MT-TL1. The minimum prevalence of MD in the NCDR was therefore 1.1%. By searching the Norwegian MODY Registry, we found 24 children with glucokinase-MODY, 15 of whom were not present in the NCDR. We estimated the minimum prevalence of MD among Norwegian children to be 3.1/100,000.

CONCLUSIONS/INTERPRETATION

This is the first prevalence study of the common forms of MD in a nationwide, population-based registry of childhood diabetes. We found that 1.1% of patients in the Norwegian Childhood Diabetes Registry had MD.

GCK-MODY diabetes associated with protein misfolding, cellular self-association and degradation

GCK-MODY, dominantly inherited mild fasting hyperglycemia, has been associated with >600 different mutations in the glucokinase (GK)-encoding gene (GCK). When expressed as recombinant pancreatic proteins, some mutations result in enzymes with normal/near-normal catalytic properties. The molecular mechanism(s) of GCK-MODY due to these mutations has remained elusive. Here, we aimed to explore the molecular mechanisms for two such catalytically 'normal' GCK mutations (S263P and G264S) in the F260-L270 loop of GK. When stably overexpressed in HEK293 cells and MIN6 β-cells, the S263P- and G264S-encoded mutations generated misfolded proteins with an increased rate of degradation (S263P>G264S) by the protein quality control machinery, and a propensity to self-associate (G264S>S263P) and form dimers (SDS resistant) and aggregates (partly Triton X-100 insoluble), as determined by pulse-chase experiments and subcellular fractionation. Thus, the GCK-MODY mutations S263P and G264S lead to protein misfolding causing destabilization, cellular dimerization/aggregation and enhanced rate of degradation. In silico predicted conformational changes of the F260-L270 loop structure are considered to mediate the dimerization of both mutant proteins by a domain swapping mechanism. Thus, similar properties may represent the molecular mechanisms for additional unexplained GCK-MODY mutations, and may also contribute to the disease mechanism in other previously characterized GCK-MODY inactivating mutations.

To test, or not to test: time for a MODY calculator?

To test, or not to test, that is often the question in diabetes genetics. This is why the paper of Shields et al in the current issue of Diabetologia is so warmly welcomed. MODY is the most common form of monogenic diabetes. Nevertheless, the optimal way of identifying MODY families still poses a challenge both for researchers and clinicians. Hattersley's group in Exeter, UK, have developed an easy-to-use MODY prediction model that can help to identify cases appropriate for genetic testing. By answering eight simple questions on the internet ( www.diabetesgenes.org/content/mody-probability-calculator ), the doctor receives a positive predictive value in return: the probability that the patient has MODY. Thus, the classical binary (yes/no) assessment provided by clinical diagnostic criteria has been substituted by a more rational, quantitative estimate. The model appears to discriminate well between MODY and type 1 and type 2 diabetes when diabetes is diagnosed before the age of 35 years. However, the performance of the MODY probability calculator should now be validated in other settings than where it was developed-and, as always, there is room for some improvements and modifications.

Binding of ATP at the active site of human pancreatic glucokinase--nucleotide-induced conformational changes with possible implications for its kinetic cooperativity

Glucokinase (GK) is the central player in glucose-stimulated insulin release from pancreatic β-cells, and catalytic activation by α-d-glucose binding has a key regulatory function. Whereas the mechanism of this activation is well understood, on the basis of crystal structures of human GK, there are no similar structural data on ATP binding to the ligand-free enzyme and how it affects its conformation. Here, we report on a conformational change induced by the binding of adenine nucleotides to human pancreatic GK, as determined by intrinsic tryptophan fluorescence, using the catalytically inactive mutant form T228M to correct for the inner filter effect. Adenosine-5′-(β,γ-imido)triphosphate and ATP bind to the wild-type enzyme with apparent [L]_0.5 (ligand concentration at half-maximal effect) values of 0.27 ± 0.02 mm and 0.78 ± 0.14 mm, respectively. The change in protein conformation was further supported by ATP inhibition of the binding of the fluorescent probe 8-anilino-1-naphthalenesulfonate and limited proteolysis by trypsin, and by molecular dynamic simulations. The simulations provide a first insight into the dynamics of the binary complex with ATP, including motion of the flexible surface/active site loop and partial closure of the active site cleft. In the complex, the adenosine moiety is packed between two α-helices and stabilized by hydrogen bonds (with Thr228, Thr332, and Ser336) and hydrophobic interactions (with Val412 and Leu415). Combined with enzyme kinetic analyses, our data indicate that the ATP-induced changes in protein conformation may have implications for the kinetic cooperativity of the enzyme.

Pancreatic function in carboxyl-ester lipase knockout mice

BACKGROUND/AIMS

CEL-MODY is a monogenic form of diabetes and exocrine pancreatic insufficiency due to mutations in the carboxyl-ester lipase (CEL) gene. We aimed to investigate endocrine and exocrine pancreatic function in CEL knockout mice (CELKO).

METHODS

A knockout mouse model with global targeted deletion of CEL was investigated physiologically and histopathologically, and compared to littermate control CEL+/+ mice at 7 and 12 months on normal chow and high-fat diets (HFD), i.e. 42 and 60% fat by calories.

RESULTS

CELKO+/+ and -/- mice showed normal growth and development and normal glucose metabolism on a chow diet. Female CEL-/- mice on 60% HFD, on the other hand, had increased random blood glucose compared to littermate controls (p = 0.02), and this was accompanied by a reduction in glucose tolerance that did not reach statistical significance. In these mice there was also islet hyperplasia, however, α- and β-islet cells appeared morphologically normal and pancreatic exocrine function was also normal.

CONCLUSION

Although we observed mild glucose intolerance in female mice with whole-body knockout of CEL, the full phenotype of human CEL-MODY was not reproduced, suggesting that the pathogenic mechanisms involved are more complex than a simple loss of CEL function. and IAP.

Target setting in intensive insulin management is associated with metabolic control: the Hvidoere childhood diabetes study group centre differences study 2005

OBJECTIVE

To evaluate glycaemic targets set by diabetes teams, their perception by adolescents and parents, and their influence on metabolic control.

METHODS

Clinical data and questionnaires were completed by adolescents, parents/carers and diabetes teams in 21 international centres. HbA1c was measured centrally.

RESULTS

A total of 2062 adolescents completed questionnaires (age 14.4 +/- 2.3 yr; diabetes duration 6.1 +/- 3.5 yr). Mean HbA 1c = 8.2 +/- 1.4% with significant differences between centres (F = 12.3; p < 0.001) range from 7.4 to 9.1%. There was a significant correlation between parent (r = 0.20) and adolescent (r = 0.21) reports of their perceived ideal HbA1c and their actual HbA1c result (p < 0.001), and a stronger association between parents' (r = 0.39) and adolescents' (r = 0.4) reports of the HbA1c they would be happy with and their actual HbA1c result. There were significant differences between centres on parent and adolescent reports of ideal and happy with HbA1c (8.1 < F > 17.4;p < 0.001). A lower target HbA1c and greater consistency between members of teams within centres were associated with lower centre HbA1c (F = 16.0; df = 15; p < 0.001).

CONCLUSIONS

Clear and consistent setting of glycaemic targets by diabetes teams is strongly associated with HbA1c outcome in adolescents. Target setting appears to play a significant role in explaining the differences in metabolic outcomes between centres.

Clinical spectrum of SIX3-associated mutations in holoprosencephaly: correlation between genotype, phenotype and function

BACKGROUND

Holoprosencephaly (HPE) is the most common structural malformation of the human forebrain. There are several important HPE mutational target genes, including the transcription factor SIX3, which encodes an early regulator of Shh, Wnt, Bmp and Nodal signalling expressed in the developing forebrain and eyes of all vertebrates.

OBJECTIVE

To characterise genetic and clinical findings in patients with SIX3 mutations.

METHODS

Patients with HPE and their family members were tested for mutations in HPE-associated genes and the genetic and clinical findings, including those for additional cases found in the literature, were analysed. The results were correlated with a mutation-specific functional assay in zebrafish.

RESULTS

In a cohort of patients (n = 800) with HPE, SIX3 mutations were found in 4.7% of probands and additional cases were found through testing of relatives. In total, 138 cases of HPE were identified, 59 of whom had not previously been clinically presented. Mutations in SIX3 result in more severe HPE than in other cases of non-chromosomal, non-syndromic HPE. An over-representation of severe HPE was found in patients whose mutations confer greater loss of function, as measured by the functional zebrafish assay. The gender ratio in this combined set of patients was 1.5:1 (F:M) and maternal inheritance was almost twice as common as paternal. About 14% of SIX3 mutations in probands occur de novo. There is a wide intrafamilial clinical range of features and classical penetrance is estimated to be at least 62%.

CONCLUSIONS

Our data suggest that SIX3 mutations result in relatively severe HPE and that there is a genotype-phenotype correlation, as shown by functional studies using animal models.

A TLR2 polymorphism is associated with type 1 diabetes and allergic asthma

Type 1 diabetes (T1D) and allergic asthma are immune-mediated diseases. Pattern recognition receptors are proteins expressed by cells in the immune system to identify microbial pathogens and endogenous ligands. Toll-like receptors (TLRs) and CD14 are members of this family and could represent a molecular link between microbial infections and immune-mediated diseases. Diverging hypotheses regarding whether there exists a common or inverse genetic etiology behind these immune-mediated diseases have been presented. We aimed to test whether there exist common or inverse associations between polymorphisms in the pattern recognition receptors TLR2, TLR4 and CD14 and T1D and allergic asthma. Eighteen single nucleotide polymorphisms (SNPs) were genotyped in TLR2 (2), TLR4 (12) and CD14 (4) in 700 T1D children, 357 nuclear families with T1D children and 796 children from the 'Environment and Childhood Asthma' study. Allele and haplotype frequencies were analyzed in relation to diseases and in addition transmission disequilibrium test analyses were performed in the family material. Both T1D and allergic asthma were significantly associated with the TLR2 rs3804100 T allele and further associated with the haplotype including this SNP, possibly representing a susceptibility locus common for the two diseases. Neither TLR4 nor CD14 were associated with T1D or allergic asthma.

Wolcott-Rallison syndrome with 3-hydroxydicarboxylic aciduria and lethal outcome

Wolcott-Rallison syndrome (WRS) (OMIM 226980) is a rare, autosomal recessive disorder with infancy-onset diabetes mellitus, multiple epiphyseal dysplasia, osteopenia, mental retardation or developmental delay, and hepatic and renal dysfunction as main clinical findings. Patients with WRS have mutations in the EIF2AK3 gene, which encodes the pancreatic eukaryotic translation initiation factor 2-alpha kinase 3. We report a female patient who developed insulin-requiring diabetes at 2.5 months of age. Multiple epiphyseal dysplasia was diagnosed at age 2 years. At age 5.5 years she developed a Reye-like syndrome with hypoketotic hypoglycaemia and renal and hepatic insufficiency and died. A partial autopsy showed fat infiltration in the liver and kidneys. Examination of urine by gas chromatography and mass spectrometry showed large amounts of C(6)-dicarboxylic acid (adipic acid), 3-hydroxy-C(8)-dicarboxylic acid, 3-hydroxy-C(10)-dicarboxylic acid, and 3-hydroxydecenedioic acid. Acetoacetate and 3-hydroxybutyrate were absent. The findings suggested a metabolic block in mitochondrial fatty acid oxidation, but lack of material precluded enzyme analyses. The clinical diagnosis of WRS was suggested in retrospect, and confirmed by sequencing of DNA extracted from stored autopsy material. The patient was compound heterozygous for the novel EIF2AK3 mutations c.1694_1695delAT (Y565X) and c.3044T > C (F1015S). Our data suggest that disruption of the EIF2AK3 gene may lead to defective mitochondrial fatty acid oxidation and hypoglycaemia, thus adding to the heterogeneous phenotype of WRS.

A coding polymorphism in NALP1 confers risk for autoimmune Addison's disease and type 1 diabetes

Variants in the gene encoding NACHT leucine-rich-repeat protein 1 (NALP1), an important molecule in innate immunity, have recently been shown to confer risk for vitiligo and associated autoimmunity. We hypothesized that sequence variants in this gene may be involved in susceptibility to a wider spectrum of autoimmune diseases. Investigating large patient cohorts from six different autoimmune diseases, that is autoimmune Addison's disease (n=333), type 1 diabetes (n=1086), multiple sclerosis (n=502), rheumatoid arthritis (n=945), systemic lupus erythematosus (n=156) and juvenile idiopathic arthritis (n=505), against 3273 healthy controls, we analyzed four single nucleotide polymorphisms (SNPs) in NALP1. The major allele of the coding SNP rs12150220 revealed significant association with autoimmune Addison's disease compared with controls (OR=1.25, 95% CI: 1.06-1.49, P=0.007), and with type 1 diabetes (OR=1.15, 95% CI: 1.04-1.27, P=0.005). Trends toward the same associations were seen in rheumatoid arthritis, systemic lupus erythematosus and, although less obvious, multiple sclerosis. Patients with juvenile idiopathic arthritis did not show association with NALP1 gene variants. The results indicate that NALP1 and the innate immune system may be implicated in the pathogenesis of many autoimmune disorders, particularly organ-specific autoimmune diseases.

Prevalence of HNF1A (MODY3) mutations in a Norwegian population (the HUNT2 Study)

AIMS

Previous reports have indicated that maturity-onset diabetes of the young (MODY) caused by hepatocyte nuclear factor 1A (HNF1A) mutations (MODY3) is the most common MODY subtype in Northern Europe, but population-based prevalence estimates are lacking. We sought to determine the prevalence of HNF1A-MODY in diabetic subjects of a defined Norwegian population (the HUNT2 Study).

METHODS

Of the 1972 diabetic HUNT2 subjects, we identified a subgroup of 43 suspected MODY cases based on information on family history, disease onset and anti-glutamic acid decarboxylase autoantibody status. These cases were considered a discovery group for HNF1A mutations and underwent full DNA sequencing. Subsequently, the entire cohort of diabetic HUNT2 subjects was screened for three selected HNF1A mutations. Possible founder effects were examined using the Norwegian MODY Registry.

RESULTS

Three subjects from the discovery group harboured HNF1A mutations. Two subjects had the previously described R229Q mutation, one had a novel S6N alteration, whereas the HNF1A hot-spot mutation P291fsinsC was not identified. Genotyping the cohort of diabetic HUNT2 subjects identified five additional R229Q-positive subjects. Microsatellite analysis performed for all R229Q-positive probands of the Norwegian MODY Registry and those found in the HUNT2 population revealed that 17 of 18 (94%) had genotypes consistent with a common haplotype.

CONCLUSIONS

Clinical MODY criteria were fulfilled in 2.2% of diabetic HUNT2 subjects. The minimum prevalence of HNF1A-MODY among diabetic HUNT2 subjects was 0.4%. Because of founder effects, registry-based prevalence studies probably need to be very large and they should also include prospectively collected phenotypes and extensive mutation screening to establish the true prevalence of MODY.

Lack of pancreatic body and tail in HNF1B mutation carriers

AIMS

Hepatocyte nuclear factor 1B (HNF1B) gene mutation carriers have a systemic disease characterized by congenital malformations in the urogenital tract, diabetes mellitus of maturity-onset diabetes of the young type and dysfunction of the liver and exocrine pancreas. We aimed to investigate pancreatic structure and exocrine function in carriers of HNF1B mutations.

METHODS

We studied five subjects from two families with the previously reported mutation R137_K161del and the novel mutation F148L in HNF1B. All patients underwent computed tomography (CT) and magnetic resonance cholangiopancreatography (MRCP). We measured faecal elastase and serum vitamins D and E.

RESULTS

One of the mutation carriers reported abdominal symptoms. All five subjects had faecal elastase deficiency, three had vitamin D deficiency and two had vitamin E deficiency. Neither CT nor MRCP depicted tissue corresponding to the pancreatic body and tail in the five mutation carriers, indicating agenesis of the dorsal pancreas. The head of the pancreas was slightly atrophic but had normal X-ray attenuation at CT in all patients.

CONCLUSIONS

Agenesis of the pancreatic body and tail and pancreatic exocrine dysfunction are parts of the phenotype in HNF1B mutation carriers. This strengthens the evidence for a critical role of HNF1B in development and differentiation of at least the dorsal pancreas.

Genetic analysis of recently identified type 2 diabetes loci in 1,638 unselected patients with type 2 diabetes and 1,858 control participants from a Norwegian population-based cohort (the HUNT study)

AIMS/HYPOTHESIS

Recent genome-wide association studies performed in selected patients and control participants have provided strong support for several new type 2 diabetes susceptibility loci. To get a better estimation of the true risk conferred by these novel loci, we tested a completely unselected population of type 2 diabetes patients from a Norwegian health survey (the HUNT study).

METHODS

We genotyped single nucleotide polymorphisms (SNPs) in PKN2, IGFBP2, FLJ39370 (also known as C4ORF32), CDKAL1, SLC30A8, CDKN2B, HHEX and FTO using a Norwegian population-based sample of 1,638 patients with type 2 diabetes and 1,858 non-diabetic control participants (the HUNT Study), for all of whom data on BMI, WHR, cholesterol and triacylglycerol levels were available. We used diabetes, measures of obesity and lipid values as phenotypes in case-control and quantitative association study designs.

RESULTS

We replicated the association with type 2 diabetes for rs10811661 in the vicinity of CDKN2B (OR 1.20, 95% CI: 1.06-1.37, p=0.004), rs9939609 in FTO (OR 1.14, 95% CI: 1.04-1.25, p=0.006) and rs13266634 in SLC30A8 (OR 1.20, 95% CI: 1.09-1.33, p=3.9 x 10(-4)). We found borderline significant association for the IGFBP2 SNP rs4402960 (OR 1.10, 95% CI: 0.99-1.22). Results for the HHEX SNP (rs1111875) and the CDKAL1 SNP (rs7756992) were non-significant, but the magnitude of effect was similar to previous estimates. We found no support for an association with the less consistently replicated FLJ39370 or PKN2 SNPs. In agreement with previous studies, FTO was most strongly associated with BMI (p=8.4 x 10(-4)).

CONCLUSIONS/INTERPRETATION

Our data show that SNPs near IGFBP2, CDKAL1, SLC30A8, CDKN2B, HHEX and FTO are also associated with diabetes in non-selected patients with type 2 diabetes.