A low-frequency GLIS3 variant associated with resistance to Japanese type 1 diabetes

Partners in Crime: Beta-Cells and Autoimmune Responses Complicit in Type 1 Diabetes Pathogenesis

Type 1 diabetes (T1D) is an autoimmune disease characterized by autoreactive T cell-mediated destruction of insulin-producing pancreatic beta-cells. Loss of beta-cells leads to insulin insufficiency and hyperglycemia, with patients eventually requiring lifelong insulin therapy to maintain normal glycemic control. Since T1D has been historically defined as a disease of immune system dysregulation, there has been little focus on the state and response of beta-cells and how they may also contribute to their own demise. Major hurdles to identifying a cure for T1D include a limited understanding of disease etiology and how functional and transcriptional beta-cell heterogeneity may be involved in disease progression. Recent studies indicate that the beta-cell response is not simply a passive aspect of T1D pathogenesis, but rather an interplay between the beta-cell and the immune system actively contributing to disease. Here, we comprehensively review the current literature describing beta-cell vulnerability, heterogeneity, and contributions to pathophysiology of T1D, how these responses are influenced by autoimmunity, and describe pathways that can potentially be exploited to delay T1D.

β-Cell failure in diabetes: Common susceptibility and mechanisms shared between type 1 and type 2 diabetes

Diabetes mellitus is etiologically classified into type 1, type 2 and other types of diabetes. Despite distinct etiologies and pathogenesis of these subtypes, many studies have suggested the presence of shared susceptibilities and underlying mechanisms in β-cell failure among different types of diabetes. Understanding these susceptibilities and mechanisms can help in the development of therapeutic strategies regardless of the diabetes subtype. In this review, we discuss recent evidence indicating the shared genetic susceptibilities and common molecular mechanisms between type 1, type 2 and other types of diabetes, and highlight the future prospects as well.

Molecular and clinical genetics of the transcription factor GLIS3 in Chinese congenital hypothyroidism

The transcription factor GLIS3 is an important factor in hormone biosynthesis and thyroid development, and mutations in GLIS3 are relatively rare. Deletions of more than one of the 11 exons of GLIS3 occur in most patients with various extrathyroidal abnormalities and congenital hypothyroidism (CH), and only 18 missense variants of GLIS3 related to thyroid disease have been reported. The aim of this study was to report the family history and molecular basis of patients with CH who carry GLIS3 variants. Three hundred and fifty-three non-consanguineous infants with CH were recruited and subjected to targeted exome sequencing of CH-related genes. The transcriptional activity and cellular localization of the variants in GLIS3 were investigated in vitro. We identified 20 heterozygous GLIS3 exonic missense variants, including eight novel sites, in 19 patients with CH. One patient carried compound heterozygous GLIS3 variants (p.His34Arg and p.Pro835Leu). None of the variants affected the nuclear localization. However, three variants (p.His34Arg, p.Pro835Leu, and p.Ser893Phe) located in the N-terminal and C-terminal regions of the GLIS3 protein downregulated the transcriptional activation of several genes required for thyroid hormone (TH) biosynthesis. This study of patients with CH extends the current knowledge surrounding the spectrum of GLIS3 variants and the mechanisms by which they cause TH biosynthesis defects.

Association study of the functional variants of the GLIS3 gene with risk of knee osteoarthritis

INTRODUCTION/OBJECTIVES

Osteoarthritis (OA) ranks the most common joint disorder and the leading cause of disability. Growing evidence has revealed that OA has a strong genetic background, except for aging and obesity. The aim of this study is to determine the associations between potential functional variants of the GLIS3 and GLIS3-AS1 gene and risk of knee OA among a Chinese population.

METHODS

In this case-control study with 810 knee OA cases and 900 healthy controls, seven selected functional SNPs of the GLIS3 and GLIS3-AS1 gene were evaluated.

RESULTS

We found minor alleles of rs10116772 (OR: 0.80, 95% CI: 0.69-0.92, P = 0.002), rs7045410 (OR: 0.74, 95% CI: 0.61-0.92, P = 0.005), and rs7032713 (OR: 0.76, 95% CI: 0.63-0.93, P = 0.006) were significantly associated with decreased risk of knee OA. Results of the dominant and recessive model, stratified analyses using Kellgren-Lawrence (KL) grading presented that the significant associations were not materially changed. Haplotype analysis indicated that haplotype CGT (OR: 0.66, 95% CI: 0.46-0.96, P = 0.031) and ATT (OR: 0.76, 95% CI: 0.6-0.95, P = 0.017) were significantly associated with decreased risk of knee OA. Further, they were also significantly associated with lower expression level of GLIS3, as well as higher expression level of GLIS3-AS1 in the articular cartilage specimens. Genotype-tissue expression (GTEX) data also validated that minor alleles of rs7045410 and rs7032713 were significantly associated with higher expression level of GLIS3-AS1 in thyroid and pituitary tissues (P < 0.001).

CONCLUSIONS

These findings revealed the essential role of genetic variants of the GLIS3 and GLIS3-AS1 gene in the occurrence of knee OA together. Key Point • Functional variants of the GLIS3 and GLIS3-AS1 gene were significantly associated with decreased risk of knee OA.

Sequencing reveals protective and pathogenic effects on development of diabetes of rare GLIS3 variants

Background

Based on the association of common GLIS3 variants with various forms of diabetes and the biological role of GLIS3 in beta-cells, we sequenced GLIS3 in non-diabetic and diabetic Danes to investigate the effect of rare missense variants on glucose metabolism.

Methods

We sequenced 53 patients with maturity-onset diabetes of the young (MODY), 5,726 non-diabetic participants, 2,930 patients with newly diagnosed type 2 diabetes and 206 patients with glutamic acid decarboxylase antibody (GADA) -positive diabetes.

Results

In total we identified 86 rare (minor allele frequency < 0.1%) missense variants. None was considered causal for the presence of MODY. Among patients with type 2 diabetes, we observed a higher prevalence of rare GLIS3 missense variants (2.5%) compared to non-diabetic individuals (1.8%) (odds ratio of 1.37 (interquartile range:1.01–1.88, p = 0.04)). A significantly increased HbA1c was found among patients with type 2 diabetes and with GADA-positive diabetes carrying rare GLIS3 variants compared to non-carriers of rare GLIS3 variants with diabetes (p = 0.02 and p = 0.004, respectively). One variant (p.I28V) was found to have a minor allele frequency of only 0.03% among patients with type 2 diabetes compared to 0.2% among non-diabetic individuals suggesting a protective function (odds ratio of 0.20 (interquartile range: 0.005–1.4, p = 0.1)), an effect which was supported by publically available data. This variant was also associated with a lower level of fasting plasma glucose among non-diabetic individuals (p = 0.046).

Conclusion

Rare missense variants in GLIS3 associates nominally with increased level of HbA1c and increased risk of developing type 2 diabetes. In contrast, the rare p.I28V variant associate with reduced level of fasting plasma glucose and may be protective against type 2 diabetes.

HisCoM-GGI: Hierarchical structural component analysis of gene-gene interactions

Although genome-wide association studies (GWAS) have successfully identified thousands of single nucleotide polymorphisms (SNPs) associated with common diseases, these observations are limited for fully explaining "missing heritability". Determining gene-gene interactions (GGI) are one possible avenue for addressing the missing heritability problem. While many statistical approaches have been proposed to detect GGI, most of these focus primarily on SNP-to-SNP interactions. While there are many advantages of gene-based GGI analyses, such as reducing the burden of multiple-testing correction, and increasing power by aggregating multiple causal signals across SNPs in specific genes, only a few methods are available. In this study, we proposed a new statistical approach for gene-based GGI analysis, "Hierarchical structural CoMponent analysis of Gene-Gene Interactions" (HisCoM-GGI). HisCoM-GGI is based on generalized structured component analysis, and can consider hierarchical structural relationships between genes and SNPs. For a pair of genes, HisCoM-GGI first effectively summarizes all possible pairwise SNP-SNP interactions into a latent variable, from which it then performs GGI analysis. HisCoM-GGI can evaluate both gene-level and SNP-level interactions. Through simulation studies, HisCoM-GGI demonstrated higher statistical power than existing gene-based GGI methods, in analyzing a GWAS of a Korean population for identifying GGI associated with body mass index. Resultantly, HisCoM-GGI successfully identified 14 potential GGI, two of which, (NCOR2 <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mo>×</mml:mo></mml:math> SPOCK1) and (LINGO2 <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mo>×</mml:mo></mml:math> ZNF385D) were successfully replicated in independent datasets. We conclude that HisCoM-GGI method may be a valuable tool for genome to identify GGI in missing heritability, allowing us to better understand the biological genetic mechanisms of complex traits. We conclude that HisCoM-GGI method may be a valuable tool for genome to identify GGI in missing heritability, allowing us to better understand biological genetic mechanisms of complex traits. An implementation of HisCoM-GGI can be downloaded from the website ( http://statgen.snu.ac.kr/software/hiscom-ggi ).

GLIS1-3 transcription factors: critical roles in the regulation of multiple physiological processes and diseases

Krüppel-like zinc finger proteins form one of the largest families of transcription factors. They function as key regulators of embryonic development and a wide range of other physiological processes, and are implicated in a variety of pathologies. GLI-similar 1-3 (GLIS1-3) constitute a subfamily of Krüppel-like zinc finger proteins that act either as activators or repressors of gene transcription. GLIS3 plays a critical role in the regulation of multiple biological processes and is a key regulator of pancreatic β cell generation and maturation, insulin gene expression, thyroid hormone biosynthesis, spermatogenesis, and the maintenance of normal kidney functions. Loss of GLIS3 function in humans and mice leads to the development of several pathologies, including neonatal diabetes and congenital hypothyroidism, polycystic kidney disease, and infertility. Single nucleotide polymorphisms in GLIS3 genes have been associated with increased risk of several diseases, including type 1 and type 2 diabetes, glaucoma, and neurological disorders. GLIS2 plays a critical role in the kidney and GLIS2 dysfunction leads to nephronophthisis, an end-stage, cystic renal disease. In addition, GLIS1-3 have regulatory functions in several stem/progenitor cell populations. GLIS1 and GLIS3 greatly enhance reprogramming efficiency of somatic cells into induced embryonic stem cells, while GLIS2 inhibits reprogramming. Recent studies have obtained substantial mechanistic insights into several physiological processes regulated by GLIS2 and GLIS3, while a little is still known about the physiological functions of GLIS1. The localization of some GLIS proteins to the primary cilium suggests that their activity may be regulated by a downstream primary cilium-associated signaling pathway. Insights into the upstream GLIS signaling pathway may provide opportunities for the development of new therapeutic strategies for diabetes, hypothyroidism, and other diseases.

Genome-Wide Study of Subcutaneous and Visceral Adipose Tissue Reveals Novel Sex-Specific Adiposity Loci in Mexican Americans

OBJECTIVE

This study aimed to explore the genetic mechanisms of regional fat deposition, which is a strong risk factor for metabolic diseases beyond total adiposity.

METHODS

A genome-wide association study of 7,757,139 single-nucleotide polymorphisms (SNPs) in 983 Mexican Americans (n_male  = 403; n_female  = 580) from the Insulin Resistance Atherosclerosis Family Study was performed. Association analyses were performed with and without sex stratification for subcutaneous adipose tissue, visceral adipose tissue (VAT), and visceral-subcutaneous ratio (VSR) obtained from computed tomography.

RESULTS

The strongest signal identified was SNP rs2185405 (minor allele frequencies [MAF] = 40%; P_VAT = 1.98 × 10^-8 ) with VAT. It is an intronic variant of the GLIS family zinc finger 3 gene (GLIS3). In addition, SNP rs12657394 (MAF = 19%) was associated with VAT in males (P_male = 2.39×10^-8 ; P_female = 2.5 × 10^-3 ). It is located intronically in the serum response factor binding protein 1 gene (SRFBP1). On average, male carriers of the variant had 24.6 cm² increased VAT compared with noncarriers. Subsequently, genome-wide SNP-sex interaction analysis was performed. SNP rs10913233 (MAF = 14%; P_int = 3.07 × 10^-8 ) in PAPPA2 and rs10923724 (MAF = 38%; P_int = 2.89 × 10^-8 ) upstream of TBX15 were strongly associated with the interaction effect for VSR.

CONCLUSIONS

Six loci were identified with genome-wide significant associations with fat deposition and interactive effects. These results provided genetic evidence for a differential basis of fat deposition between genders.

GLIS1-3: emerging roles in reprogramming, stem and progenitor cell differentiation and maintenance

Recent studies have provided evidence for a regulatory role of GLI-similar (GLIS) transcription factors in reprogramming, maintenance and differentiation of several stem and progenitor cell populations. GLIS1, in conjunction with several other reprogramming factors, was shown to markedly increase the efficiency of generating induced pluripotent stem cells (iPSC) from somatic cells. GLIS2 has been reported to contribute to the maintenance of the pluripotent state in hPSCs. In addition, GLIS2 has a function in regulating self-renewal of hematopoietic progenitors and megakaryocytic differentiation. GLIS3 plays a critical role during the development of several tissues. GLIS3 is able to promote reprogramming of human fibroblasts into retinal pigmented epithelial (RPE) cells. Moreover, GLIS3 is essential for spermatogonial stem cell renewal and spermatogonial progenitor cell differentiation. During pancreas development, GLIS3 protein is first detectable in bipotent pancreatic progenitors and pro-endocrine progenitors and plays a critical role in the generation of pancreatic beta cells. Here, we review the current status of the roles of GLIS proteins in the maintenance and differentiation of these different stem and progenitor cells.

Molecular genetics of the transcription factor GLIS3 identifies its dual function in beta cells and neurons

The GLIS family zinc finger 3 isoform (GLIS3) is a risk gene for Type 1 and Type 2 diabetes, glaucoma and Alzheimer's disease endophenotype. We identified GLIS3 binding sites in insulin secreting cells (INS1) (FDR q<0.05; enrichment range 1.40-9.11 fold) sharing the motif wrGTTCCCArTAGs, which were enriched in genes involved in neuronal function and autophagy and in risk genes for metabolic and neuro-behavioural diseases. We confirmed experimentally Glis3-mediated regulation of the expression of genes involved in autophagy and neuron function in INS1 and neuronal PC12 cells. Naturally-occurring coding polymorphisms in Glis3 in the Goto-Kakizaki rat model of type 2 diabetes were associated with increased insulin production in vitro and in vivo, suggestive alteration of autophagy in PC12 and INS1 and abnormal neurogenesis in hippocampus neurons. Our results support biological pleiotropy of GLIS3 in pathologies affecting β-cells and neurons and underline the existence of trans‑nosology pathways in diabetes and its co-morbidities.

GLIS3 rs7020673 and rs10758593 polymorphisms interact in the susceptibility for type 1 diabetes mellitus

AIMS

The transcription factor Gli-similar 3 (GLIS3) plays a key role in the development and maintenance of pancreatic beta cells as well as in the regulation of Insulin gene expression in adults. Accordingly, genome-wide association studies identified GLIS3 as a susceptibility locus for type 1 diabetes mellitus (T1DM) and glucose metabolism traits. Therefore, the aim of this study was to replicate the association of the rs10758593 and rs7020673 single nucleotide polymorphisms (SNPs) in the GLIS3 gene with T1DM in a Brazilian population.

METHODS

Frequencies of the rs7020673 (G/C) and rs10758593 (A/G) SNPs were analyzed in 503 T1DM patients (cases) and in 442 non-diabetic subjects (controls). Haplotypes constructed from the combination of these SNPs were inferred using a Bayesian statistical method.

RESULTS

Genotype and allele frequencies of rs7020673 and rs10758593 SNPs did not differ significantly between case and control groups. However, the frequency of ≥3 minor alleles of the analyzed SNPs in haplotypes was higher in T1DM patients compared to non-diabetic subjects (6.2 vs. 1.6%; P = 0.001). The presence of ≥3 minor alleles remained independently associated with risk of T1DM after adjustment for T1DM high-risk HLA DR/DQ haplotypes, age and ethnicity (OR = 3.684 95% CI 1.220-11.124). Moreover, levels of glycated hemoglobin seem to be higher in T1DM patients with rs10758593 A/A genotype than patients carrying the G allele of this SNP (P = 0.038), although this association was not kept after Bonferroni correction.

CONCLUSIONS

Our results indicate that individually the rs7020673 and rs10758593 SNPs are not significantly associated with T1DM but seem to interact in the predisposition for this disease.

Assessing variation across 8 established East Asian loci for type 2 diabetes mellitus in American Indians: Suggestive evidence for new sex-specific diabetes signals in GLIS3 and ZFAND3

BACKGROUND

Eight new loci for type 2 diabetes mellitus (T2DM) were identified in an East Asian genome-wide association study meta-analysis. We assess tag SNPs across these loci for associations with T2DM in American Indians.

METHODS

A total of 435 SNPs that tag (R²  ≥ .85) common variation across the 8 loci were analyzed for association with T2DM (n = 7710), early onset T2DM (n = 1060), body mass index (n = 6839), insulin sensitivity (n = 555), and insulin secretion (n = 298).

RESULTS

Tag SNPs within FITM2-R3HDML-HNF4A, GLIS3, KCNK16, and ZFAND3 associated with T2DM after accounting for locus-wide multiple testing. The T2DM association in FITM2-R3HDML-HNF4A (rs3212183; P = .0002; OR = 1.19 [1.09-1.30]) was independent from the East Asian lead SNP (rs6017317), which did not associate with T2DM in American Indians. The top signals in GLIS3 (rs7875253; P = .0004; OR = 1.23 [1.10-1.38]) and KCNK16 (rs1544050; P = .002; OR = 1.16 [1.06-1.27]) were attenuated after adjustment for the East Asian lead SNPs (rs7041847 in GLIS3; rs1535500 in KCNK16), both of which also associated with T2DM in American Indians (P = .02; OR = 1.11 [1.01-1.21]; P = .007; OR = 1.19 [1.05-1.36] respectively). The top SNP in ZFAND3 (rs9470794; P = .002; OR = 1.43 [1.14-1.80]) was the identical East Asian lead SNP. Additional SNPs in GLIS3 (rs180867004) and ZFAND3 (rs4714120 and rs9470701) had significant genotype × sex interactions (P ≤ .008). The GLIS3 SNP (rs180867004) associated with T2DM only in men (P = .00006, OR = 1.94 [1.40-2.68]). The ZFAND3 SNPs (rs4714120 and rs9470701) associated with T2DM only in women (P = .0002, OR = 1.35 [1.16-1.59]; P = .0003, OR = 1.37 [1.16-1.63] respectively).

CONCLUSIONS

Replication of lead T2DM SNPs in GLIS3, KCNK16, and ZFAND3 was observed in American Indians. Sex-specific T2DM signals in GLIS3 and ZFAND3, which are distinct from the East Asian GWAS signals, were also identified.

The role of GLIS3 in thyroid disease as part of a multisystem disorder

Congenital hypothyroidism is the most common hereditary endocrine disorder. In a small number of cases, mutations have been identified that are associated with maldevelopment and maldescent of the thyroid. Some of these mutations present as syndromes with a multisystem phenotype such as NKX2-1, PAX8, and FOXE. The association of permanent neonatal diabetes and congenital hypothyroidism was first reported in 2003 and subsequently led to the identification GLIS3 as the mutation responsible for this presentation. GLIS3 is a member of the GLI-similar zinc finger protein family encoding for a nuclear protein with five zinc finger domains and maps to chromosome 9p24. Given the role of GLIS3 in transcriptional activation and repression during embryogenesis, in humans, GLIS3 mutations present with multisystem involvement that also includes renal cystic dysplasia, progressive liver fibrosis and osteopenia. Thyroid findings in GLIS3 patients include thyroid aplasia, diminished colloid with interstitial fibrosis at post-mortem, and apparently normal gross thyroid anatomy on ultrasonography but with temporary TSH resistance on treatment. To date no biological mechanism has explained this variable presentation.

Emerging roles of GLIS3 in neonatal diabetes, type 1 and type 2 diabetes

GLI-similar 3 (GLIS3), a member of the Krüppel-like zinc finger protein subfamily, is predominantly expressed in the pancreas, thyroid and kidney. Glis3 mRNA can be initially detected in mouse pancreas at embryonic day 11.5 and is largely restricted to β cells, pancreatic polypeptide-expressing cells, as well as ductal cells at later stage of pancreas development. Mutations in GLIS3 cause a neonatal diabetes syndrome, characterized by neonatal diabetes, congenital hypothyroidism and polycystic kidney. Importantly, genome-wide association studies showed that variations of GLIS3 are strongly associated with both type 1 diabetes (T1D) and type 2 diabetes (T2D) in multiple populations. GLIS3 cooperates with pancreatic and duodenal homeobox 1 (PDX1), v-maf musculoaponeurotic fibrosarcoma oncogene family, protein A (MAFA), as well as neurogenic differentiation 1 (NEUROD1) and potently controls insulin gene transcription. GLIS3 also plays a role in β cell survival and likely in insulin secretion. Any perturbation of these functions may underlie all three forms of diabetes. GLIS3, synergistically with hepatocyte nuclear factor 6 (HNF6) and forkhead box A2 (FOXA2), controls fetal islet differentiation via transactivating neurogenin 3 (NGN3) and impairment of this function leads to neonatal diabetes. In addition, GLIS3 is also required for the compensatory β cell proliferation and mass expansion in response to insulin resistance, which if disrupted may predispose to T2D. The increasing understanding of the mechanisms of GLIS3 in β cell development, survival and function maintenance will provide new insights into disease pathogenesis and potential therapeutic target identification to combat diabetes.

A common variant association study reveals novel susceptibility loci for low HDL-cholesterol levels in ethnic Arabs

The genetic susceptibility to acquiring low high density lipoprotein-cholesterol (LHDLC) levels is not completely elucidated yet. In this study, we performed a common variant association study for harboring this trait in ethnic Arabs. We employed the Affymetrix high-density Axiom Genome-Wide ASI Array (Asian population) providing a coverage of 598,000 single nucleotide variations (SNPs) to genotype 5495 individuals in a two-phase study involving discovery and validation sets of experiments. The rs1800775 [1.31 (1.22-1.42); p = 3.41E-12] in the CETP gene and rs359027 [1.26 (1.16-1.36); p = 2.55E-08] in the LMCD1 gene were significantly associated with LHDLC levels. Furthermore, rs3104435 [1.26 (1.15-1.38); p = 1.19E-06] at the MATN1 locus, rs9835344 [1.16 (1.08-1.26); p = 8.75E-06] in the CNTN6 gene, rs1559997 [1.3 (1.14-1.47); p = 9.48E-06] in the SDS gene and rs1670273 [1.2 (1.1-1.31); p = 4.81E-06] in the DMN/SYNM gene exhibited suggestive association with the disorder. Seven other variants including rs1147169 in the PLCL1 gene, rs10248618 in the DNAH11, rs476155 in the GLIS3, rs7024300 in the ABCA1, intergenic rs10836699, rs11603691 in P2RX3 and rs750134 in CORO1C gene exhibited borderline protective properties. Validation and joint meta-analysis resulted in rs1800775, rs3104435 and rs359027 retaining their predisposing properties, while rs10836699 and rs11603691 showed protective properties. Our data show several predisposing variants across the genome for LHDLC levels in ethnic Arabs.

A genome-wide association study for diabetic retinopathy in a Japanese population: potential association with a long intergenic non-coding RNA

Elucidation of the genetic susceptibility factors for diabetic retinopathy (DR) is important to gain insight into the pathogenesis of DR, and may help to define genetic risk factors for this condition. In the present study, we conducted a three-stage genome-wide association study (GWAS) to identify DR susceptibility loci in Japanese patients, which comprised a total of 837 type 2 diabetes patients with DR (cases) and 1,149 without DR (controls). From the stage 1 genome-wide scan of 446 subjects (205 cases and 241 controls) on 614,216 SNPs, 249 SNPs were selected for the stage 2 replication in 623 subjects (335 cases and 288 controls). Eight SNPs were further followed up in a stage 3 study of 297 cases and 620 controls. The top signal from the present association analysis was rs9362054 in an intron of RP1-90L14.1 showing borderline genome-wide significance (P_met = 1.4×10⁻⁷, meta-analysis of stage 1 and stage 2, allele model). RP1-90L14.1 is a long intergenic non-coding RNA (lincRNA) adjacent to KIAA1009/QN1/CEP162 gene; CEP162 plays a critical role in ciliary transition zone formation before ciliogenesis. The present study raises the possibility that the dysregulation of ciliary-associated genes plays a role in susceptibility to DR.