Exome array analysis identifies new loci and low-frequency variants influencing insulin processing and secretion

TBC1D30 regulates proinsulin and insulin secretion and is the target of a genomic association signal for proinsulin

AIMS/HYPOTHESIS

Components of the insulin processing and secretion pathways remain incompletely understood. Here, we examined a genome-wide association study (GWAS) signal for plasma proinsulin levels. Lead GWAS variant rs150781447-T encodes an Arg279Cys substitution in TBC1 domain family member 30 (TBC1D30), but no role for this protein in insulin processing or secretion has been established previously. This study aimed to evaluate whether TBC1D30 drives the GWAS association signal by determining whether TBC1D30 is involved in proinsulin secretion and, if so, to examine the effects of variant alleles and potential mechanisms.

METHODS

Using CRISPR/Cas9 genome editing to create double-strand breaks and prime editing to install substitutions in INS1 832/13 insulinoma cells, we generated clonal cell lines with altered TBC1D30, as well as homozygous and heterozygous lines carrying the lead GWAS variant. We characterised lines by Sanger sequencing, quantitative PCR and ELISAs to measure glucose-stimulated proinsulin and insulin secretion. We also tested the effects of TBC1D30 knockdown on proinsulin and insulin secretion in human islets. We further assessed TBC1D30's contribution to secretory pathways by examining the effects of altered gene function on intracellular proinsulin and insulin content and insulin localisation, and by identifying potential proteins that interact with TBC1D30 using affinity purification mass spectrometry.

RESULTS

Compared with mock-edited cells, cell lines with reduced TBC1D30 expression or altered Rab GTPase-activating protein (RabGAP) domain had significantly more secreted proinsulin, 1.8- and 2.6-fold more than controls, respectively. Similarly, cells expressing the variant substitution demonstrated increased proinsulin secretion. Cell lines with a partial deletion of a critical functional domain showed 1.8-fold higher expression of Tbc1d30 and at least 2.0-fold less secreted proinsulin. Cells with altered RabGAP domain sequence also demonstrated, to a lesser extent, changes in secreted insulin levels. TBC1D30 knockdown in human islets resulted in increased insulin secretion with no significant effect on proinsulin secretion. The effects of altered TBC1D30 on mislocalisation of insulin, intracellular proinsulin and insulin content and the identities of interacting proteins are consistent with a role for TBC1D30 in proinsulin and insulin secretion.

CONCLUSIONS/INTERPRETATION

These findings suggest that effects on TBC1D30 are responsible for the GWAS signal and that TBC1D30 plays a critical role in the secretion of mature insulin.

Deficiency of Peptidylglycine-alpha-amidating Monooxygenase, a Cause of Sarcopenic Diabetes Mellitus

CONTEXT

Peptidylglycine-α-amidating monooxygenase (PAM) is a critical enzyme in the endocrine system responsible for activation, by amidation, of bioactive peptides.

OBJECTIVE

To define the clinical phenotype of carriers of genetic mutations associated with impaired PAM-amidating activity (PAM-AMA).

DESIGN

We used genetic and phenotypic data from cohort studies: the Malmö Diet and Cancer (MDC; 1991-1996; reexamination in 2002-2012), the Malmö Preventive Project (MPP; 2002-2006), and the UK Biobank (UKB; 2012).

SETTING

Exome-wide association analysis was used to identify loss-of-function (LoF) variants associated with reduced PAM-AMA and subsequently used for association with the outcomes.

PATIENTS OR OTHER PARTICIPANTS

This study included n∼4500 participants from a subcohort of the MDC (MDC-Cardiovascular cohort), n∼4500 from MPP, and n∼300,000 from UKB.

MAIN OUTCOME MEASURES

Endocrine-metabolic traits suggested by prior literature, muscle mass, muscle function, and sarcopenia.

RESULTS

Two LoF variants in the PAM gene, Ser539Trp (minor allele frequency: 0.7%) and Asp563Gly (5%), independently contributed to a decrease of 2.33 [95% confidence interval (CI): 2.52/2.15; P = 2.5E-140] and 0.98 (1.04/0.92; P = 1.12E-225) SD units of PAM-AMA, respectively. The cumulative effect of the LoF was associated with diabetes, reduced insulin secretion, and higher levels of GH and IGF-1. Moreover, carriers had reduced muscle mass and function, followed by a higher risk of sarcopenia. Indeed, the Ser539Trp mutation increased the risk of sarcopenia by 30% (odds ratio 1.31; 95% CI: 1.16/1.47; P = 9.8E-06), independently of age and diabetes.

CONCLUSION

PAM-AMA genetic deficiency results in a prediabetic sarcopenic phenotype. Early identification of PAM LoF carriers would allow targeted exercise interventions and calls for novel therapies that restore enzymatic activity.

Insulin-related traits and prostate cancer: A Mendelian randomization study

Investigating the causal relationship between insulin secretion and prostate cancer (PCa) development is challenging due to the multifactorial nature of PCa, which complicates the isolation of the specific impact of insulin-related factors. We conducted a Mendelian randomization (MR) study to investigate the associations between insulin secretion-related traits and PCa. We used 36, 60, 56, 23, 48, and 49 single nucleotide polymorphisms (SNPs) as instrumental variables for fasting insulin, insulin sensitivity, proinsulin, and proinsulin in nondiabetic individuals, individuals with diabetes, and individuals receiving exogenous insulin, respectively. These SNPs were selected from various genome-wide association studies. To clarify the causal relationship between insulin-related traits and PCa, we utilized a multivariable MR analysis to adjust for obesity and body fat percentage. Additionally, two-step Mendelian randomization was conducted to assess the role of insulin-like growth factor 1 (IGF-1) in the relationship between proinsulin and PCa. Two-sample MR analysis revealed strong associations between genetically predicted fasting insulin, insulin sensitivity, proinsulin, and proinsulin in nondiabetic individuals and the development of PCa. After adjustment for obesity and body fat percentage using multivariable MR analysis, proinsulin remained significantly associated with PCa, whereas other factors were not. Furthermore, two-step MR analysis demonstrated that proinsulin acts as a negative factor in prostate carcinogenesis, largely independent of IGF-1. This study provides evidence suggesting that proinsulin may act as a negative factor contributing to the development of PCa. Novel therapies targeting proinsulin may have potential benefits for PCa patients, potentially reducing the need for unnecessary surgical treatments.

Histone lactylation mediated by Fam172a in POMC neurons regulates energy balance

Glycolysis-derived lactate was identified as substrate for histone lactylation, which has been regarded as a significant role in transcriptional regulation in many tissues. However, the role of histone lactylation in the metabolic center, the hypothalamus, is still unknown. Here, we show that hypothalamic pro-opiomelanocortin (POMC) neuron-specific deletion of family with sequence similarity 172, member A (Fam172a) can increase histone lactylation and protect mice against diet-induced obesity (DIO) and related metabolic disorders. Conversely, overexpression of Fam172a in POMC neurons led to an obesity-like phenotype. Using RNA-seq and CUT&Tag chromatin profiling analyzes, we find that knockdown of Fam172a activates the glycolytic process and increases peptidylglycine α-amidating monooxygenase (PAM), which affects the synthesis of α-MSH, via H4K12la (histone lactylation). In addition, pharmacological inhibition of lactate production clearly abrogates the anti-obesity effect of PFKO (POMC-Cre, Fam172aloxP/loxP, POMC neurons Fam172a knockout). These findings highlight the importance of Fam172a and lactate in the development of obesity and our results mainly concern male mice.

Multi-omic human pancreatic islet endoplasmic reticulum and cytokine stress response mapping provides type 2 diabetes genetic insights

Endoplasmic reticulum (ER) and inflammatory stress responses contribute to islet dysfunction in type 2 diabetes (T2D). Comprehensive genomic understanding of these human islet stress responses and whether T2D-associated genetic variants modulate them is lacking. Here, comparative transcriptome and epigenome analyses of human islets exposed ex vivo to these stressors revealed 30% of expressed genes and 14% of islet cis-regulatory elements (CREs) as stress responsive, modulated largely in an ER- or cytokine-specific fashion. T2D variants overlapped 86 stress-responsive CREs, including 21 induced by ER stress. We linked the rs6917676-T T2D risk allele to increased islet ER-stress-responsive CRE accessibility and allele-specific β cell nuclear factor binding. MAP3K5, the ER-stress-responsive putative rs6917676 T2D effector gene, promoted stress-induced β cell apoptosis. Supporting its pro-diabetogenic role, MAP3K5 expression correlated inversely with human islet β cell abundance and was elevated in T2D β cells. This study provides genome-wide insights into human islet stress responses and context-specific T2D variant effects.

Genetic architecture of oral glucose-stimulated insulin release provides biological insights into type 2 diabetes aetiology

The genetics of β-cell function (BCF) offer valuable insights into the aetiology of type 2 diabetes (T2D)1,2. Previous studies have expanded the catalogue of BCF genetic associations through candidate gene studies3-7, large-scale genome-wide association studies (GWAS) of fasting BCF8,9 or functional islet studies on T2D risk variants10-14. Nonetheless, GWAS focused on BCF traits derived from oral glucose tolerance test (OGTT) data have been limited in sample size15,16 and have often overlooked the potential for related traits to capture distinct genetic features of insulin-producing β-cells17,18. We reasoned that investigating the genetic basis of multiple BCF estimates could provide a broader understanding of β-cell physiology. Here, we aggregate GWAS data of eight OGTT-based BCF traits from ~26,000 individuals of European descent, identifying 55 independent genetic associations at 44 loci. By examining the effects of BCF genetic signals on related phenotypes, we uncover diverse disease mechanisms whereby genetic regulation of BCF may influence T2D risk. Integrating BCF-GWAS data with pancreatic islet transcriptomic and epigenomic datasets reveals 92 candidate effector genes. Gene silencing in β-cell models highlights ACSL1 and FAM46C as key regulators of insulin secretion. Overall, our findings yield insights into the biology of insulin release and the molecular processes linking BCF to T2D risk, shedding light on the heterogeneity of T2D pathophysiology.

cNPAS2 induced β cell dysfunction by regulating KANK1 expression in type 2 diabetes

BACKGROUND

Diabetes mellitus type 2 (T2DM) is formed by defective insulin secretion with the addition of peripheral tissue resistance of insulin action. It has been affecting over 400 million people all over the world.

AIM

To explore the pathogenesis of T2DM and to develop and implement new prevention and treatment strategies for T2DM.

METHODS

Receiver operating characteristic (ROC) curve analysis was used to conduct diagnostic markers. The expression level of genes was determined by reverse transcription-PCR as well as Western blot. Cell proliferation assays were performed by cell counting kit-8 (CCK-8) tests. At last, T2DM mice underwent Roux-en-Y gastric bypass surgery.

RESULTS

We found that NPAS2 was significantly up-regulated in islet β cell apoptosis of T2DM. The ROC curve revealed that NPAS2 was capable of accurately diagnosing T2DM. NPAS2 overexpression did increase the level of KANK1. In addition, the CCK-8 test revealed knocking down NPAS2 and KANK1 increased the proliferation of MIN6 cells. At last, we found that gastric bypass may treat type 2 diabetes by down-regulating NPAS2 and KANK1.

CONCLUSION

This study demonstrated that NPAS2 induced β cell dysfunction by regulating KANK1 expression in type 2 diabetes, and it may be an underlying therapy target of T2DM.

Characterizing Common and Rare Variations in Nontraditional Glycemic Biomarkers Using Multivariate Approaches on Multiancestry ARIC Study

Genetic studies of nontraditional glycemic biomarkers, glycated albumin and fructosamine, can shed light on unknown aspects of type 2 diabetes genetics and biology. We performed a multiphenotype genome-wide association study of glycated albumin and fructosamine from 7,395 White and 2,016 Black participants in the Atherosclerosis Risk in Communities (ARIC) study on common variants from genotyped/imputed data. We discovered two genome-wide significant loci, one mapping to a known type 2 diabetes gene (ARAP1/STARD10) and another mapping to a novel region (UGT1A complex of genes), using multiomics gene-mapping strategies in diabetes-relevant tissues. We identified additional loci that were ancestry- and sex-specific (e.g., PRKCA in African ancestry, FCGRT in European ancestry, TEX29 in males). Further, we implemented multiphenotype gene-burden tests on whole-exome sequence data from 6,590 White and 2,309 Black ARIC participants. Ten variant sets annotated to genes across different variant aggregation strategies were exome-wide significant only in multiancestry analysis, of which CD1D, EGFL7/AGPAT2, and MIR126 had notable enrichment of rare predicted loss of function variants in African ancestry despite smaller sample sizes. Overall, 8 of 14 discovered loci and genes were implicated to influence these biomarkers via glycemic pathways, and most of them were not previously implicated in studies of type 2 diabetes. This study illustrates improved locus discovery and potential effector gene discovery by leveraging joint patterns of related biomarkers across the entire allele frequency spectrum in multiancestry analysis. Future investigation of the loci and genes potentially acting through glycemic pathways may help us better understand the risk of developing type 2 diabetes.

ARTICLE HIGHLIGHTS

A splice site variant in MADD affects hormone expression in pancreatic β cells and pituitary gonadotropes

MAPK activating death domain (MADD) is a multifunctional protein regulating small GTPases RAB3 and RAB27, MAPK signaling, and cell survival. Polymorphisms in the MADD locus are associated with glycemic traits, but patients with biallelic variants in MADD manifest a complex syndrome affecting nervous, endocrine, exocrine, and hematological systems. We identified a homozygous splice site variant in MADD in 2 siblings with developmental delay, diabetes, congenital hypogonadotropic hypogonadism, and growth hormone deficiency. This variant led to skipping of exon 30 and in-frame deletion of 36 amino acids. To elucidate how this mutation causes pleiotropic endocrine phenotypes, we generated relevant cellular models with deletion of MADD exon 30 (dex30). We observed reduced numbers of β cells, decreased insulin content, and increased proinsulin-to-insulin ratio in dex30 human embryonic stem cell-derived pancreatic islets. Concordantly, dex30 led to decreased insulin expression in human β cell line EndoC-βH1. Furthermore, dex30 resulted in decreased luteinizing hormone expression in mouse pituitary gonadotrope cell line LβT2 but did not affect ontogeny of stem cell-derived GnRH neurons. Protein-protein interactions of wild-type and dex30 MADD revealed changes affecting multiple signaling pathways, while the GDP/GTP exchange activity of dex30 MADD remained intact. Our results suggest MADD-specific processes regulate hormone expression in pancreatic β cells and pituitary gonadotropes.

Unravelling the Genetic Landscape of Hemiplegic Migraine: Exploring Innovative Strategies and Emerging Approaches

Migraine is a severe, debilitating neurovascular disorder. Hemiplegic migraine (HM) is a rare and debilitating neurological condition with a strong genetic basis. Sequencing technologies have improved the diagnosis and our understanding of the molecular pathophysiology of HM. Linkage analysis and sequencing studies in HM families have identified pathogenic variants in ion channels and related genes, including CACNA1A, ATP1A2, and SCN1A, that cause HM. However, approximately 75% of HM patients are negative for these mutations, indicating there are other genes involved in disease causation. In this review, we explored our current understanding of the genetics of HM. The evidence presented herein summarises the current knowledge of the genetics of HM, which can be expanded further to explain the remaining heritability of this debilitating condition. Innovative bioinformatics and computational strategies to cover the entire genetic spectrum of HM are also discussed in this review.

Disentangling archaic introgression and genomic signatures of selection at human immunity genes

Pathogens and infectious diseases have imposed exceptionally strong selective pressure on ancient and modern human genomes and contributed to the current variation in many genes. There is evidence that modern humans acquired immune variants through interbreeding with ancient hominins, but the impact of such variants on human traits is not fully understood. The main objectives of this research were to infer the genetic signatures of positive selection that may be involved in adaptation to infectious diseases and to investigate the function of Neanderthal alleles identified within a set of 50 Lithuanian genomes. Introgressed regions were identified using the machine learning tool ArchIE. Recent positive selection signatures were analysed using iHS. We detected high-scoring signals of positive selection at innate immunity genes (EMB, PARP8, HLAC, and CDSN) and evaluated their interactions with the structural proteins of pathogens. Interactions with human immunodeficiency virus (HIV) 1 and severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) were identified. Overall, genomic regions introgressed from Neanderthals were shown to be enriched in genes related to immunity, keratinocyte differentiation, and sensory perception.

Characterizing common and rare variations in non-traditional glycemic biomarkers using multivariate approaches on multi-ancestry ARIC study

Glycated hemoglobin, fasting glucose, glycated albumin, and fructosamine are biomarkers that reflect different aspects of the glycemic process. Genetic studies of these glycemic biomarkers can shed light on unknown aspects of type 2 diabetes genetics and biology. While there exists several GWAS of glycated hemoglobin and fasting glucose, very few GWAS have focused on glycated albumin or fructosamine. We performed a multi-phenotype GWAS of glycated albumin and fructosamine from 7,395 White and 2,016 Black participants in the Atherosclerosis Risk in Communities (ARIC) study on the common variants from genotyped/imputed data. We found 2 genome-wide significant loci, one mapping to known type 2 diabetes gene (ARAP1/STARD10, p = 2.8 × 10-8) and another mapping to a novel gene (UGT1A, p = 1.4 × 10-8) using multi-omics gene mapping strategies in diabetes-relevant tissues. We identified additional loci that were ancestry-specific (e.g., PRKCA from African ancestry individuals, p = 1.7 × 10-8) and sex-specific (TEX29 locus in males only, p = 3.0 × 10-8). Further, we implemented multi-phenotype gene-burden tests on whole-exome sequence data from 6,590 White and 2,309 Black ARIC participants. Eleven genes across different rare variant aggregation strategies were exome-wide significant only in multi-ancestry analysis. Four out of 11 genes had notable enrichment of rare predicted loss of function variants in African ancestry participants despite smaller sample size. Overall, 8 out of 15 loci/genes were implicated to influence these biomarkers via glycemic pathways. This study illustrates improved locus discovery and potential effector gene discovery by leveraging joint patterns of related biomarkers across entire allele frequency spectrum in multi-ancestry analyses. Most of the loci/genes we identified have not been previously implicated in studies of type 2 diabetes, and future investigation of the loci/genes potentially acting through glycemic pathways may help us better understand risk of developing type 2 diabetes.

Genetics of diabetes

Diabetes mellitus is a complicated disease characterized by a complex interplay of genetic, epigenetic, and environmental variables. It is one of the world's fastest-growing diseases, with 783 million adults expected to be affected by 2045. Devastating macrovascular consequences (cerebrovascular disease, cardiovascular disease, and peripheral vascular disease) and microvascular complications (like retinopathy, nephropathy, and neuropathy) increase mortality, blindness, kidney failure, and overall quality of life in individuals with diabetes. Clinical risk factors and glycemic management alone cannot predict the development of vascular problems; multiple genetic investigations have revealed a clear hereditary component to both diabetes and its related complications. In the twenty-first century, technological advancements (genome-wide association studies, next-generation sequencing, and exome-sequencing) have led to the identification of genetic variants associated with diabetes, however, these variants can only explain a small proportion of the total heritability of the condition. In this review, we address some of the likely explanations for this "missing heritability", for diabetes such as the significance of uncommon variants, gene-environment interactions, and epigenetics. Current discoveries clinical value, management of diabetes, and future research directions are also discussed.

Type 2 Diabetes risk alleles in Peptidyl-glycine Alpha-amidating Monooxygenase influence GLP-1 levels and response to GLP-1 Receptor Agonists

Patients with type 2 diabetes vary in their response to currently available therapeutic agents (including GLP-1 receptor agonists) leading to suboptimal glycemic control and increased risk of complications. We show that human carriers of hypomorphic T2D-risk alleles in the gene encoding peptidyl-glycine alpha-amidating monooxygenase (PAM), as well as Pam-knockout mice, display increased resistance to GLP-1 in vivo. Pam inactivation in mice leads to reduced gastric GLP-1R expression and faster gastric emptying: this persists during GLP-1R agonist treatment and is rescued when GLP-1R activity is antagonized, indicating resistance to GLP-1's gastric slowing properties. Meta-analysis of human data from studies examining GLP-1R agonist response (including RCTs) reveals a relative loss of 44% and 20% of glucose lowering (measured by glycated hemoglobin) in individuals with hypomorphic PAM alleles p.S539W and p.D536G treated with GLP-1R agonist. Genetic variation in PAM has effects on incretin signaling that alters response to medication used commonly for treatment of T2D.

The potential of integrating human and mouse discovery platforms to advance our understanding of cardiometabolic diseases

Cardiometabolic diseases encompass a range of interrelated conditions that arise from underlying metabolic perturbations precipitated by genetic, environmental, and lifestyle factors. While obesity, dyslipidaemia, smoking, and insulin resistance are major risk factors for cardiometabolic diseases, individuals still present in the absence of such traditional risk factors, making it difficult to determine those at greatest risk of disease. Thus, it is crucial to elucidate the genetic, environmental, and molecular underpinnings to better understand, diagnose, and treat cardiometabolic diseases. Much of this information can be garnered using systems genetics, which takes population-based approaches to investigate how genetic variance contributes to complex traits. Despite the important advances made by human genome-wide association studies (GWAS) in this space, corroboration of these findings has been hampered by limitations including the inability to control environmental influence, limited access to pertinent metabolic tissues, and often, poor classification of diseases or phenotypes. A complementary approach to human GWAS is the utilisation of model systems such as genetically diverse mouse panels to study natural genetic and phenotypic variation in a controlled environment. Here, we review mouse genetic reference panels and the opportunities they provide for the study of cardiometabolic diseases and related traits. We discuss how the post-GWAS era has prompted a shift in focus from discovery of novel genetic variants to understanding gene function. Finally, we highlight key advantages and challenges of integrating complementary genetic and multi-omics data from human and mouse populations to advance biological discovery.

Loci for insulin processing and secretion provide insight into type 2 diabetes risk

Insulin secretion is critical for glucose homeostasis, and increased levels of the precursor proinsulin relative to insulin indicate pancreatic islet beta-cell stress and insufficient insulin secretory capacity in the setting of insulin resistance. We conducted meta-analyses of genome-wide association results for fasting proinsulin from 16 European-ancestry studies in 45,861 individuals. We found 36 independent signals at 30 loci (p value < 5 × 10-8), which validated 12 previously reported loci for proinsulin and ten additional loci previously identified for another glycemic trait. Half of the alleles associated with higher proinsulin showed higher rather than lower effects on glucose levels, corresponding to different mechanisms. Proinsulin loci included genes that affect prohormone convertases, beta-cell dysfunction, vesicle trafficking, beta-cell transcriptional regulation, and lysosomes/autophagy processes. We colocalized 11 proinsulin signals with islet expression quantitative trait locus (eQTL) data, suggesting candidate genes, including ARSG, WIPI1, SLC7A14, and SIX3. The NKX6-3/ANK1 proinsulin signal colocalized with a T2D signal and an adipose ANK1 eQTL signal but not the islet NKX6-3 eQTL. Signals were enriched for islet enhancers, and we showed a plausible islet regulatory mechanism for the lead signal in the MADD locus. These results show how detailed genetic studies of an intermediate phenotype can elucidate mechanisms that may predispose one to disease.

Genome-wide association studies for thyroid physiology and diseases

Genetic factors are involved in the etiology of most diseases, but prior to 2000, the methods for identifying such factors were very limited. Genome-wide association study (GWAS), developed in the 2000s, is an analytical method that can be applied to most diseases, including endocrine disorders. GWAS has provided a wealth of information on disease risks and the molecular pathogenesis of many human diseases. This review summarizes key findings from GWAS for thyroid physiology and diseases, and illustrates how GWAS is a powerful research tool to elucidate the molecular mechanisms of the diseases.

The SFT2D2 gene is associated with the autoimmune pathology of schizophrenia in a Chinese population

Background

The relative risk of GWAS-confirmed loci strongly associated with schizophrenia may be underestimated due to the decay of linkage disequilibrium between index SNPs and causal variants. This study is aimed to investigate schizophrenia-associated signals detected in the 1q24-25 region in order to identify a causal variant in LD with GWAS index SNPs, and the potential biological functions of the risk gene.

Methods

Re-genotyping analysis was performed in the 1q24-25 region that harbors three GWAS index SNPs associated with schizophrenia (rs10489202, rs11586522, and rs6670165) in total of 9801 case-control subjects of Chinese Han origin. Circulating autoantibody levels were assessed using an in-house ELISA against a protein derived fragment encoded by SFT2D2 in total of 682 plasma samples.

Results

A rare variant (rs532193193) in the SFT2D2 locus was identified to be strongly associated with schizophrenia. Compared with control subjects, patients with schizophrenia showed increased anti-SFT2D2 IgG levels. Receiver operating characteristic (ROC) analysis revealed an area under the ROC curve (AUC) of 0.803 with sensitivity of 28.57% against specificity of 95% for the anti-SFT2D2 IgG assay.

Discussion

Our findings indicate that SFT2D2 is a novel gene for risk of schizophrenia, while endogenous anti-SFT2D2 IgG may underlie the pathophysiology of the immunological aspects of schizophrenia.

The contribution of common and rare genetic variants to variation in metabolic traits in 288,137 East Asians

Metabolic traits are heritable phenotypes widely-used in assessing the risk of various diseases. We conduct a genome-wide association analysis (GWAS) of nine metabolic traits (including glycemic, lipid, liver enzyme levels) in 125,872 Korean subjects genotyped with the Korea Biobank Array. Following meta-analysis with GWAS from Biobank Japan identify 144 novel signals (MAF ≥ 1%), of which 57.0% are replicated in UK Biobank. Additionally, we discover 66 rare (MAF < 1%) variants, 94.4% of them co-incident to common loci, adding to allelic series. Although rare variants have limited contribution to overall trait variance, these lead, in carriers, substantial loss of predictive accuracy from polygenic predictions of disease risk from common variant alone. We capture groups with up to 16-fold variation in type 2 diabetes (T2D) prevalence by integration of genetic risk scores of fasting plasma glucose and T2D and the I349F rare protective variant. This study highlights the need to consider the joint contribution of both common and rare variants on inherited risk of metabolic traits and related diseases.

The Genetic Architecture of the Etiology of Lower Extremity Peripheral Artery Disease: Current Knowledge and Future Challenges in the Era of Genomic Medicine

Lower extremity artery disease (LEAD), caused by atherosclerotic obstruction of the arteries of the lower limb extremities, has exhibited an increase in mortality and morbidity worldwide. The phenotypic variability of LEAD is correlated with its complex, multifactorial etiology. In addition to traditional risk factors, it has been shown that the interaction between genetic factors (epistasis) or between genes and the environment potentially have an independent role in the development and progression of LEAD. In recent years, progress has been made in identifying genetic variants associated with LEAD, by Genome-Wide Association Studies (GWAS), Whole Exome Sequencing (WES) studies, and epigenetic profiling. The aim of this review is to present the current knowledge about the genetic factors involved in the etiopathogenic mechanisms of LEAD, as well as possible directions for future research. We analyzed data from the literature, starting with candidate gene-based association studies, and then continuing with extensive association studies, such as GWAS and WES. The results of these studies showed that the genetic architecture of LEAD is extremely heterogeneous. In the future, the identification of new genetic factors will allow for the development of targeted molecular therapies, and the use of polygenic risk scores (PRS) to identify individuals at an increased risk of LEAD will allow for early prophylactic measures and personalized therapy to improve their prognosis.

Long non-coding RNAs: a valuable biomarker for metabolic syndrome

Long non-coding RNAs (lncRNAs) have become important regulators of gene expression because they affect a wide range of biological processes, such as cell growth, death, differentiation, and aging. More and more evidence suggests that lncRNAs play a role in maintaining metabolic homeostasis. When certain lncRNAs are out of balance, metabolic disorders like diabetes, obesity, and heart disease get worse. In this review, we talk about what we know about how lncRNAs control metabolism, with a focus on diseases caused by long-term inflammation and the characteristics of the metabolic syndrome. We looked at lncRNAs and their molecular targets in the pathogenesis of signaling pathways. We also talked about how lncRNAs are becoming more and more interesting as diagnostic and therapeutic targets for improving metabolic homeostasis.

Genetics of Type 2 Diabetes: Past, Present, and Future

Diabetes has reached epidemic proportions worldwide. Currently, approximately 537 million adults (20–79 years) have diabetes, and the total number of people with diabetes is continuously increasing. Diabetes includes several subtypes. About 80% of all cases of diabetes are type 2 diabetes (T2D). T2D is a polygenic disease with an inheritance ranging from 30 to 70%. Genetic and environment/lifestyle factors, especially obesity and sedentary lifestyle, increase the risk of T2D. In this review, we discuss how studies on the genetics of diabetes started, how they expanded when genome-wide association studies and exome and whole-genome sequencing became available, and the current challenges in genetic studies of diabetes. T2D is heterogeneous with respect to clinical presentation, disease course, and response to treatment, and has several subgroups which differ in pathophysiology and risk of micro- and macrovascular complications. Currently, genetic studies of T2D focus on these subgroups to find the best diagnoses and treatments for these patients according to the principles of precision medicine.

PAM variants were associated with type 2 diabetes mellitus risk in the Chinese population

This study aimed to assess the association between PAM single-nucleotide polymorphisms (SNPs) and T2DM risk in the Chinese population. We performed the genotype of PAM SNPs using Agena MassARRAY in 1002 subjects. The effect of PAM polymorphisms on T2DM occurrence was evaluated by logistic regression analysis. False-positive report probability (FPRP) was utilized to assess the noteworthiness of the significant results. This study showed that PAM rs406761, rs17154889, and rs6889592 were related to an increased risk of T2DM. The similar results were also in subjects with ≤ 60 years. Rs2431320 and rs406761 were related to an increased risk of T2DM in males, and rs6889592 was only found to be associated with T2DM risk in females. Rs2431320 and rs406761 increased T2DM risk in people with BMI > 24, and rs6889592 and rs26431 significantly correlated with T2DM risk in people with BMI ≤ 24. By comparing patients with no retinopathy with controls, the correlation between PAM rs406761 and rs17154889 and T2DM risk was observed. The significant association between T2DM risk and PAM SNPs was remarkable by FPRP values. PAM SNPs were correlated with T2DM risk in the Chinese population, illustrating the importance of PAM SNPs in the pathogenesis of T2DM.

Integration of Transcriptome and Exome Genotyping Identifies Significant Variants with Autism Spectrum Disorder

Autism is a complex disease with genetic predisposition factors. Real factors for treatment and early diagnosis are yet to be defined. This study integrated transcriptome and exome genotyping for identifying functional variants associated with autism spectrum disorder and their impact on gene expression to find significant variations. More than 1800 patients were screened, and 70 (47 male/23 female) with an average age of 7.56 ± 3.68 years fulfilled the DSM-5 criteria for autism. Analysis revealed 682 SNPs of 589 genes significantly (p < 0.001) associated with autism among the putative functional exonic variants (n = 243,345) studied. Olfactory receptor genes on chromosome 6 were significant after Bonferroni correction (α = 0.05/243345 = 2.05 × 10−7) with a high degree of linkage disequilibrium on 6p22.1 (p = 6.71 × 10−9). The differentially expressed gene analysis of autistic patients compared to controls in whole RNA sequencing identified significantly upregulated (foldchange ≥ 0.8 and p-value ≤ 0.05; n = 125) and downregulated (foldchange ≤ −0.8 and p-value ≤ 0.05; n = 117) genes. The integration of significantly up- and downregulated genes and genes of significant SNPs identified regulatory variants (rs6657480, rs3130780, and rs1940475) associated with the up- (ITGB3BP) and downregulation (DDR1 and MMP8) of genes in autism spectrum disorder in people of Arab ancestries. The significant variants could be a biomarker of interest for identifying early autism among Arabs and helping to characterize the genes involved in the susceptibility mechanisms for autistic subjects.

Coding and regulatory variants are associated with serum protein levels and disease

Circulating proteins can be used to diagnose and predict disease-related outcomes. A deep serum proteome survey recently revealed close associations between serum protein networks and common disease. In the current study, 54,469 low-frequency and common exome-array variants were compared to 4782 protein measurements in the serum of 5343 individuals from the AGES Reykjavik cohort. This analysis identifies a large number of serum proteins with genetic signatures overlapping those of many diseases. More specifically, using a study-wide significance threshold, we find that 2021 independent exome array variants are associated with serum levels of 1942 proteins. These variants reside in genetic loci shared by hundreds of complex disease traits, highlighting serum proteins' emerging role as biomarkers and potential causative agents of a wide range of diseases.

Subcutaneous adipose tissue splice quantitative trait loci reveal differences in isoform usage associated with cardiometabolic traits

Alternate splicing events can create isoforms that alter gene function, and genetic variants associated with alternate gene isoforms may reveal molecular mechanisms of disease. We used subcutaneous adipose tissue of 426 Finnish men from the METSIM study and identified splice junction quantitative trait loci (sQTLs) for 6,077 splice junctions (FDR < 1%). In the same individuals, we detected expression QTLs (eQTLs) for 59,443 exons and 15,397 genes (FDR < 1%). We identified 595 genes with an sQTL and exon eQTL but no gene eQTL, which could indicate potential isoform differences. Of the significant sQTL signals, 2,114 (39.8%) included at least one proxy variant (linkage disequilibrium r2 > 0.8) located within an intron spanned by the splice junction. We identified 203 sQTLs that colocalized with 141 genome-wide association study (GWAS) signals for cardiometabolic traits, including 25 signals for lipid traits, 24 signals for body mass index (BMI), and 12 signals for waist-hip ratio adjusted for BMI. Among all 141 GWAS signals colocalized with an sQTL, we detected 26 that also colocalized with an exon eQTL for an exon skipped by the sQTL splice junction. At a GWAS signal for high-density lipoprotein cholesterol colocalized with an NR1H3 sQTL splice junction, we show that the alternative splice product encodes an NR1H3 transcription factor that lacks a DNA binding domain and fails to activate transcription. Together, these results detect splicing events and candidate mechanisms that may contribute to gene function at GWAS loci.

Stem Cell-Derived β Cells: A Versatile Research Platform to Interrogate the Genetic Basis of β Cell Dysfunction

Pancreatic β cell dysfunction is a central component of diabetes progression. During the last decades, the genetic basis of several monogenic forms of diabetes has been recognized. Genome-wide association studies (GWAS) have also facilitated the identification of common genetic variants associated with an increased risk of diabetes. These studies highlight the importance of impaired β cell function in all forms of diabetes. However, how most of these risk variants confer disease risk, remains unanswered. Understanding the specific contribution of genetic variants and the precise role of their molecular effectors is the next step toward developing treatments that target β cell dysfunction in the era of personalized medicine. Protocols that allow derivation of β cells from pluripotent stem cells, represent a powerful research tool that allows modeling of human development and versatile experimental designs that can be used to shed some light on diabetes pathophysiology. This article reviews different models to study the genetic basis of β cell dysfunction, focusing on the recent advances made possible by stem cell applications in the field of diabetes research.

Novel insights into peptide amidation and amidating activity in the human circulation

C-terminal α-amidation is the final and essential step in the biosynthesis of several peptide hormones. Peptidylglycine α-amidating monooxygenase (PAM) is the only known enzyme to catalyse this reaction. PAM amidating activity (AMA) is known to be present in human circulation, but its physiological role and significance as a clinical biomarker remains unclear. We developed a PAM-specific amidation assay that utilizes the naturally occurring substrate Adrenomedullin-Gly (ADM-Gly, 1-53). Using our amidation assay we quantified serum amidating activities in a large population-based cohort of more than 4900 individuals. A correlation of serum amidating activity with several clinical parameters including high blood pressure was observed. Increasing PAM-AMA was an independent predictor of hard outcomes related to hemodynamic stress such as cardiovascular mortality, atrial fibrillation and heart failure during long-term follow-up (8.8 ± 2.5 years). Moreover, results from an animal study in rats utilizing recombinant human PAM provide novel insights into the physiological role of circulating PAM and show its potential significance in circulating peptide amidation.

Genetics of Type 2 Diabetes: Opportunities for Precision Medicine: JACC Focus Seminar

Type 2 diabetes (T2D) is highly prevalent and is a strong contributor for cardiovascular disease. However, there is significant heterogeneity in disease pathogenesis and the risk of complications. Enormous progress has been made in our ability to catalog genetic variation associated with T2D risk and variation in disease-relevant quantitative traits. These discoveries hold the potential to shed light on tractable targets and pathways for safe and effective therapeutic development, but the promise of precision medicine has been slow to be realized. Recent studies have identified subgroups of individuals with differential risk for intermediate phenotypes (eg, lipid levels, fasting insulin, body mass index) that contribute to T2D risk, helping to account for the observed clinical heterogeneity. These "partitioned genetic risk scores" not only have the potential to identify patients at greatest risk of cardiovascular disease and rapid disease progression, but also could aid patient stratification bridging the gap toward precision medicine for T2D.

GIGYF1 loss of function is associated with clonal mosaicism and adverse metabolic health

Mosaic loss of chromosome Y (LOY) in leukocytes is the most common form of clonal mosaicism, caused by dysregulation in cell-cycle and DNA damage response pathways. Previous genetic studies have focussed on identifying common variants associated with LOY, which we now extend to rarer, protein-coding variation using exome sequences from 82,277 male UK Biobank participants. We find that loss of function of two genes—CHEK2 and GIGYF1—reach exome-wide significance. Rare alleles in GIGYF1 have not previously been implicated in any complex trait, but here loss-of-function carriers exhibit six-fold higher susceptibility to LOY (OR = 5.99 [3.04–11.81], p = 1.3 × 10⁻¹⁰). These same alleles are also associated with adverse metabolic health, including higher susceptibility to Type 2 Diabetes (OR = 6.10 [3.51–10.61], p = 1.8 × 10⁻¹²), 4 kg higher fat mass (p = 1.3 × 10⁻⁴), 2.32 nmol/L lower serum IGF1 levels (p = 1.5 × 10⁻⁴) and 4.5 kg lower handgrip strength (p = 4.7 × 10⁻⁷) consistent with proposed GIGYF1 enhancement of insulin and IGF-1 receptor signalling. These associations are mirrored by a common variant nearby associated with the expression of GIGYF1. Our observations highlight a potential direct connection between clonal mosaicism and metabolic health.

Mosaic loss of chromosome Y (LOY) is a common form of clonal mosaicism in leukocytes. Here, the authors extend genetic association analyses to rare variation using exome-sequence data from 82,277 males, finding that loss-of-function alleles in GIGYF1 are associated with six-fold higher susceptibility to both LOY and Type 2 Diabetes.

Emerging roles of rare and low-frequency genetic variants in type 1 diabetes mellitus

Type 1 diabetes mellitus (T1DM) is defined as an autoimmune disorder and has enormous complexity and heterogeneity. Although its precise pathogenic mechanisms are obscure, this disease is widely acknowledged to be precipitated by environmental factors in individuals with genetic susceptibility. To date, the known susceptibility loci, which have mostly been identified by genome-wide association studies, can explain 80%–85% of the heritability of T1DM. Researchers believe that at least a part of its missing genetic component is caused by undetected rare and low-frequency variants. Most common variants have only small to modest effect sizes, which increases the difficulty of dissecting their functions and restricts their potential clinical application. Intriguingly, many studies have indicated that rare and low-frequency variants have larger effect sizes and play more significant roles in susceptibility to common diseases, including T1DM, than common variants do. Therefore, better recognition of rare and low-frequency variants is beneficial for revealing the genetic architecture of T1DM and for providing new and potent therapeutic targets for this disease. Here, we will discuss existing challenges as well as the great significance of this field and review current knowledge of the contributions of rare and low-frequency variants to T1DM.

Genetic variants associated with expression of TCF19 contribute to the risk of head and neck cancer in Chinese population

BACKGROUND

Squamous cell carcinoma of the head and neck (SCCHN) is one of the most common cancers worldwide and includes cancers arising from the oral cavity, pharynx and larynx. Genome-wide association studies have found several genetic variants related to the risk of SCCHN; however, they could only explain a small fraction of the heritability. Thus, more susceptibility loci associated with SCCHN need to be identified.

METHODS

An association study was conducted by genotyping 555 patients with SCCHN and 1367 controls in a Chinese population. Single-variant association analysis was conducted on 63 373 SNPs, and the promising variants were then confirmed by a two-stage validation with 1875 SCCHN cases and 4637 controls. Bioinformatics analysis and functional assays were applied to uncover the potential pathogenic mechanism of the promising variants and genes associated with SCCHN.

RESULTS

We first identified three novel genetic variants significantly associated with the risk of SCCHN (p=7.45×10^-7 for rs2517611 at 6p22.1, p=1.76×10^-9 for rs2524182 at 6p21.33 and p=2.17×10^-10 for rs3131018 at 6p21.33). Further analysis and biochemical assays showed that rs3094187, which was in a region in high linkage disequilibrium with rs3131018, could modify TCF19 expression by regulating the binding affinity of the transcription factor SREBF1 to the promoter of TCF19. In addition, experiments revealed that the inhibition of TCF19 may affect several important pathways involved in tumourigenesis and attenuate the cell proliferation and migration of SCCHN.

CONCLUSION

These findings offer important evidence that functional genetic variants could contribute to development of SCCHN and that TCF19 may function as a putative susceptibility gene for SCCHN.

Exome Chip Analysis of 14,026 Koreans Reveals Known and Newly Discovered Genetic Loci Associated with Type 2 Diabetes Mellitus

BACKGROUND

Most loci associated with type 2 diabetes mellitus (T2DM) discovered to date are within noncoding regions of unknown functional significance. By contrast, exonic regions have advantages for biological interpretation.

METHODS

We analyzed the association of exome array data from 14,026 Koreans to identify susceptible exonic loci for T2DM. We used genotype information of 50,543 variants using the Illumina exome array platform.

RESULTS

In total, 7 loci were significant with a Bonferroni adjusted P=1.03×10-6. rs2233580 in paired box gene 4 (PAX4) showed the highest odds ratio of 1.48 (P=1.60×10-10). rs11960799 in membrane associated ring-CH-type finger 3 (MARCH3) and rs75680863 in transcobalamin 2 (TCN2) were newly identified loci. When we built a model to predict the incidence of diabetes with the 7 loci and clinical variables, area under the curve (AUC) of the model improved significantly (AUC=0.72, P<0.05), but marginally in its magnitude, compared with the model using clinical variables (AUC=0.71, P<0.05). When we divided the entire population into three groups-normal body mass index (BMI; <25 kg/m2), overweight (25≤ BMI <30 kg/m2), and obese (BMI ≥30 kg/m2) individuals-the predictive performance of the 7 loci was greatest in the group of obese individuals, where the net reclassification improvement was highly significant (0.51; P=8.00×10-5).

CONCLUSION

We found exonic loci having a susceptibility for T2DM. We found that such genetic information is advantageous for predicting T2DM in a subgroup of obese individuals.

The transcriptome-wide association search for genes and genetic variants which associate with BMI and gestational weight gain in women with type 1 diabetes

BACKGROUND

Clinical data suggest that BMI and gestational weight gain (GWG) are strongly interconnected phenotypes; however, the genetic basis of the latter is rather unclear. Here we aim to find genes and genetic variants which influence BMI and/or GWG.

METHODS

We have genotyped 316 type 1 diabetics using Illumina Infinium Omni Express Exome-8 v1.4 arrays. The GIANT, ARIC and T2D-GENES summary statistics were used for TWAS (performed with PrediXcan) in adipose tissue. Next, the analysis of association of imputed expression with BMI in the general and diabetic cohorts (Analysis 1 and 2) or GWG (Analysis 3 and 4) was performed, followed by variant association analysis (1 Mb around identified loci) with the mentioned phenotypes.

RESULTS

In Analysis 1 we have found 175 BMI associated genes and 19 variants (p < 10-4) which influenced GWG, with the strongest association for rs11465293 in CCL24 (p = 3.18E-06). Analysis 2, with diabetes included in the model, led to discovery of 1812 BMI associated loci and 207 variants (p < 10-4) influencing GWG, with the strongest association for rs9690213 in PODXL (p = 9.86E-07). In Analysis 3, among 648 GWG associated loci, 2091 variants were associated with BMI (FDR < 0.05). In Analysis 4, 7 variants in GWG associated loci influenced BMI in the ARIC cohort.

CONCLUSIONS

Here, we have shown that loci influencing BMI might have an impact on GWG and GWG associated loci might influence BMI, both in the general and T1DM cohorts. The results suggest that both phenotypes are related to insulin signaling, glucose homeostasis, mitochondrial metabolism, ubiquitinoylation and inflammatory responses.

Evaluation of Polygenic Risk Scores for Prediction of Prostate Cancer in Korean Men

Aims

The purpose of this study is to evaluate an aggregate influence of prostate cancer (PCa) susceptibility variants on the development of PCa in Korean men by using the polygenic risk score (PRS) approach.

Methods

An analysis of 1,001 cases of PCa and 2,641 controls was performed to: (i) identify potential PCa-related risk loci in Koreans and (ii) validate the cumulative association between these loci and PCa using the PRS. Subgroup analyses based on risk stratification were conducted to better characterize the potential correlation to key PCa-related clinical outcomes (e.g., Gleason score, prostate-specific antigen levels). The results were replicated using 514 cases of PCa and 548 controls from an independent cohort.

Results

Genome-wide association analysis from our discovery cohort revealed 11 candidate single-nucleotide polymorphisms (SNPs) associated with PCa showing statistical significance of p < 5.0 × 10^–5. Seven variants were located at 8q24.21 (rs1016343, rs16901979, and rs13252298 in PRNCR1; rs4242384, rs7837688, and rs1447295 in CASC8; and rs1512268 in NKX3). Two variants located within HNF1B (rs7501939 and rs4430796) had a significant negative association with PCa risk [odds ratio (OR) = 0.717 and 0.747, p = 6.42 × 10^–7 and 3.67 × 10^–6, respectively]. Of the six independent SNPs that remained after linkage disequilibrium (LD) pruning, the top four SNPs best predicted PCa risk with an area under the receiver operating characteristic curve (AUC) of 0.637 (95% CI: 0.582–0.692). Those with top 25% polygenic risk had a 4.2-fold increased risk of developing PCa compared with those with low risk.

Conclusion

Eleven PCa risk variants in Korean men were identified; PRSs of a subset of these variants could help predict PCa susceptibility.

Functional Genomics in Pancreatic β Cells: Recent Advances in Gene Deletion and Genome Editing Technologies for Diabetes Research

The inheritance of variants that lead to coding changes in, or the mis-expression of, genes critical to pancreatic beta cell function can lead to alterations in insulin secretion and increase the risk of both type 1 and type 2 diabetes. Recently developed clustered regularly interspaced short palindromic repeats (CRISPR/Cas9) gene editing tools provide a powerful means of understanding the impact of identified variants on cell function, growth, and survival and might ultimately provide a means, most likely after the transplantation of genetically "corrected" cells, of treating the disease. Here, we review some of the disease-associated genes and variants whose roles have been probed up to now. Next, we survey recent exciting developments in CRISPR/Cas9 technology and their possible exploitation for β cell functional genomics. Finally, we will provide a perspective as to how CRISPR/Cas9 technology may find clinical application in patients with diabetes.

A Genome-Wide Pharmacogenetic Study of Growth Hormone Responsiveness

CONTEXT

Individual patients vary in their response to growth hormone (GH). No large-scale genome-wide studies have looked for genetic predictors of GH responsiveness.

OBJECTIVE

To identify genetic variants associated with GH responsiveness.

DESIGN

Genome-wide association study (GWAS).

SETTING

Cohorts from multiple academic centers and a clinical trial.

PATIENTS

A total of 614 individuals from 5 short stature cohorts receiving GH: 297 with idiopathic short stature, 276 with isolated GH deficiency, and 65 born small for gestational age.

INTERVENTION

Association of more than 2 million variants was tested.

MAIN OUTCOME MEASURES

Primary analysis: individual single nucleotide polymorphism (SNP) association with first-year change in height standard deviation scores. Secondary analyses: SNP associations in clinical subgroups adjusted for clinical variables; association of polygenic score calculated from 697 genome-wide significant height SNPs with GH responsiveness.

RESULTS

No common variant associations reached genome-wide significance in the primary analysis. The strongest suggestive signals were found near the B4GALT4 and TBCE genes. After meta-analysis including replication data, signals at several loci reached or retained genome-wide significance in secondary analyses, including variants near ST3GAL6. There was no significant association with variants previously reported to be associated with GH response nor with a polygenic predicted height score.

CONCLUSIONS

We performed the largest GWAS of GH responsiveness to date. We identified 2 loci with a suggestive effect on GH responsiveness in our primary analysis and several genome-wide significant associations in secondary analyses that require further replication. Our results are consistent with a polygenic component to GH responsiveness, likely distinct from the genetic regulators of adult height.

Whole-exome sequencing reveals ANO8 as a genetic risk factor for intrahepatic cholestasis of pregnancy

Background

Intrahepatic cholestasis of pregnancy (ICP) is characterized by pruritus and cholestasis in late pregnancy and results in adverse pregnancy outcomes, including preterm delivery and birth weight, which are affected by the genetic and environmental background. However, until now, the genetic architecture of ICP has remained largely unclear.

Methods

Twenty-six clinical data points were recorded for 151 Chinese ICP patients. The data generated from whole-exome sequencing (WES) using the BGISEQ-500 platform were further analyzed by Burrows-Wheeler Aligner (BWA) software, Genome Analysis Toolkit (GATK), ANNOVAR tool, etc. R packages were used to conduct t-test, Fisher’s test and receiver operating characteristic (ROC) curve analyses.

Results

We identified eighteen possible pathogenic loci associated with ICP disease in known genes, covering ABCB4, ABCB11, ATP8B1 and TJP2. The loci Lys386Gln, Gly527Gln and Trp708Ter in ABCB4, Leu589Met, Gln605Pro and Gln1194Ter in ABCB11, and Arg189Ser in TJP2 were novel discoveries. In addition, WES analysis indicated that the gene ANO8 involved in the transport of bile salts is newly identified as associated with ICP. The functional network of the ANO8 gene confirmed this finding. ANO8 contained 8 rare missense mutations that were found in eight patients among the 151 cases and were absent from 1029 controls. Out of the eight SNPs, 3 were known, and the remaining five are newly identified. These variants have a low frequency, ranging from 0.000008 to 0.00001 in the ExAC, gnomAD – Genomes and TOPMED databases. Bioinformatics analysis showed that the sites and their corresponding amino acids were both highly conserved among vertebrates. Moreover, the influences of all the mutations on protein function were predicted to be damaging by the SIFT tool. Combining clinical data, it was found that the mutation group (93.36 µmol/L) had significantly (P = 0.038) higher total bile acid (TBA) levels than the wild-type group (40.81 µmol/L).

Conclusions

To the best of our knowledge, this is the first study to employ WES technology to detect genetic loci for ICP. Our results provide new insights into the genetic basis of ICP and will benefit the final identification of the underlying mutations.

A SNP-mediated lncRNA (LOC146880) and microRNA (miR-539-5p) interaction and its potential impact on the NSCLC risk

BACKGROUND

Many cancer-associated single nucleotide polymorphisms (SNPs) are located in the genomic regions of long non-coding RNAs (lncRNAs). Mechanisms of these SNPs in connection to cancer risk are not fully understood.

METHODS

Association of SNP (rs140618127) in lncRNA LOC146880 with non-small cell lung cancer (NSCLC) was evaluated in a case-control study of 2707 individuals. The mechanism of the SNP's biologic influence was explored with in vitro and in vivo experiments, including plasmid transfection, siRNA knockdown, flow cytometry assessment, and assays of cell proliferation, migration, invasion, and colony formation.

RESULTS

Association analysis showed that A allele of SNP rs140618127 was associated with low risk of NSCLC in the Chinese population. Lab experiments indicated that SNP rs140618127 contained a binding site for miR-539-5p and the binding between miR-539-5p and LOC146880 resulted in declined phosphorylation of an oncogene, ENO1. The reduced phosphorylation of ENO1 led to decreased phosphorylation of PI3K and Akt, which is further linked to the decline in cell proliferation and tumor progression.

CONCLUSION

The study demonstrates that SNP rs140618127 in lncRNA loc146880 provides an alternate binding site for microRNA miR-539-5p which affects the phosphorylation of ENO1 and activation of the PI3K and Akt pathway.

Genetic predisposition in type 2 diabetes: A promising approach toward a personalized management of diabetes

Diabetes mellitus, also known simply as diabetes, has been described as a chronic and complex endocrine metabolic disorder that is a leading cause of death across the globe. It is considered a key public health problem worldwide and one of four important non-communicable diseases prioritized for intervention through world health campaigns by various international foundations. Among its four categories, Type 2 diabetes (T2D) is the commonest form of diabetes accounting for over 90% of worldwide cases. Unlike monogenic inherited disorders that are passed on in a simple pattern, T2D is a multifactorial disease with a complex etiology, where a mixture of genetic and environmental factors are strong candidates for the development of the clinical condition and pathology. The genetic factors are believed to be key predisposing determinants in individual susceptibility to T2D. Therefore, identifying the predisposing genetic variants could be a crucial step in T2D management as it may ameliorate the clinical condition and preclude complications. Through an understanding the unique genetic and environmental factors that influence the development of this chronic disease individuals can benefit from personalized approaches to treatment. We searched the literature published in three electronic databases: PubMed, Scopus and ISI Web of Science for the current status of T2D and its associated genetic risk variants and discus promising approaches toward a personalized management of this chronic, non-communicable disorder.

Insights into pancreatic islet cell dysfunction from type 2 diabetes mellitus genetics

Type 2 diabetes mellitus (T2DM) is an increasingly prevalent multifactorial disease that has both genetic and environmental risk factors, resulting in impaired glucose homeostasis. Genome-wide association studies (GWAS) have identified over 400 genetic signals that are associated with altered risk of T2DM. Human physiology and epigenomic data support a central role for the pancreatic islet in the pathogenesis of T2DM. This Review focuses on the promises and challenges of moving from genetic associations to molecular mechanisms and highlights efforts to identify the causal variant and effector transcripts at T2DM GWAS susceptibility loci. In addition, we examine current human models that are used to study both β-cell development and function, including EndoC-β cell lines and human induced pluripotent stem cell-derived β-like cells. We use examples of four T2DM susceptibility loci (CDKAL1, MTNR1B, SLC30A8 and PAM) to emphasize how a holistic approach involving genetics, physiology, and cellular and developmental biology can disentangle disease mechanisms at T2DM GWAS signals.

PAM haploinsufficiency does not accelerate the development of diet- and human IAPP-induced diabetes in mice

AIMS/HYPOTHESIS

Peptide hormones are first synthesised as larger, inactive precursors that are converted to their active forms by endopeptidase cleavage and post-translational modifications, such as amidation. Recent, large-scale genome-wide studies have suggested that two coding variants of the amidating enzyme, peptidylglycine α-amidating monooxygenase (PAM), are associated with impaired insulin secretion and increased type 2 diabetes risk. We aimed to elucidate the role of PAM in modulating beta cell peptide amidation, beta cell function and the development of diabetes.

METHODS

PAM transcript and protein levels were analysed in mouse islets following induction of endoplasmic reticulum (ER) or cytokine stress, and PAM expression patterns were examined in human islets. To study whether haploinsufficiency of PAM accelerates the development of diabetes, Pam^+/- and Pam^+/+ mice were fed a low-fat diet (LFD) or high-fat diet (HFD) and glucose homeostasis was assessed. Since aggregates of the PAM substrate human islet amyloid polypeptide (hIAPP) lead to islet inflammation and beta cell failure, we also investigated whether PAM haploinsufficiency accelerated hIAPP-induced diabetes and islet amyloid formation in Pam^+/- and Pam^+/+ mice with beta cell expression of hIAPP.

RESULTS

Immunostaining revealed high expression of PAM in alpha, beta and delta cells in human pancreatic islets. Pam mRNA and PAM protein expression were reduced in mouse islets following administration of an HFD, and in isolated islets following induction of ER stress with thapsigargin, or cytokine stress with IL-1β, IFN-γ and TFN-α. Despite Pam^+/- only having 50% PAM expression and enzyme activity as compared with Pam^+/+ mice, glucose tolerance and body mass composition were comparable in the two models. After 24 weeks of HFD, both Pam^+/- and Pam^+/+ mice had insulin resistance and impaired glucose tolerance, but no differences in glucose tolerance, insulin sensitivity or plasma insulin levels were observed in PAM haploinsufficient mice. Islet amyloid formation and beta cell function were also similar in Pam^+/- and Pam^+/+ mice with beta cell expression of hIAPP.

CONCLUSIONS/INTERPRETATION

Haploinsufficiency of PAM in mice does not accelerate the development of diet-induced obesity or hIAPP transgene-induced diabetes.

Genome-wide meta-analysis associates GPSM1 with type 2 diabetes, a plausible gene involved in skeletal muscle function

Genome-wide association studies (GWASs) have identified many genetic variations associated with type 2 diabetes mellitus (T2DM) in Asians, but understanding the functional genetic variants that influence traits is often a complex process. In this study, fine mapping and other analytical strategies were performed to investigate the effects of G protein signaling modulator 1 (GPSM1) on insulin resistance in skeletal muscle. A total of 128 single-nucleotide polymorphisms (SNPs) within GPSM1 were analysed in 21,897 T2DM cases and 32,710 healthy controls from seven GWASs. The SNP rs28539249 in intron 9 of GPSM1 showed a nominally significant association with T2DM in Asians (OR = 1.07, 95% CI = 1.04-1.10, P < 10^-4). The GPSM1 mRNA was increased in skeletal muscle and correlated with T2DM traits across obese mice model. An eQTL for the cis-acting regulation of GPSM1 expression in human skeletal muscle was identified for rs28539249, and the increased GPSM1 expression related with T2DM traits within GEO datasets. Another independent Asian cohort showed that rs28539249 is associated with the skeletal muscle expression of CACFD1, GTF3C5, SARDH, and FAM163B genes, which are functionally enriched for endoplasmic reticulum stress (ERS) and unfolded protein response (UPR) pathways. Moreover, rs28539249 locus was predicted to disrupt regulatory regions in human skeletal muscle with enriched epigenetic marks and binding affinity for CTCF. Supershift EMSA assays followed luciferase assays demonstrated the CTCF specifically binding to rs28539249-C allele leading to decreased transcriptional activity. Thus, the post-GWAS annotation confirmed the Asian-specific association of genetic variant in GPSM1 with T2DM, suggesting a role for the variant in the regulation in skeletal muscle.

Population and evolutionary genetics of the PAH locus to uncover overdominance and adaptive mechanisms in phenylketonuria: Results from a multiethnic study

BACKGROUND

Phenylketonuria (PKU) is the most common inborn error of amino acid metabolism in Europe. The reasons underlying the high prevalence of heterozygous carriers are not clearly understood. We aimed to look for pathogenic PAH variant enrichment according to geographical areas and patients' ethnicity using a multiethnic nationwide cohort of patients with PKU in France. We subsequently appraised the population differentiation, balancing selection and the molecular evolutionary history of the PAH locus.

METHODS

The French nationwide PKU study included patients who have been referred at the national level to the University Hospital of Nancy, and for whom a molecular diagnosis of phenylketonuria was made by Sanger sequencing. We performed enrichment analyses by comparing alternative allele frequencies using Fisher's exact test with Bonferroni adjustment. We estimated the amount of genetic differentiation among populations using Wright's fixation index (Fst). To estimate the molecular evolutionary history of the PAH gene, we performed phylogenetic and evolutionary analyses using whole-genome and exome-sequencing data from healthy individuals and non-PKU patients, respectively. Finally, we used exome-wide association study to decipher potential genetic loci associated with population divergence on PAH.

FINDINGS

The study included 696 patients and revealed 132 pathogenic PAH variants. Three geographical areas showed significant enrichment for a pathogenic PAH variant: North of France (p.Arg243Leu), North-West of France (p.Leu348Val), and Mediterranean coast (p.Ala403Val). One PAH variant (p.Glu280Gln) was significantly enriched among North-Africans (OR = 23·23; 95% CI: 9·75-55·38). PAH variants exhibiting a strong genetic differentiation were significantly enriched in the 'Biopterin_H' domain (OR = 6·45; 95% CI: 1·99-20·84), suggesting a balancing selection pressure on the biopterin function of PAH. Phylogenetic and timetree analyses were consistent with population differentiation events on European-, African-, and Asian-ancestry populations. The five PAH variants most strongly associated with a high selection pressure were phylogenetically close and were located within the biopterin domain coding region of PAH or in its vicinity. Among the non-PAH loci potentially associated with population divergence, two reached exome-wide significance: SSPO (SCO-spondin) and DBH (dopamine beta-hydroxylase), involved in neuroprotection and metabolic adaptation, respectively.

INTERPRETATION

Our data provide evidence on the combination of evolutionary and adaptive events in populations with distinct ancestries, which may explain the overdominance of some genetic variants on PAH.

FUNDING

French National Institute of Health and Medical Research (INSERM) UMR_S 1256.

Gene Expression Analysis of the Pre-Diabetic Pancreas to Identify Pathogenic Mechanisms and Biomarkers of Type 1 Diabetes

Type 1 Diabetes (T1D) occurs as a result of the autoimmune destruction of pancreatic β-cells by self-reactive T cells. The etiology of this disease is complex and difficult to study due to a lack of disease-relevant tissues from pre-diabetic individuals. In this study, we performed gene expression analysis on human pancreas tissues obtained from the Network of Pancreatic Organ Donors with Diabetes (nPOD), and showed that 155 genes were differentially expressed by ≥2-fold in the pancreata of autoantibody-positive (AA+) at-risk individuals compared to healthy controls. Only 48 of these genes remained changed by ≥2-fold in the pancreata of established T1D patients. Pathway analysis of these genes showed a significant association with various immune pathways. We were able to validate the differential expression of eight disease-relevant genes by QPCR analysis: A significant upregulation of CADM2, and downregulation of TRPM5, CRH, PDK4, ANGPL4, CLEC4D, RSG16, and FCGR2B was confirmed in the pancreata of AA+ individuals versus controls. Studies have already implicated FCGR2B in the pathogenesis of disease in non-obese diabetic (NOD) mice. Here we showed that CADM2, TRPM5, PDK4, and ANGPL4 were similarly changed in the pancreata of pre-diabetic 12-week-old NOD mice compared to NOD.B10 controls, suggesting a possible role for these genes in the pathogenesis of both T1D and NOD disease. The loss of the leukocyte-specific gene, FCGR2B, in the pancreata of AA+ individuals, is particularly interesting, as it may serve as a potential whole blood biomarker of disease progression. To test this, we quantified FCGR2B expression in peripheral blood samples of T1D patients, and AA+ and AA- first-degree relatives of T1D patients enrolled in the TrialNet Pathway to Prevention study. We showed that FCGR2B was significantly reduced in the peripheral blood of AA+ individuals compared to AA- controls. Together, these findings demonstrate that gene expression analysis of pancreatic tissue and peripheral blood samples can be used to identify disease-relevant genes and pathways and potential biomarkers of disease progression in T1D.

Bioinformatics analysis of hepatic gene expression profiles in type 2 diabetes mellitus

Type 2 diabetes mellitus (T2DM) is characterized by hyperglycemia. The liver has a critical role in regulating glucose homeostasis. The present study aimed to analyze hepatic gene expression profiles and to identify the key genes and pathways involved in T2DM. Gene expression profiles of 10 patients with T2DM and 7 subjects with normal glucose tolerance were downloaded from the Gene Expression Omnibus database. Subsequently, differentially expressed genes (DEGs) were identified and functional enrichment analysis was performed. In addition, a protein-protein interaction network was built and hub genes were identified. In total, 1,320 DEGs were identified, including 698 up- and 622 downregulated genes, and these were mainly enriched in positive regulation of transcription from RNA polymerase II promoter, cell adhesion, inflammatory response, positive regulation of apoptotic process, signal transduction and the Tolllike receptor signaling pathway. A total of 8 hub genes (G-protein subunit gamma transducin 2, ubiquitinconjugating enzyme E2 D1, glutamate metabotropic receptor 1, G-protein signaling modulator 1, C-X-C motif chemokine ligand 9, neurotensin, purinergic receptor P2Y1 and ring finger protein 41) were screened from the network. The present study may contribute to the elucidation of the hepatic pathology of T2DM.

Exome-chip association analysis of intracranial aneurysms

OBJECTIVE

To investigate to what extent low-frequency genetic variants (with minor allele frequencies <5%) affect the risk of intracranial aneurysms (IAs).

METHODS

One thousand fifty-six patients with IA and 2,097 population-based controls from the Netherlands were genotyped with the Illumina HumanExome BeadChip. After quality control (QC) of samples and single nucleotide variants (SNVs), we conducted a single variant analysis using the Fisher exact test. We also performed the variable threshold (VT) test and the sequence kernel association test (SKAT) at different minor allele count (MAC) thresholds of >5 and >0 to test the hypothesis that multiple variants within the same gene are associated with IA risk. Significant results were tested in a replication cohort of 425 patients with IA and 311 controls, and results of the 2 cohorts were combined in a meta-analysis.

RESULTS

After QC, 995 patients with IA and 2,080 controls remained for further analysis. The single variant analysis comprising 46,534 SNVs did not identify significant loci at the genome-wide level. The gene-based tests showed a statistically significant association for fibulin 2 (FBLN2) (best p = 1 × 10^-6 for the VT test, MAC >5). Associations were not statistically significant in the independent but smaller replication cohort (p > 0.57) but became slightly stronger in a meta-analysis of the 2 cohorts (best p = 4.8 × 10^-7 for the SKAT, MAC ≥1).

CONCLUSION

Gene-based tests indicated an association for FBLN2, a gene encoding an extracellular matrix protein implicated in vascular wall remodeling, but independent validation in larger cohorts is warranted. We did not identify any significant associations for single low-frequency genetic variants.

Effects of variants of 50 genes on diabetes risk among the Chinese population born in the early 1960s

BACKGROUND

Genome-wide association studies have identified loci that significantly increase diabetes risk. This study explored the genetic susceptibility in relation to diabetes risk in adulthood among a Chinese population born in the early 1960s.

METHODS

In all, 2129 subjects (833 males, 1296 females) were selected from the cross-sectional 2010 to 2012 China National Nutrition and Health Survey. Fifty diabetes-related single nucleotide polymorphisms (SNPs) were detected. Two diabetes genetic risk scores (GRSs) based on the 50 diabetes-predisposing variants were developed to examine the association of these SNPs with diabetes risk.

RESULTS

Associations were found between diabetes risk and SNPs in the MTNR1B (rs10830963), KLHDC5 (rs10842994), GRK5 (rs10886471), cyclindependentkinase 5 regulatory subunit associated protein 1 (rs10946398), adaptorrelated protein complex 3 subunit sigma 2 (rs2028299), diacylglycerol kinase beta/transmembrane protein 195 (rs2191349), SREBF chaperone (rs4858889), ankyrin1 (rs516946), RAS guanyl releasing protein 1 (rs7403531), and zinc finger AN1-type containing 3 (rs9470794) genes. As a continuous variable, with a 1-point increase in the GRS or weighted (w) GRS, fasting plasma glucose (FPG) increased 0.045 and 0.044 mM, respectively (P < 0.001 for both), after adjusting for confounders. Both GRS and wGRS showed an association with diabetes, with a multivariable-adjusted odds ratio (95% confidence interval) of 1.09 (1.00-1.19) and 1.12 (1.03-1.22), respectively, among all subjects. No significant associations were found between the GRS or wGRS and impaired fasting glucose or impaired glucose tolerance.

CONCLUSIONS

The data suggest the association of 10 SNPs and the GRS or wGRS with diabetes risk. Genetic susceptibility to diabetes may synergistically affect the risk of diabetes in adulthood.