This article has Open Peer Review reports available.
Comprehensive genomic analysis identifies pathogenic variants in maturity-onset diabetes of the young (MODY) patients in South India
- Viswanathan Mohan†1Email author,
- Venkatesan Radha†1,
- Thong T. Nguyen†2,
- Eric W. Stawiski†2, 3,
- Kanika Bajaj Pahuja2,
- Leonard D. Goldstein2, 3,
- Jennifer Tom3,
- Ranjit Mohan Anjana1,
- Monica Kong-Beltran2,
- Tushar Bhangale3, 4,
- Suresh Jahnavi1,
- Radhakrishnan Chandni5,
- Vijay Gayathri1,
- Paul George6,
- Na Zhang2,
- Sakthivel Murugan6,
- Sameer Phalke6,
- Subhra Chaudhuri2,
- Ravi Gupta6,
- Jingli Zhang2,
- Sam Santhosh6,
- Jeremy Stinson2,
- Zora Modrusan2,
- V. L. Ramprasad6,
- Somasekar Seshagiri2Email authorView ORCID ID profile and
- Andrew S. Peterson2Email author
© The Author(s). 2018
Received: 25 August 2017
Accepted: 19 January 2018
Published: 13 February 2018
Maturity-onset diabetes of the young (MODY) is an early-onset, autosomal dominant form of non-insulin dependent diabetes. Genetic diagnosis of MODY can transform patient management. Earlier data on the genetic predisposition to MODY have come primarily from familial studies in populations of European origin.
In this study, we carried out a comprehensive genomic analysis of 289 individuals from India that included 152 clinically diagnosed MODY cases to identify variants in known MODY genes. Further, we have analyzed exome data to identify putative MODY relevant variants in genes previously not implicated in MODY. Functional validation of MODY relevant variants was also performed.
We found MODY 3 (HNF1A; 7.2%) to be most frequently mutated followed by MODY 12 (ABCC8; 3.3%). They together account for ~ 11% of the cases. In addition to known MODY genes, we report the identification of variants in RFX6, WFS1, AKT2, NKX6–1 that may contribute to development of MODY. Functional assessment of the NKX6–1 variants showed that they are functionally impaired.
Our findings showed HNF1A and ABCC8 to be the most frequently mutated MODY genes in south India. Further we provide evidence for additional MODY relevant genes, such as NKX6–1, and these require further validation.
Maturity-onset diabetes of the young (MODY) refers to a group of monogenic, non-insulin dependent forms of diabetes that affect non-obese young adults with a prevalence of ~ 100 per million in European populations [1, 2]. MODY has a high inheritance rate and represents one end of a continuum of monogenic forms of diabetes that includes neonatal diabetes [1–3]. Seven of the known MODY genes, the largest functional class, encode transcription factors which are involved in development, maturation and maintenance of β-cells while the remaining are primarily involved in the process of insulin release in response to glucose [1–3]. Although molecular genetic studies have identified heterozygous causal variants in 14 MODY genes, it is believed that additional MODY genes remain to be identified [1, 4].
MODY is often misdiagnosed as either type 1 or type 2 diabetes . Thus, an accurate molecular genetic diagnosis of monogenic diabetes forms like MODY can have a significant impact on patient care. For example, individuals with mutations in GCK do not ordinarily need pharmacological treatment whereas those with mutations in KCNJ11 or ABCC8 are often best managed with high dose sulfonylureas rather than with insulin, and those with mutations in HNF1A or HNF4A benefit from treatment with low dose sulfonylureas [6, 7]. A recent study on diabetes in youth reported that MODY was misdiagnosed in 36% of the cases as type 1 (T1DM) and in 51% of the cases as type 2 (T2DM), highlighting the need for the inclusion of molecular genetic diagnosis in the management of young diabetic patients .
Given the opportunity for accurate diagnosis and informed treatment, genetic testing for MODY is an ideal test case for effective implementation of precision health care [1, 7]. For those individuals whose age of onset is less than 30 years, molecular genetic testing for variants in known MODY genes has the potential to provide significant patient benefit . For those with later onset or for those who are young and asymptomatic, it is not clear whether genotyping would be beneficial since most studies of MODY have focused on affected individuals and closely related family members, leading to inflated estimates of effect sizes. A recent study of 7 MODY genes in ~ 4000 individuals of European and African-American ancestry identified variants that are predicted to be damaging at a general population frequency of 0.5–1.5% . Damaging variants were however substantially enriched in cohorts of MODY patients or in young, lean subjects with diabetes thus reiterating the potential value of genetic diagnosis in guiding treatment for those diagnosed with diabetes by conventional criteria. A significant caveat to the interpretations outlined above is that the overwhelming majority of sequencing studies that seek to understand disease-causing variants and estimate disease burden have focused on individuals of Northern European descent . These studies have found that ~ 70–80% of causal variants are in either GCK or HNF1A. Small-scale studies of South Asian MODY patients have indicated that substantially fewer MODY cases can be explained by variants in HNF1A [11, 12], HNF4A , GCK  and HNF1B . Given the distinct clinical profile of diabetes in the Indian subcontinent [16, 17] it is not surprising that a different pattern of MODY inheritance may exist and indeed South Asians may harbor variants in yet to be identified MODY loci. The large population, the significant burden of Mendelian disorders, the very high incidence of diabetes and the earlier age of onset in the Indian subcontinent  make it particularly important to understand the pattern of MODY gene variation in India.
In this study, we have performed a comprehensive characterization of the spectrum of MODY gene variation in a cohort of patients from South India. Our genomic analysis and functional testing indicate NKX6–1 to be a potential MODY gene.
Samples, DNA and RNA preps
Diabetic subjects were recruited from a large diabetes center in Chennai (formerly Madras) city in southern India. All patients underwent a structured assessment including detailed family history and met MODY clinical criteria of Fajans and Tattersal : age at diagnosis of diabetes 30 years or less, control of hyperglycemia for a minimum period of 2 years without insulin, negative for auto antibodies, absence of ketonuria at any time and evidence of autosomal dominant inheritance including a three-generation family history of diabetes.
Clinical characteristics of the MODY patients and those of the normal glucose tolerant controls are given in Additional file 1: Table S1a. Written informed consent was obtained from all the study participants. The study was approved by the Institutional Ethics Committee of the Madras Diabetes Research Foundation. The reported investigations have been carried out in accordance with the principles of the Declaration of Helsinki.
Consented and de-identified patient blood samples were used for extraction of DNA. EDTA anti-coagulated venous blood samples were collected from all study subjects, and the genomic DNA was isolated from whole blood by proteinase K digestion followed by phenol-chloroform extraction. Subsequently genomic DNA was precipitated in ethanol. The quality and quantity were assessed spectrophotometrically.
A total of 289 samples were used in this study. This consisted of 152 clinically diagnosed MODY cases and 137 normal glucose tolerant (NGT) subjects used for variant filtering.
Whole genome, exome and targeted sequencing
Whole genome libraries were constructed using TruSeqNano DNA Library Preparation Kit (Illumina, CA) and sequenced on Illumina HiSeq2500 (Illumina, CA) to generate 2 × 75 bp read.
The whole exome analysis was performed using Agilent SureSelect (Santa Clara, CA) Human All Exome kit v5 (50 Mb). Exome capture libraries were sequenced on HiSeq 2500 (Illumina, CA). Targeted exome sequencing of 258 samples (includes 121 MODY samples and 137 normal glucose tolerant subjects) was performed using custom probes corresponding to 1965 genes that included established MODY genes and those that were implicated in pancreatic cell biology and/or diabetes (Additional file 1: Table S3).
Sequence data processing
All sequencing reads were evaluated for quality using the Bioconductor ShortRead package . An all-against-all sample comparison was done on germline variants to confirm the patient exome-WGS data pairing and establish sample relatedness prior to additional data analysis.
Whole genome and whole exome sequencing data were combined and then processed using the Genome Analysis Toolkit (GATK) (version v3.5–0-g36282e4) best practices recommendations [21, 22]. Briefly, reads were mapped to the human reference genome GRCh37 using BWA-MEM (version 0.7.10; http://arxiv.org/abs/1303.3997). Duplicate alignments were marked and removed using Picard tool (version 1.126) (http://broadinstitute.github.io/picard/) followed by indel realignment and base quality score recalibration. The GATK Haplotype Caller algorithm was used to generate gVCFs for all samples. Joint variant calling was performed separately for discovery cohort samples and validation cohort samples using GATK Genotype GVCFs. Variant quality score recalibration was carried out to estimate the confidence of variants called in the discovery cohort. For the validation cohort, standard hard filters were applied using a set of criteria and parameters as recommended by GATK . Variant annotation was carried out using SnpEff (version 4.1) .
Our variant filtering strategy is outlined in Fig. 2a. It involved sequential filtering steps followed by manual review. We filtered for variants that meet the following three criteria: variants that (1) were rare (MAF < = 1% in 1000 genomes  and the NHLBI Exome Sequencing Project (ESP) ; (2) were protein-altering or potentially protein-altering; (3) and were in the curated list of 35 genes that included known MODY genes and others implicated in early onset diabetes including neonatal diabetes (Additional file 1: Table S3). In the manual review process, the set of variants obtained following filtering were manually checked for read evidence using Integrative Genomics Viewer (IGV) . We retained ultra-rare variants, which have an allele-frequency < 0.01% in ExAC and were present in at the most only one sample in the ExAC South Asian cohort. Finally, we checked if any of our filtered variants were reported in previous MODY studies as disease relevant and annotated them.
Plasmids, stable cell lines and expression studies
Mutant NKX6-1 constructs were generated using a site-directed mutagenesis kit (Stratagene, USA) in TOPO-TA vector backbone. The mutant and wildtype sequences were sub-cloned from TOPO-TA vector into pLVX-TetOne-Puro vector (Clonetech, USA; Cat No. 631847) or pRK5 vector (Genentech, CA). Lenti-X 293T packaging cells (Clontech, USA; Cat No. 632180) were transfected with clones carrying variants in the pLVX-TetOne-Puro vector using Lenti-X Packaging Single Shots (VSV-G) (Clontech, USA; Cat No. 631275). At 48h post transfection, lentivirus containing supernatants were used to infect β-TC-6 cells according to the standard protocol in the Lenti-X Tet-One Inducible Expression Systems User Manual (Clontech, USA) for 24 h. 2 μg/mL puromycin (Gibco, USA; Cat No. A11138–03) containing β-TC-6 complete growth media with Tet-free FBS (Clontech, USA; Cat No. 631101) was used to select infected cells for 7 days. Post selection, β-TC-6 stable cell lines were maintained in complete growth media supplemented with 0.5 μg/mL puromycin. All the stable cell lines used in the subsequent studies were passaged for at least 3 times and tested with Lenti-X p24 Rapid Titer Kit (Clontech, Cat No. 632200) for the clearance of lentivirus. Stable cell lines were expanded and cultured in complete growth media with 15% Tet-free FBS for 24 h. We seeded 50,000 cells in a 12 well plate and induced gene expression using 100 ng/ml doxycycline (dox). Cells were harvested immediately following induction (0 h) and at 48 h post induction. Cells were lysed and total RNA from cell lines was extracted using RNeasy Plus Mini kit (Qiagen, CA). For each stable cell line expressing a particular gene we prepared three biological replicates for each condition (wt and variants) at 0 h (no Dox) and 48 h (with Dox) time point.
RNA-Seq and data analysis
We used 0.5μg of total RNA to generate RNA-seq libraries using TruSeq RNA Sample Preparation kit (Illumina, CA). The libraries were multiplexed and sequenced on HiSeq 2500 to obtain on average ~ 30 million single end 50 bp reads per sample.
RNA-seq data from mouse cell lines with human expression constructs were mapped to the mouse genome (GRCm38/mm10) and respective wildtype human construct. Mapping was performed with GSNAP (2013–10-10) . Only uniquely mapped reads were considered for downstream analysis. Gene models were based on RefSeq transcripts and NCBI gene annotation. Differential expression was performed with the R/Bioconductor package DESeq2 . For clustering, count data were transformed by variance stabilization and genes were centered to have mean zero. Clustering was performed using 1-Pearson correlation as distance metric and average linkage. Heat maps were generated using the R/Bioconductor package NMF.
Western blot analysis
Expression of NKX6–1 was tested using western blot as previously described . Briefly, β-TC-6 stable cells (50,000) expressing NKX6–1 were harvested at 48 h post induction with 100 ng/ml dox. Lysates were prepared for western blotting as previously described . We used anti-NKX6–1 (Novus Biologicals; NBP-149672) and anti-Hsp90 (Santa Cruz Biotechnology; sc-7947) antibodies for protein detection. Appropriate HRP-conjugated secondary antibodies (Pierce Biotechnology, IL) were used along with the primary antibodies in the western blot assay.
Identification of MODY relevant variants
We performed whole or targeted exome sequencing of our 152 MODY cases to assess the presence of MODY relevant variation. All samples were confirmed to be unrelated (Additional file 2: Figure S2). We initially performed whole exome sequencing on 31 MODY probands (Additional file 1: Table S1a). We also performed low-pass whole genome sequencing (mean coverage ~ 6×) of the 31 probands. We combined reads from exome and low-pass WGS and achieved > 46× coverage for the 31 samples (Additional file 1: Table S2a). Principal component analysis  using the genotypes called from low-pass WGS, merged with those from 1000 Genomes  project phase 3 samples, showed that the 31 samples overlap with Indian reference samples, as expected (Fig. 1d) confirming the ethnicity of the samples used in the study. We assessed the remaining 121 MODY patient samples for presence of pathogenic variants using targeted exome sequencing (see Methods). We obtained an average coverage of 82× for the 121 samples (Additional file 1: Table S2b).
Targeted exome sequencing and analysis of the 121 MODY samples resulted 2,792,168 variants. We applied the same filtering strategy as done for the whole exome samples above and finally arrived at 29 variants (Fig. 2a; Additional file 1: Table S4).
Combined together, we identified 39 candidate variants in 39 of the 152 MODY samples (Fig. 2b; Additional file 1: Table S4). Of these, variants R272H and V103M found in HNF1A were observed in two and three patient samples, respectively (Fig. 2b). Of the 14 known MODY genes , 9 were found mutated in our samples (Fig. 2b). HNF1A (MODY 3) was most frequently mutated at 7.2% (11/152), followed by ABCC8 at 3.3% (MODY 12; 5/152). Other MODY genes were found to be mutated in one sample and this includes BLK (MODY 11), CEL (MODY 8), GCK (MODY 2), HNF1B (MODY 5), KCNJ11 (MODY 13), KLF11 (MODY 7), and PDX1 (MODY 4). In addition to the known MODY genes, we found candidate variants in 10 MODY-relevant genes. This includes RFX6 (n = 3), WFS1 (n = 3), NKX6–1 (n = 2), AKT2 (n = 2) (Fig. 2b; Additional file 1: Table S4). All variants identified, except for WFS1 P320R, were heterozygous, consistent with the dominant inheritance pattern of MODY.
NKX6–1 variants found in MODY patients are functionally impaired
For each of the two NKX6–1 variants, and the functionally impaired control mutant, we found significant differences in gene expression between cells expressing the variant as compared to cells expressing wild type NKX6–1 (adjusted p < 0.05, fold-change > 1.5; Additional file 1: Table S5). Notably, clustering of samples based on genes identified as differentially expressed in at least one of the three comparisons showed that the transcriptional activity of the two MODY relevant variants was similar to the functionally impaired control mutant (Fig. 4c). This suggests that S317L and P329L NKX6–1 variants were functionally impaired and that disruption of the carboxy-terminal domain of NKX6–1 is likely responsible for MODY in these individuals (Fig. 4c). While the function of the carboxy terminus of NKX6–1 is poorly understood, deletions and substitutions have been shown to disrupt the ability of the nearby homeobox domain to bind DNA with appropriate specificity .
We found 39 genes that were differentially regulated in common between the two NKX6–1 MODY mutants (Additional file 1: Table S5), suggesting that some of these may be related to the development of diabetes in these patients. Further, inspection of the genes that are dysregulated by the two variants and comparison to NKX6–1 binding sites in the β-TC-6 cell line (MKB, JZ and ASP, unpublished) revealed a number of genes which may be related to the development of diabetes in these patients (Additional file 1: Table S5). For example, BHLAHA15 (MIST1), with a NKX6–1 binding site upstream of the transcription start site, is significantly down-regulated by both variants. BHLHA15 knockout mice have decreased expression of Glut2 in pancreatic β-cells and age-related impairment in glucose clearance , suggesting that its down-regulation could contribute to the development of diabetes in these patients. Interestingly, the most significantly up-regulated gene for all three variants was MAP3K15 (ASK3). The function of MAP3K15 has not been studied in any detail but it is highly homologous to Apoptosis-Signal-regulated Kinase 1 (ASK1), a gene that stimulates apoptosis in response to ER stress and other noxious stimuli , suggesting that an increased rate of apoptotic β cell loss may also be involved in the development of diabetes in patients carrying the NKX6–1 variants.
Overall, we analyzed 152 clinically diagnosed unrelated MODY patients from southern India for the presence of variants in known MODY genes and additional loci of potential MODY-relevance. We found relevant variants in previously described MODY genes in ~ 15% (23/152) of the cases analyzed (Fig. 2b). As has been previously reported, the locus diversity of MODY in India is unlike that in European populations where up to 80% of patients carry variants in HNF1A and GCK. We found HNF1A variants in ~ 7% of the patients and GCK alterations in < 1% of the patient samples, a substantially lower rate than in Europeans and consistent with previous reports of Indian MODY patients [11, 12, 14, 37–39]. We found ~ 3% of the samples carried an ABCC8 variant (MODY 12). The ABCC8 variants have important treatment implications, as it has been suggested that these patients may be best managed with the use of sulfonylureas [40, 41], similar to the situation with cases of neonatal diabetes caused by mutations in KCNJ11 or ABCC8.
Importantly, we found 6 variants, including three PDX1 variants P33T , E224K [12, 43], P242L , HNF4A V169I [12, 45], BLK A71T [12, 46] and NEUROD1 H241Q [12, 47], previously reported to be linked to MODY, to occur at a similar frequencies in the MODY cases and the general population. This suggests that these variants most likely have little to no effect on risk for developing MODY. The findings underscore the difficulty in interpreting relevance of variants where population specific common variant data are not available and highlights the need for such data to enable the application of genomics in the clinic.
In addition to known MODY genes we found rare variants of functional consequence in genes implicated in other forms of diabetes or beta-cell biology. This includes variants in NKX6-1 (S317L, P329L). A role for NKX6–1 in MODY is consistent with its role in β-cell function, maintenance and/or development and the phenotypes described in knockout mice studies [48, 49]. The NKX6-1 S317 and P329 are highly conserved residues across orthologs from multiple organisms and substitution at these positions is likely to impair its normal function. Consistent with this, both S317L and P329 were found to be functionally impaired in their transcriptional activity (Fig. 4). These findings while suggest a plausible role NKX6–1 in MODY further studies would be needed to confirm its role as a MODY gene.
In addition to NKX6–1 we found variants in WFS1, RFX6 and 7 other genes associated with monogenic forms of diabetes (AKT2, EIF2AK3, GLIS3, HADH, MNX1, NKX2–2, and PTF1A) and they may be relevant to development of MODY. Additional genetic or functional studies will be required to evaluate and establish the relevance of these genes in MODY.
In conclusion, our study reports that HNF1A and ABCC8 are among the most frequently mutated MODY genes in south India. Further we provide evidence for additional MODY relevant genes, such as NKX6–1, and these require further validation. Our findings have important implications for MODY gene testing.
The authors would like to acknowledge Genentech DNA Sequencing, Oligo and Bioinformatics groups for their help with the project and members of the Departments of Molecular Biology and Human Genetics for helpful discussions and review of the manuscript. The authors would like to thank Ying Littman for technical support. The authors wish to thank the Indian Council of Medical Research (ICMR), New Delhi, India, for financial support to VR for Sanger sequencing under the grant schemes ‘Advanced Centre for Genomics of Type 2 Diabetes and Genetic Analysis of MODY and Neonatal diabetes in India’.
Availability of data and materials
Sequencing and genotype data for patients specifically consenting to this have been deposited at the European Genome-phenome Archive, which is hosted by the European Bioinformatics Institute (EBI), under accession EGAS00001001699 (URL: https://www.ebi.ac.uk/ega/studies/EGAS00001001699).
ASP, SS, TTN, VR and VM designed the study. VR, VG and CR helped in collection of samples.VM and RMA carried out clinical evaluations of MODY patients and controls. VR and SJ performed DNA extraction of study samples. PG and JT managed and analysed clinical data. TB carried out genetic analyses of data. SM, SP, SSanthosh, SJ, ZM, JAS, and VLR designed the plan for generation of sequence data and participated in its analysis. TTN, ES, RG,SJ, VR, ZM and ASP managed and analysed sequence data. KBP, MKB, NZ, SC and JZ carried out functional studies. LG analysed RNA-seq data. SS and ASP co-wrote the manuscript with the assistance and contributions of all of the other authors. All authors read and approved the final manuscript.
Ethics approval and consent to participate
Written informed consent was obtained from all the study participants. Consent from the Parents or legal guardian was obtained wherever samples from minors were included in the study. The study was approved by the Institutional Ethics Committee of the Madras Diabetes Research Foundation. The reported investigations have been carried out in accordance with the principles of the Declaration of Helsinki.
Consent for publication
All Genentech Inc. authors hold shares in Roche.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Anik A, Catli G, Abaci A, Bober E. Maturity-onset diabetes of the young (MODY): an update. J Pediatr Endocrinol Metab. 2015;28(3–4):251–63.PubMedGoogle Scholar
- Schwitzgebel VM. Many faces of monogenic diabetes. J Diabetes Investig. 2014;5(2):121–33.View ArticlePubMedPubMed CentralGoogle Scholar
- Ashcroft FM, Rorsman P. Diabetes mellitus and the beta cell: the last ten years. Cell. 2012;148(6):1160–71.View ArticlePubMedGoogle Scholar
- Prudente S, Jungtrakoon P, Marucci A, Ludovico O, Buranasupkajorn P, Mazza T, Hastings T, Milano T, Morini E, Mercuri L, et al. Loss-of-function mutations in APPL1 in familial diabetes mellitus. Am J Hum Genet. 2015;97(1):177–85.View ArticlePubMedPubMed CentralGoogle Scholar
- Murphy R, Ellard S, Hattersley AT. Clinical implications of a molecular genetic classification of monogenic beta-cell diabetes. Nat Clin Pract Endocrinol Metab. 2008;4(4):200–13.View ArticlePubMedGoogle Scholar
- Naylor R, Philipson LH. Who should have genetic testing for maturity-onset diabetes of the young? Clin Endocrinol. 2011;75(4):422–6.View ArticleGoogle Scholar
- Thanabalasingham G, Pal A, Selwood MP, Dudley C, Fisher K, Bingley PJ, Ellard S, Farmer AJ, McCarthy MI, Owen KR. Systematic assessment of etiology in adults with a clinical diagnosis of young-onset type 2 diabetes is a successful strategy for identifying maturity-onset diabetes of the young. Diabetes Care. 2012;35(6):1206–12.View ArticlePubMedPubMed CentralGoogle Scholar
- Pihoker C, Gilliam LK, Ellard S, Dabelea D, Davis C, Dolan LM, Greenbaum CJ, Imperatore G, Lawrence JM, Marcovina SM, et al. Prevalence, characteristics and clinical diagnosis of maturity onset diabetes of the young due to mutations in HNF1A, HNF4A, and glucokinase: results from the SEARCH for diabetes in youth. J Clin Endocrinol Metab. 2013;98(10):4055–62.View ArticlePubMedPubMed CentralGoogle Scholar
- Flannick J, Beer NL, Bick AG, Agarwala V, Molnes J, Gupta N, Burtt NP, Florez JC, Meigs JB, Taylor H, et al. Assessing the phenotypic effects in the general population of rare variants in genes for a dominant Mendelian form of diabetes. Nat Genet. 2013;45(11):1380–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Johansson S, Irgens H, Chudasama KK, Molnes J, Aerts J, Roque FS, Jonassen I, Levy S, Lima K, Knappskog PM, et al. Exome sequencing and genetic testing for MODY. PLoS One. 2012;7(5):e38050.View ArticlePubMedPubMed CentralGoogle Scholar
- Radha V, Ek J, Anuradha S, Hansen T, Pedersen O, Mohan V. Identification of novel variants in the hepatocyte nuclear factor-1alpha gene in south Indian patients with maturity onset diabetes of young. J Clin Endocrinol Metab. 2009;94(6):1959–65.View ArticlePubMedGoogle Scholar
- Chapla A, Mruthyunjaya MD, Asha HS, Varghese D, Varshney M, Vasan SK, Venkatesan P, Nair V, Mathai S, Paul TV, et al. Maturity onset diabetes of the young in India - a distinctive mutation pattern identified through targeted next-generation sequencing. Clin Endocrinol. 2015;82(4):533–42.View ArticleGoogle Scholar
- Anuradha S, Radha V, Mohan V. Association of novel variants in the hepatocyte nuclear factor 4A gene with maturity onset diabetes of the young and early onset type 2 diabetes. Clin Genet. 2011;80(6):541–9.View ArticlePubMedGoogle Scholar
- Kanthimathi S, Jahnavi S, Balamurugan K, Ranjani H, Sonya J, Goswami S, Chowdhury S, Mohan V, Radha V. Glucokinase gene mutations (MODY 2) in Asian Indians. Diabetes Technol Ther. 2014;16(3):180–5.View ArticlePubMedGoogle Scholar
- Kanthimathi S, Balamurugan K, Mohan V, Shanthirani CS, Gayathri V, Radha V. Identification and molecular characterization of HNF1B gene mutations in Indian diabetic patients with renal abnormalities. Ann Hum Genet. 2015;79(1):10–9.View ArticlePubMedGoogle Scholar
- Reich D, Thangaraj K, Patterson N, Price AL, Singh L. Reconstructing Indian population history. Nature. 2009;461(7263):489–94.View ArticlePubMedPubMed CentralGoogle Scholar
- Majumder PP, Basu A. A genomic view of the peopling and population structure of India. Cold Spring Harb Perspect Biol. 2015;7(4):a008540.View ArticlePubMed CentralGoogle Scholar
- Anjana RM, Pradeepa R, Deepa M, Datta M, Sudha V, Unnikrishnan R, Bhansali A, Joshi SR, Joshi PP, Yajnik CS, et al. Prevalence of diabetes and prediabetes (impaired fasting glucose and/or impaired glucose tolerance) in urban and rural India: phase I results of the Indian Council of Medical Research-INdia DIABetes (ICMR-INDIAB) study. Diabetologia. 2011;54(12):3022–7.View ArticlePubMedGoogle Scholar
- Tattersall RB, Fajans SS. A difference between the inheritance of classical juvenile-onset and maturity-onset type diabetes of young people. Diabetes. 1975;24(1):44–53.View ArticlePubMedGoogle Scholar
- Morgan M, Anders S, Lawrence M, Aboyoun P, Pages H, Gentleman R. ShortRead: a bioconductor package for input, quality assessment and exploration of high-throughput sequence data. Bioinformatics. 2009;25(19):2607–8.View ArticlePubMedPubMed CentralGoogle Scholar
- DePristo MA, Banks E, Poplin R, Garimella KV, Maguire JR, Hartl C, Philippakis AA, del Angel G, Rivas MA, Hanna M, et al. A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat Genet. 2011;43(5):491–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Van der Auwera GA, Carneiro MO, Hartl C, Poplin R, Del Angel G, Levy-Moonshine A, Jordan T, Shakir K, Roazen D, Thibault J, et al. From FastQ data to high confidence variant calls: the genome analysis toolkit best practices pipeline. Curr Protoc Bioinformatics. 2013;43:11 10 11–33.Google Scholar
- Cingolani P, Platts A, Wang le L, Coon M, Nguyen T, Wang L, Land SJ, Lu X, Ruden DM. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila Melanogaster strain w1118; iso-2; iso-3. Fly (Austin). 2012;6(2):80–92.View ArticleGoogle Scholar
- 1000 Genomes Project Consortium, Auton A, Brooks LD, Durbin RM, Garrison EP, Kang HM, Korbel JO, Marchini JL, McCarthy S, McVean GA, et al. A global reference for human genetic variation. Nature. 2015;526(7571):68–74.View ArticleGoogle Scholar
- Fu W, O'Connor TD, Jun G, Kang HM, Abecasis G, Leal SM, Gabriel S, Rieder MJ, Altshuler D, Shendure J, et al. Analysis of 6,515 exomes reveals the recent origin of most human protein-coding variants. Nature. 2013;493(7431):216–20.View ArticlePubMedGoogle Scholar
- Robinson JT, Thorvaldsdottir H, Winckler W, Guttman M, Lander ES, Getz G, Mesirov JP. Integrative genomics viewer. Nat Biotechnol. 2011;29(1):24–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Wu TD, Nacu S. Fast and SNP-tolerant detection of complex variants and splicing in short reads. Bioinformatics. 2010;26(7):873–81.View ArticlePubMedPubMed CentralGoogle Scholar
- Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.View ArticlePubMedPubMed CentralGoogle Scholar
- Jaiswal BS, Kljavin NM, Stawiski EW, Chan E, Parikh C, Durinck S, Chaudhuri S, Pujara K, Guillory J, Edgar KA, et al. Oncogenic ERBB3 mutations in human cancers. Cancer Cell. 2013;23(5):603–17.View ArticlePubMedGoogle Scholar
- Price AL, Patterson NJ, Plenge RM, Weinblatt ME, Shadick NA, Reich D. Principal components analysis corrects for stratification in genome-wide association studies. Nat Genet. 2006;38(8):904–9.View ArticlePubMedGoogle Scholar
- Flanagan SE, De Franco E, Lango Allen H, Zerah M, Abdul-Rasoul MM, Edge JA, Stewart H, Alamiri E, Hussain K, Wallis S, et al. Analysis of transcription factors key for mouse pancreatic development establishes NKX2-2 and MNX1 mutations as causes of neonatal diabetes in man. Cell Metab. 2014;19(1):146–54.View ArticlePubMedPubMed CentralGoogle Scholar
- Sansbury FH, Kirel B, Caswell R, Lango Allen H, Flanagan SE, Hattersley AT, Ellard S, Shaw-Smith CJ. Biallelic RFX6 mutations can cause childhood as well as neonatal onset diabetes mellitus. Eur J Hum Genet. 2015;23(12):1744–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Taylor DG, Babu D, Mirmira RG. The C-terminal domain of the beta cell homeodomain factor Nkx6.1 enhances sequence-selective DNA binding at the insulin promoter. Biochemistry. 2005;44(33):11269–78.View ArticlePubMedGoogle Scholar
- Mirmira RG, Watada H, German MS. Beta-cell differentiation factor Nkx6.1 contains distinct DNA binding interference and transcriptional repression domains. J Biol Chem. 2000;275(19):14743–51.View ArticlePubMedGoogle Scholar
- Fazio EN, Everest M, Colman R, Wang R, Pin CL. Altered Glut-2 accumulation and beta-cell function in mice lacking the exocrine-specific transcription factor, Mist1. J Endocrinol. 2005;187(3):407–18.View ArticlePubMedGoogle Scholar
- Shiizaki S, Naguro I, Ichijo H. Activation mechanisms of ASK1 in response to various stresses and its significance in intracellular signaling. Advances in biological regulation. 2013;53(1):135–44.View ArticlePubMedGoogle Scholar
- Anuradha S, Radha V, Deepa R, Hansen T, Carstensen B, Pedersen O, Mohan V. A prevalent amino acid polymorphism at codon 98 (Ala98Val) of the hepatocyte nuclear factor-1alpha is associated with maturity-onset diabetes of the young and younger age at onset of type 2 diabetes in Asian Indians. Diabetes Care. 2005;28(10):2430–5.View ArticlePubMedGoogle Scholar
- Chavali S, Mahajan A, Tabassum R, Dwivedi OP, Chauhan G, Ghosh S, Tandon N, Bharadwaj D. Association of variants in genes involved in pancreatic beta-cell development and function with type 2 diabetes in north Indians. J Hum Genet. 2011;56(10):695–700.View ArticlePubMedGoogle Scholar
- Sahu RP, Aggarwal A, Zaidi G, Shah A, Modi K, Kongara S, Aggarwal S, Talwar S, Chu S, Bhatia V, et al. Etiology of early-onset type 2 diabetes in Indians: islet autoimmunity and mutations in hepatocyte nuclear factor 1alpha and mitochondrial gene. J Clin Endocrinol Metab. 2007;92(7):2462–7.View ArticlePubMedGoogle Scholar
- Artuso R, Provenzano A, Mazzinghi B, Giunti L, Palazzo V, Andreucci E, Blasetti A, Chiuri RM, Gianiorio FE, Mandich P, et al. Therapeutic implications of novel mutations of the RFX6 gene associated with early-onset diabetes. Pharmacogenomics J. 2015;15(1):49–54.View ArticlePubMedGoogle Scholar
- Bowman P, Flanagan SE, Edghill EL, Damhuis A, Shepherd MH, Paisey R, Hattersley AT, Ellard S. Heterozygous ABCC8 mutations are a cause of MODY. Diabetologia. 2012;55(1):123–7.View ArticlePubMedGoogle Scholar
- Gragnoli C, Stanojevic V, Gorini A, Von Preussenthal GM, Thomas MK, Habener JF. IPF-1/MODY4 gene missense mutation in an Italian family with type 2 and gestational diabetes. Metabolism. 2005;54(8):983–8.View ArticlePubMedGoogle Scholar
- Cockburn BN, Bermano G, Boodram LL, Teelucksingh S, Tsuchiya T, Mahabir D, Allan AB, Stein R, Docherty K, Bell GI. Insulin promoter factor-1 mutations and diabetes in Trinidad: identification of a novel diabetes-associated mutation (E224K) in an indo-Trinidadian family. J Clin Endocrinol Metab. 2004;89(2):971–8.View ArticlePubMedGoogle Scholar
- Ang SF, Lim SC, Tan C, Fong JC, Kon WY, Lian JX, Subramanium T, Sum CF. A preliminary study to evaluate the strategy of combining clinical criteria and next generation sequencing (NGS) for the identification of monogenic diabetes among multi-ethnic Asians. Diabetes Res Clin Pract. 2016;119:13–22.View ArticlePubMedGoogle Scholar
- Gragnoli C, Menzinger Von Preussenthal G, Habener JF. Triple genetic variation in the HNF-4alpha gene is associated with early-onset type 2 diabetes mellitus in a philippino family. Metabolism. 2004;53(8):959–63.View ArticlePubMedGoogle Scholar
- Borowiec M, Liew CW, Thompson R, Boonyasrisawat W, Hu J, Mlynarski WM, El Khattabi I, Kim SH, Marselli L, Rich SS, et al. Mutations at the BLK locus linked to maturity onset diabetes of the young and beta-cell dysfunction. Proc Natl Acad Sci U S A. 2009;106(34):14460–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Gonsorcikova L, Pruhova S, Cinek O, Ek J, Pelikanova T, Jorgensen T, Eiberg H, Pedersen O, Hansen T, Lebl J. Autosomal inheritance of diabetes in two families characterized by obesity and a novel H241Q mutation in NEUROD1. Pediatr Diabetes. 2008;9(4 Pt 2):367–72.View ArticlePubMedGoogle Scholar
- Schaffer AE, Taylor BL, Benthuysen JR, Liu J, Thorel F, Yuan W, Jiao Y, Kaestner KH, Herrera PL, Magnuson MA, et al. Nkx6.1 controls a gene regulatory network required for establishing and maintaining pancreatic Beta cell identity. PLoS Genet. 2013;9(1):e1003274.View ArticlePubMedPubMed CentralGoogle Scholar
- Sander M, Sussel L, Conners J, Scheel D, Kalamaras J, Dela Cruz F, Schwitzgebel V, Hayes-Jordan A, German M: Homeobox gene Nkx6.1 lies downstream of Nkx2.2 in the major pathway of beta-cell formation in the pancreas. Development 2000, 127(24):5533-5540.Google Scholar