Patient-derived organoid biobank identifies epigenetic dysregulation of intestinal epithelial MHC-I as a novel mechanism in severe Crohn's Disease.

OBJECTIVE: Epigenetic mechanisms, including DNA methylation (DNAm), have been proposed to play a key role in Crohn's disease (CD) pathogenesis. However, the specific cell types and pathways affected as well as their potential impact on disease phenotype and outcome remain unknown. We set out to investigate the role of intestinal epithelial DNAm in CD pathogenesis. DESIGN: We generated 312 intestinal epithelial organoids (IEOs) from mucosal biopsies of 168 patients with CD (n=72), UC (n=23) and healthy controls (n=73). We performed genome-wide molecular profiling including DNAm, bulk as well as single-cell RNA sequencing. Organoids were subjected to gene editing and the functional consequences of DNAm changes evaluated using an organoid-lymphocyte coculture and a nucleotide-binding oligomerisation domain, leucine-rich repeat and CARD domain containing 5 (NLRC5) dextran sulphate sodium (DSS) colitis knock-out mouse model. RESULTS: We identified highly stable, CD-associated loss of DNAm at major histocompatibility complex (MHC) class 1 loci including NLRC5 and cognate gene upregulation. Single-cell RNA sequencing of primary mucosal tissue and IEOs confirmed the role of NLRC5 as transcriptional transactivator in the intestinal epithelium. Increased mucosal MHC-I and NLRC5 expression in adult and paediatric patients with CD was validated in additional cohorts and the functional role of MHC-I highlighted by demonstrating a relative protection from DSS-mediated mucosal inflammation in NLRC5-deficient mice. MHC-I DNAm in IEOs showed a significant correlation with CD disease phenotype and outcomes. Application of machine learning approaches enabled the development of a disease prognostic epigenetic molecular signature. CONCLUSIONS: Our study has identified epigenetically regulated intestinal epithelial MHC-I as a novel mechanism in CD pathogenesis.

An expanded clinical spectrum of hypoinsulinaemic hypoketotic hypoglycaemia.

BACKGROUND: Hypoketotic hypoglycaemia with suppressed plasma fatty acids and detectable insulin suggests congenital hyperinsulinism (CHI). Severe hypoketotic hypoglycaemia mimicking hyperinsulinism but without detectable insulin has recently been described in syndromic individuals with mosaic genetic activation of post-receptor insulin signalling. We set out to expand understanding of this entity focusing on metabolic phenotypes. METHODS: Metabolic profiling, candidate gene and exome sequencing were performed in six infants with hypoketotic, hypoinsulinaemic hypoglycaemia, with or without syndromic features. Additional signalling studies were carried out in dermal fibroblasts from two individuals. RESULTS: Two infants had no syndromic features. One was mistakenly diagnosed with CHI. One had mild features of megalencephaly-capillary malformation-polymicrogyria (MCAP) syndrome, one had non-specific macrosomia, and two had complex syndromes. All required intensive treatment to maintain euglycaemia, with CHI-directed therapies being ineffective. Pathogenic PIK3CA variants were found in two individuals - de novo germline c.323G>A (p.Arg108His) in one non-syndromic infant and postzygotic mosaic c.2740G>A (p.Gly914Arg) in the infant with MCAP. No causal variants were proven in the other individuals despite extensive investigation, although rare variants in mTORC components were identified in one. No increased PI3K signalling in fibroblasts of two individuals was seen. CONCLUSIONS: We expand the spectrum of PI3K-related hypoinsulinaemic hypoketotic hypoglycaemia. We demonstrate that pathogenic germline variants activating post-insulin-receptor signalling may cause non-syndromic hypoinsulinaemic hypoketotic hypoglycaemia closely resembling CHI. This distinct biochemical footprint should be sought and differentiated from CHI in infantile hypoglycaemia. To facilitate adoption of this differential diagnosis, we propose the term "pseudohyperinsulinism".

Culture-Associated DNA Methylation Changes Impact on Cellular Function of Human Intestinal Organoids.

BACKGROUND & AIMS: Human intestinal epithelial organoids (IEOs) are a powerful tool to model major aspects of intestinal development, health, and diseases because patient-derived cultures retain many features found in vivo. A necessary aspect of the organoid model is the requirement to expand cultures in vitro through several rounds of passaging. This is of concern because the passaging of cells has been shown to affect cell morphology, ploidy, and function. METHODS: Here, we analyzed 173 human IEO lines derived from the small and large bowel and examined the effect of culture duration on DNA methylation (DNAm). Furthermore, we tested the potential impact of DNAm changes on gene expression and cellular function. RESULTS: Our analyses show a reproducible effect of culture duration on DNAm in a large discovery cohort as well as 2 publicly available validation cohorts generated in different laboratories. Although methylation changes were seen in only approximately 8% of tested cytosine-phosphate-guanine dinucleotides (CpGs) and global cellular function remained stable, a subset of methylation changes correlated with altered gene expression at baseline as well as in response to inflammatory cytokine exposure and withdrawal of Wnt agonists. Importantly, epigenetic changes were found to be enriched in genomic regions associated with colonic cancer and distant to the site of replication, indicating similarities to malignant transformation. CONCLUSIONS: Our study shows distinct culture-associated epigenetic changes in mucosa-derived human IEOs, some of which appear to impact gene transcriptomic and cellular function. These findings highlight the need for future studies in this area and the importance of considering passage number as a potentially confounding factor.

Disease-associated DNA methylation signatures in esophageal biopsies of children diagnosed with Eosinophilic Esophagitis.

Eosinophilic esophagitis (EoE) is a leading cause of dysphagia and food impaction in children and adults. The diagnosis relies on histological examination of esophageal mucosal biopsies and requires the presence of > 15 eosinophils per high-powered field. Potential pitfalls include the impact of biopsy sectioning as well as regional variations of eosinophil density. We performed genome-wide DNA methylation analyses on 30 esophageal biopsies obtained from children diagnosed with EoE (n = 7) and matched controls (n = 13) at the time of diagnosis as well as following first-line treatment. Analyses revealed striking disease-associated differences in mucosal DNA methylation profiles in children diagnosed with EoE, highlighting the potential for these epigenetic signatures to be developed into clinically applicable biomarkers.

Transcription and DNA Methylation Patterns of Blood-Derived CD8+ T Cells Are Associated With Age and Inflammatory Bowel Disease But Do Not Predict Prognosis.

BACKGROUND & AIMS: Gene expression patterns of CD8+ T cells have been reported to correlate with clinical outcomes of adults with inflammatory bowel diseases (IBD). We aimed to validate these findings in independent patient cohorts. METHODS: We obtained peripheral blood samples from 112 children with a new diagnosis of IBD (71 with Crohn's disease and 41 with ulcerative colitis) and 19 children without IBD (controls) and recorded medical information on disease activity and outcomes. CD8+ T cells were isolated from blood samples by magnetic bead sorting at the point of diagnosis and during the course of disease. Genome-wide transcription (n = 192) and DNA methylation (n = 66) profiles were generated using Affymetrix and Illumina arrays, respectively. Publicly available transcriptomes and DNA methylomes of CD8+ T cells from 3 adult patient cohorts with and without IBD were included in data analyses. RESULTS: Previously reported CD8+ T-cell prognostic expression and exhaustion signatures were only found in the original adult IBD patient cohort. These signatures could not be detected in either a pediatric or a second adult IBD cohort. In contrast, an association between CD8+ T-cell gene expression with age and sex was detected across all 3 cohorts. CD8+ gene transcription was clearly associated with IBD in the 2 cohorts that included non-IBD controls. Lastly, DNA methylation profiles of CD8+ T cells from children with Crohn's disease correlated with age but not with disease outcome. CONCLUSIONS: We were unable to validate previously reported findings of an association between CD8+ T-cell gene transcription and disease outcome in IBD. Our findings reveal the challenges of developing prognostic biomarkers for patients with IBD and the importance of their validation in large, independent cohorts before clinical application.

Exome Sequencing Identifies Genes and Gene Sets Contributing to Severe Childhood Obesity, Linking PHIP Variants to Repressed POMC Transcription.

Obesity is genetically heterogeneous with monogenic and complex polygenic forms. Using exome and targeted sequencing in 2,737 severely obese cases and 6,704 controls, we identified three genes (PHIP, DGKI, and ZMYM4) with an excess burden of very rare predicted deleterious variants in cases. In cells, we found that nuclear PHIP (pleckstrin homology domain interacting protein) directly enhances transcription of pro-opiomelanocortin (POMC), a neuropeptide that suppresses appetite. Obesity-associated PHIP variants repressed POMC transcription. Our demonstration that PHIP is involved in human energy homeostasis through transcriptional regulation of central melanocortin signaling has potential diagnostic and therapeutic implications for patients with obesity and developmental delay. Additionally, we found an excess burden of predicted deleterious variants involving genes nearest to loci from obesity genome-wide association studies. Genes and gene sets influencing obesity with variable penetrance provide compelling evidence for a continuum of causality in the genetic architecture of obesity, and explain some of its missing heritability.

The metabolic syndrome- associated small G protein ARL15 plays a role in adipocyte differentiation and adiponectin secretion.

Common genetic variants at the ARL15 locus are associated with plasma adiponectin, insulin and HDL cholesterol concentrations, obesity, and coronary atherosclerosis. The ARL15 gene encodes a small GTP-binding protein whose function is currently unknown. In this study adipocyte-autonomous roles for ARL15 were investigated using conditional knockdown of Arl15 in murine 3T3-L1 (pre)adipocytes. Arl15 knockdown in differentiated adipocytes impaired adiponectin secretion but not adipsin secretion or insulin action, while in preadipocytes it impaired adipogenesis. In differentiated adipocytes GFP-tagged ARL15 localized predominantly to the Golgi with lower levels detected at the plasma membrane and intracellular vesicles, suggesting involvement in intracellular trafficking. Sequencing of ARL15 in 375 severely insulin resistant patients identified four rare heterozygous variants, including an early nonsense mutation in a proband with femorogluteal lipodystrophy and non classical congenital adrenal hyperplasia, and an essential splice site mutation in a proband with partial lipodystrophy and a history of childhood yolk sac tumour. No nonsense or essential splice site mutations were found in 2,479 controls, while five such variants were found in the ExAC database. These findings provide evidence that ARL15 plays a role in adipocyte differentiation and adiponectin secretion, and raise the possibility that human ARL15 haploinsufficiency predisposes to lipodystrophy.

Hypoinsulinaemic, hypoketotic hypoglycaemia due to mosaic genetic activation of PI3-kinase.

OBJECTIVE: Genetic activation of the insulin signal-transducing kinase AKT2 causes syndromic hypoketotic hypoglycaemia without elevated insulin. Mosaic activating mutations in class 1A phospatidylinositol-3-kinase (PI3K), upstream from AKT2 in insulin signalling, are known to cause segmental overgrowth, but the metabolic consequences have not been systematically reported. We assess the metabolic phenotype of 22 patients with mosaic activating mutations affecting PI3K, thereby providing new insight into the metabolic function of this complex node in insulin signal transduction. METHODS: Three patients with megalencephaly, diffuse asymmetric overgrowth, hypoketotic, hypoinsulinaemic hypoglycaemia and no AKT2 mutation underwent further genetic, clinical and metabolic investigation. Signalling in dermal fibroblasts from one patient and efficacy of the mTOR inhibitor Sirolimus on pathway activation were examined. Finally, the metabolic profile of a cohort of 19 further patients with mosaic activating mutations in PI3K was assessed. RESULTS: In the first three patients, mosaic mutations in PIK3CA (p.Gly118Asp or p.Glu726Lys) or PIK3R2 (p.Gly373Arg) were found. In different tissue samples available from one patient, the PIK3CA p.Glu726Lys mutation was present at burdens from 24% to 42%, with the highest level in the liver. Dermal fibroblasts showed increased basal AKT phosphorylation which was potently suppressed by Sirolimus. Nineteen further patients with mosaic mutations in PIK3CA had neither clinical nor biochemical evidence of hypoglycaemia. CONCLUSIONS: Mosaic mutations activating class 1A PI3K cause severe non-ketotic hypoglycaemia in a subset of patients, with the metabolic phenotype presumably related to the extent of mosaicism within the liver. mTOR or PI3K inhibitors offer the prospect for future therapy.

Human biallelic MFN2 mutations induce mitochondrial dysfunction, upper body adipose hyperplasia, and suppression of leptin expression.

MFN2 encodes mitofusin 2, a membrane-bound mediator of mitochondrial membrane fusion and inter-organelle communication. MFN2 mutations cause axonal neuropathy, with associated lipodystrophy only occasionally noted, however homozygosity for the p.Arg707Trp mutation was recently associated with upper body adipose overgrowth. We describe similar massive adipose overgrowth with suppressed leptin expression in four further patients with biallelic MFN2 mutations and at least one p.Arg707Trp allele. Overgrown tissue was composed of normal-sized, UCP1-negative unilocular adipocytes, with mitochondrial network fragmentation, disorganised cristae, and increased autophagosomes. There was strong transcriptional evidence of mitochondrial stress signalling, increased protein synthesis, and suppression of signatures of cell death in affected tissue, whereas mitochondrial morphology and gene expression were normal in skin fibroblasts. These findings suggest that specific MFN2 mutations cause tissue-selective mitochondrial dysfunction with increased adipocyte proliferation and survival, confirm a novel form of excess adiposity with paradoxical suppression of leptin expression, and suggest potential targeted therapies.

A two-stage inter-rater approach for enrichment testing of variants associated with multiple traits.

Shared genetic aetiology may explain the co-occurrence of diseases in individuals more often than expected by chance. On identifying associated variants shared between two traits, one objective is to determine whether such overlap may be explained by specific genomic characteristics (eg, functional annotation). In clinical studies, inter-rater agreement approaches assess concordance among expert opinions on the presence/absence of a complex disease for each subject. We adapt a two-stage inter-rater agreement model to the genetic association setting to identify features predictive of overlap variants, while accounting for their marginal trait associations. The resulting corrected overlap and marginal enrichment test (COMET) also assesses enrichment at the individual trait level. Multiple categories may be tested simultaneously and the method is computationally efficient, not requiring permutations to assess significance. In an extensive simulation study, COMET identifies features predictive of enrichment with high power and has well-calibrated type I error. In contrast, testing for overlap with a single-trait enrichment test has inflated type I error. COMET is applied to three glycaemic traits using a set of functional annotation categories as predictors, followed by further analyses that focus on tissue-specific regulatory variants. The results support previous findings that regulatory variants in pancreatic islets are enriched for fasting glucose-associated variants, and give insight into differences/similarities between characteristics of variants associated with glycaemic traits. Also, despite regulatory variants in pancreatic islets being enriched for variants that are marginally associated with fasting glucose and fasting insulin, there is no enrichment of shared variants between the traits.

Integrative genomic analysis implicates limited peripheral adipose storage capacity in the pathogenesis of human insulin resistance.

Insulin resistance is a key mediator of obesity-related cardiometabolic disease, yet the mechanisms underlying this link remain obscure. Using an integrative genomic approach, we identify 53 genomic regions associated with insulin resistance phenotypes (higher fasting insulin levels adjusted for BMI, lower HDL cholesterol levels and higher triglyceride levels) and provide evidence that their link with higher cardiometabolic risk is underpinned by an association with lower adipose mass in peripheral compartments. Using these 53 loci, we show a polygenic contribution to familial partial lipodystrophy type 1, a severe form of insulin resistance, and highlight shared molecular mechanisms in common/mild and rare/severe insulin resistance. Population-level genetic analyses combined with experiments in cellular models implicate CCDC92, DNAH10 and L3MBTL3 as previously unrecognized molecules influencing adipocyte differentiation. Our findings support the notion that limited storage capacity of peripheral adipose tissue is an important etiological component in insulin-resistant cardiometabolic disease and highlight genes and mechanisms underpinning this link.

Insulin resistance uncoupled from dyslipidemia due to C-terminal PIK3R1 mutations.

Obesity-related insulin resistance is associated with fatty liver, dyslipidemia, and low plasma adiponectin. Insulin resistance due to insulin receptor (INSR) dysfunction is associated with none of these, but when due to dysfunction of the downstream kinase AKT2 phenocopies obesity-related insulin resistance. We report 5 patients with SHORT syndrome and C-terminal mutations in PIK3R1, encoding the p85α/p55α/p50α subunits of PI3K, which act between INSR and AKT in insulin signaling. Four of 5 patients had extreme insulin resistance without dyslipidemia or hepatic steatosis. In 3 of these 4, plasma adiponectin was preserved, as in insulin receptor dysfunction. The fourth patient and her healthy mother had low plasma adiponectin associated with a potentially novel mutation, p.Asp231Ala, in adiponectin itself. Cells studied from one patient with the p.Tyr657X PIK3R1 mutation expressed abundant truncated PIK3R1 products and showed severely reduced insulin-stimulated association of mutant but not WT p85α with IRS1, but normal downstream signaling. In 3T3-L1 preadipocytes, mutant p85α overexpression attenuated insulin-induced AKT phosphorylation and adipocyte differentiation. Thus, PIK3R1 C-terminal mutations impair insulin signaling only in some cellular contexts and produce a subphenotype of insulin resistance resembling INSR dysfunction but unlike AKT2 dysfunction, implicating PI3K in the pathogenesis of key components of the metabolic syndrome.

Conjunctival fibrosis and the innate barriers to Chlamydia trachomatis intracellular infection: a genome wide association study.

Chlamydia trachomatis causes both trachoma and sexually transmitted infections. These diseases have similar pathology and potentially similar genetic predisposing factors. We aimed to identify polymorphisms and pathways associated with pathological sequelae of ocular Chlamydia trachomatis infections in The Gambia. We report a discovery phase genome-wide association study (GWAS) of scarring trachoma (1090 cases, 1531 controls) that identified 27 SNPs with strong, but not genome-wide significant, association with disease (5 × 10(-6) > P > 5 × 10(-8)). The most strongly associated SNP (rs111513399, P = 5.38 × 10(-7)) fell within a gene (PREX2) with homology to factors known to facilitate chlamydial entry to the host cell. Pathway analysis of GWAS data was significantly enriched for mitotic cell cycle processes (P = 0.001), the immune response (P = 0.00001) and for multiple cell surface receptor signalling pathways. New analyses of published transcriptome data sets from Gambia, Tanzania and Ethiopia also revealed that the same cell cycle and immune response pathways were enriched at the transcriptional level in various disease states. Although unconfirmed, the data suggest that genetic associations with chlamydial scarring disease may be focussed on processes relating to the immune response, the host cell cycle and cell surface receptor signalling.

Truncation of POC1A associated with short stature and extreme insulin resistance.

We describe a female proband with primordial dwarfism, skeletal dysplasia, facial dysmorphism, extreme dyslipidaemic insulin resistance and fatty liver associated with a novel homozygous frameshift mutation in POC1A, predicted to affect two of the three protein products of the gene. POC1A encodes a protein associated with centrioles throughout the cell cycle and implicated in both mitotic spindle and primary ciliary function. Three homozygous mutations affecting all isoforms of POC1A have recently been implicated in a similar syndrome of primordial dwarfism, although no detailed metabolic phenotypes were described. Primary cells from the proband we describe exhibited increased centrosome amplification and multipolar spindle formation during mitosis, but showed normal DNA content, arguing against mitotic skipping, cleavage failure or cell fusion. Despite evidence of increased DNA damage in cells with supernumerary centrosomes, no aneuploidy was detected. Extensive centrosome clustering both at mitotic spindles and in primary cilia mitigated the consequences of centrosome amplification, and primary ciliary formation was normal. Although further metabolic studies of patients with POC1A mutations are warranted, we suggest that POC1A may be added to ALMS1 and PCNT as examples of centrosomal or pericentriolar proteins whose dysfunction leads to extreme dyslipidaemic insulin resistance. Further investigation of links between these molecular defects and adipose tissue dysfunction is likely to yield insights into mechanisms of adipose tissue maintenance and regeneration that are critical to metabolic health.

Hypomorphism in human NSMCE2 linked to primordial dwarfism and insulin resistance.

Structural maintenance of chromosomes (SMC) complexes are essential for maintaining chromatin structure and regulating gene expression. Two the three known SMC complexes, cohesin and condensin, are important for sister chromatid cohesion and condensation, respectively; however, the function of the third complex, SMC5-6, which includes the E3 SUMO-ligase NSMCE2 (also widely known as MMS21) is less clear. Here, we characterized 2 patients with primordial dwarfism, extreme insulin resistance, and gonadal failure and identified compound heterozygous frameshift mutations in NSMCE2. Both mutations reduced NSMCE2 expression in patient cells. Primary cells from one patient showed increased micronucleus and nucleoplasmic bridge formation, delayed recovery of DNA synthesis, and reduced formation of foci containing Bloom syndrome helicase (BLM) after hydroxyurea-induced replication fork stalling. These nuclear abnormalities in patient dermal fibroblast were restored by expression of WT NSMCE2, but not a mutant form lacking SUMO-ligase activity. Furthermore, in zebrafish, knockdown of the NSMCE2 ortholog produced dwarfism, which was ameliorated by reexpression of WT, but not SUMO-ligase-deficient NSMCE. Collectively, these findings support a role for NSMCE2 in recovery from DNA damage and raise the possibility that loss of its function produces dwarfism through reduced tolerance of replicative stress.

Mutations disrupting the Kennedy phosphatidylcholine pathway in humans with congenital lipodystrophy and fatty liver disease.

Phosphatidylcholine (PC) is the major glycerophospholipid in eukaryotic cells and is an essential component in all cellular membranes. The biochemistry of de novo PC synthesis by the Kennedy pathway is well established, but less is known about the physiological functions of PC. We identified two unrelated patients with defects in the Kennedy pathway due to biallellic loss-of-function mutations in phosphate cytidylyltransferase 1 alpha (PCYT1A), the rate-limiting enzyme in this pathway. The mutations lead to a marked reduction in PCYT1A expression and PC synthesis. The phenotypic consequences include some features, such as severe fatty liver and low HDL cholesterol levels, that are predicted by the results of previously reported liver-specific deletion of murine Pcyt1a. Both patients also had lipodystrophy, severe insulin resistance, and diabetes, providing evidence for an additional and essential role for PCYT1A-generated PC in the normal function of white adipose tissue and insulin action.

Impact of type 2 diabetes susceptibility variants on quantitative glycemic traits reveals mechanistic heterogeneity.

Patients with established type 2 diabetes display both ß-cell dysfunction and insulin resistance. To define fundamental processes leading to the diabetic state, we examined the relationship between type 2 diabetes risk variants at 37 established susceptibility loci, and indices of proinsulin processing, insulin secretion, and insulin sensitivity. We included data from up to 58,614 nondiabetic subjects with basal measures and 17,327 with dynamic measures. We used additive genetic models with adjustment for sex, age, and BMI, followed by fixed-effects, inverse-variance meta-analyses. Cluster analyses grouped risk loci into five major categories based on their relationship to these continuous glycemic phenotypes. The first cluster (PPARG, KLF14, IRS1, GCKR) was characterized by primary effects on insulin sensitivity. The second cluster (MTNR1B, GCK) featured risk alleles associated with reduced insulin secretion and fasting hyperglycemia. ARAP1 constituted a third cluster characterized by defects in insulin processing. A fourth cluster (TCF7L2, SLC30A8, HHEX/IDE, CDKAL1, CDKN2A/2B) was defined by loci influencing insulin processing and secretion without a detectable change in fasting glucose levels. The final group contained 20 risk loci with no clear-cut associations to continuous glycemic traits. By assembling extensive data on continuous glycemic traits, we have exposed the diverse mechanisms whereby type 2 diabetes risk variants impact disease predisposition.

Mosaic overgrowth with fibroadipose hyperplasia is caused by somatic activating mutations in PIK3CA.

The phosphatidylinositol 3-kinase (PI3K)-AKT signaling pathway is critical for cellular growth and metabolism. Correspondingly, loss of function of PTEN, a negative regulator of PI3K, or activating mutations in AKT1, AKT2 or AKT3 have been found in distinct disorders featuring overgrowth or hypoglycemia. We performed exome sequencing of DNA from unaffected and affected cells from an individual with an unclassified syndrome of congenital progressive segmental overgrowth of fibrous and adipose tissue and bone and identified the cancer-associated mutation encoding p.His1047Leu in PIK3CA, the gene that encodes the p110α catalytic subunit of PI3K, only in affected cells. Sequencing of PIK3CA in ten additional individuals with overlapping syndromes identified either the p.His1047Leu alteration or a second cancer-associated alteration, p.His1047Arg, in nine cases. Affected dermal fibroblasts showed enhanced basal and epidermal growth factor (EGF)-stimulated phosphatidylinositol 3,4,5-trisphosphate (PIP(3)) generation and concomitant activation of downstream signaling relative to their unaffected counterparts. Our findings characterize a distinct overgrowth syndrome, biochemically demonstrate activation of PI3K signaling and thereby identify a rational therapeutic target.

Rare MTNR1B variants impairing melatonin receptor 1B function contribute to type 2 diabetes.

Genome-wide association studies have revealed that common noncoding variants in MTNR1B (encoding melatonin receptor 1B, also known as MT(2)) increase type 2 diabetes (T2D) risk(1,2). Although the strongest association signal was highly significant (P < 1 × 10(-20)), its contribution to T2D risk was modest (odds ratio (OR) of ∼1.10-1.15)(1-3). We performed large-scale exon resequencing in 7,632 Europeans, including 2,186 individuals with T2D, and identified 40 nonsynonymous variants, including 36 very rare variants (minor allele frequency (MAF) <0.1%), associated with T2D (OR = 3.31, 95% confidence interval (CI) = 1.78-6.18; P = 1.64 × 10(-4)). A four-tiered functional investigation of all 40 mutants revealed that 14 were non-functional and rare (MAF < 1%), and 4 were very rare with complete loss of melatonin binding and signaling capabilities. Among the very rare variants, the partial- or total-loss-of-function variants but not the neutral ones contributed to T2D (OR = 5.67, CI = 2.17-14.82; P = 4.09 × 10(-4)). Genotyping the four complete loss-of-function variants in 11,854 additional individuals revealed their association with T2D risk (8,153 individuals with T2D and 10,100 controls; OR = 3.88, CI = 1.49-10.07; P = 5.37 × 10(-3)). This study establishes a firm functional link between MTNR1B and T2D risk.