Next Article in Journal
A Perspective on the Link between Mitochondria-Associated Membranes (MAMs) and Lipid Droplets Metabolism in Neurodegenerative Diseases
Previous Article in Journal
Biological Activity of Cyclic Peptide Extracted from Sphaeranthus amaranthoides Using De Novo Sequencing Strategy by Mass Spectrometry for Cancer
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Whole-Exome Sequencing in Family Trios Reveals De Novo Mutations Associated with Type 1 Diabetes Mellitus

1
Center of Biotechnology, Khalifa University of Science and Technology, Abu Dhabi 127788, United Arab Emirates
2
Department of Biomedical Engineering, Khalifa University of Science and Technology, Abu Dhabi 127788, United Arab Emirates
3
College of Medicine and Health Sciences, Khalifa University, Abu Dhabi 127788, United Arab Emirates
4
Department of Endocrinology, Mafraq Hospital, Abu Dhabi 127788, United Arab Emirates
*
Author to whom correspondence should be addressed.
Biology 2023, 12(3), 413; https://doi.org/10.3390/biology12030413
Submission received: 27 December 2022 / Revised: 16 February 2023 / Accepted: 23 February 2023 / Published: 7 March 2023

Abstract

:

Simple Summary

Type 1 diabetes mellitus (T1DM) is a chronic autoimmune condition in which the immune system destroys insulin-making cells in the pancreas. Many advances have been made in the past decade to understand the pathophysiology of T1DM. With an estimated heritability risk of 50%, the strong genetic component plays an important role in the discovery of novel disease pathways and identification of new targets for therapeutic purposes. In this study, we aim to identify new (de novo) genetic markers for T1DM patients by sequencing the genes of the affected individual and their parents (trio family). This is a powerful approach to identify causal mutations for inherited diseases, such as T1DM, to improve our understanding of the condition. With 13 trio families, we identified 32 new (de novo) genetic mutations. Of these, 12 variants that were linked to T1DM, and the remaining 20 variants were linked to endocrine, metabolic, or autoimmune diseases. The findings of this study have allowed us to identify the genetic markers associated with the development of T1DM, to be able to improve diagnosis through therapeutic advancements.

Abstract

Type 1 diabetes mellitus (T1DM) is a chronic autoimmune disease characterized by insulin deficiency and loss of pancreatic islet β-cells. The objective of this study is to identify de novo mutations in 13 trios from singleton families that contribute to the genetic basis of T1DM through the application of whole-exome sequencing (WES). Of the 13 families sampled for this project, 12 had de novo variants, with Family 7 having the highest number (nine) of variants linked to T1DM/autoimmune pathways, whilst Family 4 did not have any variants past the filtering steps. There were 10 variants of 7 genes reportedly associated with T1DM (MST1; TDG; TYRO3; IFIHI; GLIS3; VEGFA; TYK2). There were 20 variants of 13 genes that were linked to endocrine, metabolic, or autoimmune diseases. Our findings demonstrate that trio-based WES is a powerful approach for identifying new candidate genes for the pathogenesis of T1D. Genotyping and functional annotation of the discovered de novo variants in a large cohort is recommended to ascertain their association with disease pathogenesis.

1. Introduction

Type 1 diabetes mellitus (T1DM) is a chronic autoimmune disease primarily characterized by insulin deficiency and loss of pancreatic islet β-cells [1,2]. T1DM accounts for 5–10% of total cases of diabetes worldwide, and is one of the most common endocrine and metabolic conditions occurring in childhood. The peak age of presentation of childhood-onset T1DM has a bimodal distribution, with the first peak between four to six years, and the second in early puberty [3]. The pathogenesis of T1DM has been suggested to be a continuum development, starting from the detection of autoantibodies before symptom onset to the progression of β-cell destruction, dysglycaemia, and hyperglycemia [4].
T1DM is a heterogeneous multifactorial condition, and the elucidation of its complex etiology is highly dependent on the interaction between numerous environmental factors precipitated by genetic susceptibility. Environmental influences, such as lifestyle, viral infections (especially respiratory such as rotavirus, cytomegalovirus, and mumps or coxsackie B and other enterovirus infections), and gestational events have been proposed as candidate etiological factors. In individuals with genetic predisposition, exposure to these environmental triggers leads to an autoimmune response, mediated by autoreactive CD4 and CD8 T-cells, resulting in the destruction of insulin-producing beta cells [5].
Family studies have revealed that T1DM is heritable, with an estimated 50% risk. The concordance rate among monozygotic twins is reported to be 30% within 10 years of diagnosis of the first twin, reaching as high as 65% by 60 years of age [6]. A child of an affected mother has a lifetime rise between 1 and 4%, whereas a child of an affected father has a 3–8% risk, and where both parents are affected, the risk reaches up to 30% [7]. Linkage and genetic association studies have identified >50 loci that contribute to susceptibility or resistance to developing T1DM, including its association with pancreatic beta cell autoimmunity (insulin gene (INS)), inflammatory-associated factors (interleukin-2 signaling pathway, B-cell and T-cell development, and cytokine signaling), human leucocyte antigen (HLA) genes, and shared genetic architecture with autoimmune diseases (APS1 and STAT3 poly-autoimmunity) [8,9,10].
Despite the rapid advancements in genetic methods, our ability to understand the pathogenesis of T1DM to improve therapeutic or diagnostic potential is lacking. The increasing availability of whole-exome sequencing (WES) platforms has substantially improved the identification and investigation of the role of inherited and de novo variants in T1DM-associated genes. In this study, we implemented WES to identify variants in a cohort of 50 samples comprising 13 case–parent trios from singleton families with T1DM in the UAE.

2. Materials and Methods

2.1. Study Participants and Recruitment

Thirteen families of UAE nationality were recruited, with at least one member clinically diagnosed with T1DM, as per the American Diabetic Association guidelines. There was a total of fifty individuals, in which 14 were T1DM patients. One family had two children with T1DM, while the rest only had one. Supplementary Figure S1 details the pedigree structure of the participants. All participants were provided with a questionnaire that included details on demographic information and medical history, such as age, gender, weight, height, waist circumference, systolic and diastolic blood pressure, family history of diseases, consanguineous marriages, and smoking status. For the 14 T1DM patients, biochemical tests were carried out, reporting glycosylated hemoglobin (HbA1c), hemoglobin, and white blood cell count (lymphocyte, monocyte). Saliva samples were collected from all subjects using the Oragene OGR-500 kit (DNA Genotek, Ottawa, ON, Canada).

2.2. DNA Extraction and Library Preparation

Genomic DNA was extracted from the buccal cells in the saliva samples using the prepIT®L2P system (DNA Genotek, Ottawa, ON, Canada). The extracted DNA aliquots were quantified using the DS-11 FX Fluorometer (Denovix Inc. Wilmington, DE, USA), and the quality of each sample was assessed through agarose gel electrophoresis. Using the protocol recommended by the manufacturer of the Illumina TruSeq Exome Library Prep kit (Illumina Inc., San Diego, CA, USA), the libraries were prepared from the cleaned and sheared genomic DNA (gDNA). The indexed paired-end libraries were then quantified using the Denovix DS-11 FX Fluorometer to determine the optimal loading concentration of gDNA, providing the adequate clustering density on the flow cell during library sequencing. The fragment size was confirmed using the Advances Analytical Fragment Analyzer (Ankeny, IA, USA). Upon ensuring correct sizes and repeating incorrect runs, the libraries were loaded into NextSeq 500 (Illumina Inc., San Diego, CA, USA) separately.

2.3. Bioinformatics and Data Filtering Pipeline

The data analysis pipeline was designed based on the best practices recommended by the Broad Institute’s Genome Analysis Tool Kit (GATK) instructions, v4.0.6.0 [11]. After subjection to quality assurance, assuming default parameter, the raw reads of samples in the FastQ format files were checked using FastQC software, v0.11.5 [12]. The Trimmomatic tool version 0.33.0 was used to clip all the reads containing Illumina adapters with default parameters for paired-end sequencing. The reads were trimmed from the 3’ end [13]. Using the Burrows-Wheeler Aligner (BWA) v0.7.12 (BWA-MEM), alignment results were generated for each family member by mapping the raw reads to the human reference genome (GRCh37/hg19) [14,15] with reads of 76 base pairs (bps) in length. The sorted BAM files from all the sample lanes were merged, and the duplicate reads were removed using Picard version 2.9.4 tool’s commands SortSam, MergeSamFiles, and MarkDuplicates tools, respectively [16]. The samples’ mapping quality was checked using the Qualimap v2.2.1 tool and verified for an average coverage above 45X [17]. Using BaseRecaliberator from the GATK toolkit, the qualities of the mapped bases were improved, and the variant file from the dbSNP database was used to mark the sites of known variation [18].
The variant calling was performed using Haplotypecaller by GATK in GVCF mode. The results were jointly genotyped using the GenotypeGVCF tool, yielding the final raw variant file [19]. The GATK’s variant recalibration was performed to obtain both SNPs and INDELs before combining them into a single VCF file. Variant annotation was performed using snpEff v4.3p (i.e., frameshift, nonsynonymous, splice variant, etc.). SnpSift was used to annotate the variants using population allele frequencies from gnomAD v2.1.1, 1000 Genomes, and ExAC [20]. The resulting VCF file with the single-nucleotide variants/polymorphisms (SNVs/SNPs) and insertions or deletions (INDELs) that passed the GATK filtering were stored in a GEMINI (GEnome MINIng, v0.12.2) database after decomposing and normalization steps [21]. The variants were then loaded into the GEMINI database with the family’s pedigree file as the input file. The familial analysis was performed with the GEMINI de novo and autosomal recessive commands. Each variant was annotated by comparing it to several genome annotations from the GEMINI database, including the GnomAD, ENCODE tracks, UCSC tracks, OMIM, ClinVar, dbSNP, KEGG, and HPRD. The identified variants were then visually inspected using the IGV tool. The analysis was then focused on SNVs, SNPs, and INDELs that had minor allele frequency < 5% using gnomad_all, which we defined to be nonsynonymous, splice site, or missense SNPs. The impact information comes from snpEff through the GEMINI annotation that internally uses SIFT and PolyPhen to predict functional annotation. Only the variants that were repeated in two or more families were included in the next step of the analysis. Variants were classified as de novo if they were present in the affected child, but not in either of their parents.

2.4. Biological Functional Analysis

The clinical significance, disease associations, and linked phenotypes of the variants were determined using the ClinVar database, GeneCard, NCBI, and online databases [22]. Then, the search converged on genes with relations to type 1 diabetes and all its known pathways and related genes. The genome-wide association study (GWAS) of T1DM was cross-checked with all the variants before and after being filtered [23]. The Human Protein Atlas (http://www.proteinatlas.org (accessed on 1 December 2022)) and the Genotype-Tissue Expression v7 (GTEx) tool were utilized to understand the functional roles, regulatory landscape of gene expression, and splicing variation in a broad selection of primary human tissues. The NHGRI-EBI Catalog of human GWAS (www.ebi.ac.uk/gwas (accessed on 1 December 2022)) was used to identify if an SNP had been identified in a global population, or if they were distinctive to the Emirati population.

3. Results

3.1. Demographic Factors

Of the 50 samples collected from the 13 families, there were 14 probands and 36 unaffected family members. All families had samples from an unaffected mother and unaffected father, and a subset (9 families) had samples of an unaffected sibling (Table 1). The mean age of the probands was 11.86 ± 5.08 years old, while the family members ranged from siblings to parents, and had a mean age of 31.46 ± 14.85 years old. There were nine males and five females in the proband group, whereas for the family members, 19 were male, and 17 were female. The average body mass index (BMI) of the probands was 16.94 ± 2.77, with a waist circumference of 64.92 ± 8.41, whereas the healthy family members had an average BMI of 27.46 ± 5.52, with a waist circumference of 87.82 ± 16.48. Nine out of thirteen families (69.2%) had consanguineous marriages.
Family history of the disease was also collected from the parents, and nine out of thirteen (69.23%) had a history of dyslipidemia, eight (61.54%) had a history of hypertension, five (38.46%) had a history of type 1 diabetes, and two (15.38%) had a history of hypothyroidism. The average HbA1c percentage was 9.12% ± 1.88%, which was higher than the normal range (<6%) for healthy individuals. The rest of the variables (hemoglobin: 13.13 ± 1.40; white cell count: 6.95 ± 4.00; lymphocytes: 5.87± 7.30; monocytes: 2.33 ± 3.58) were all within regular ranges and, thus, had no further significance. The lifestyle of the sample families was assessed for smoking habits, revealing 100% of the probands never smoked, and 25% of the family members never smoked.

3.2. De Novo Variants

Only 1186 variants with high impact passed the filtering pipeline, as detailed in Supplementary Tables S1 and S2. After eliminating the variants with less than 0.05 minor allele frequency (MAF) and variants that were not found in 2 or more families, we had 98 variants. Finally, we filtered out the genes with no biological relevance to T1DM or autoimmune pathways/diseases, and were left with 30 rare variants that belonged to 20 different genes (Table 2). The genes found to be related to T1DM were ten variants across seven genes (MST1, rs201139286; TDG, rs760400700, rs764159587; TYRO3, rs746533465, rs750893216, rs757748573; IFIHI, rs141469634; GLIS3, rs113076411; VEGFA, rs750060813; TYK2, 19-10475177-TA-T). There were twenty variants across thirteen genes that were linked to other forms of diabetes (type 2 diabetes (T2D), gestational diabetes, diabetes retinopathy, and nephrogenic diabetes insipidus) or other autoimmune diseases (BCR, rs372013175; CACNA1B, 9-140773611-G-GACGACACGGAGCCCTATTTCATCGGGATCTT; CNN2, 19-1036442-C-A; COLGALT1, 19-17666649-G-A; LAMA3, 18-21338476-T-G; LGALS9C, rs376412531; MBD4, 3-129155546-CT-C; MST1L, rs11260920; MUC6, rs368342230, rs376177791, rs754249101, rs761220536, rs766751467, rs766833662; PABC1, rs140822921; RNASEH2B, rs200320729; ZNF596, rs756701581).
The MST1L gene located in 1p36.13 is involved in the regulation of macrophage chemotaxis, and mainly expressed in the kidney, liver, and pancreas tissue. This will trigger an inflammatory response, which initiates insulitis and pancreatic cell death, leading to the production of IFNγ, TNFα, IL-1β, and the amplification of beta cell death cycle [24,25]. The IFIH1 gene located in 2q24.2 is an innate immune receptor that plays a major role in sensing viral infection and in the activation of a cascade of antiviral responses, including the induction of type I interferons and proinflammatory cytokines. The IFIH1 gene has been associated with the pathogenesis of diabetes (type 1 and type 2) and multiple autoimmune diseases [23,26,27,28,29,30,31,32,33,34,35,36]. MBD4, located in 3q21.3, is involved in DNA glycosylase and endonuclease activity, and is mainly expressed in lymphocytes. It is associated with obesity, BMI, sclerosis, cancer, and autoimmune disease [37,38,39,40].
The MST1 gene located in 3p21.31 is involved in the regulation of T cell selection, and its deficiency restores normoglycemia, improves beta cell function and prevents the development of diabetes [24,27,28,41]. The VEGFA gene located in 6p21.1 is mainly expressed in thyroid tissue, and has been associated with severe retinopathy in type 1 diabetes, and glomerular microvasculature in diabetes, specifically due to islet vessel density, alteration in expression of genes regulating islet blood flow, insulin deficiency, and inflammation in intra-islet endothelial cells [42,43,44]. The PAPBC1 gene located in 8q22.3 encodes a poly(A) binding protein and is associated with tumor progression [45,46]. The ZNF596 gene located in 8p23.3 is highly expressed in the brain and cerebellum [47]. Characterized by lesions in the central nervous system disseminated in time and space [48], ZNF596 has been reported to be involved in the pathogenesis of multiple sclerosis [49].
The CACNA1B gene located in 9q34.3 is involved in the N-type voltage-dependent calcium channel, which can cause beta cell dysfunction and death, and lead to both types 1 and 2 diabetes [50]. The CACNA1B gene has been associated with acute lymphoblastic leukemia and myeloid leukemia through the regulation of immune functions and leukocyte chemotaxis [51,52]. The GLIS3 gene located in 9p24.2 encodes a nuclear protein that is involved in the expression and development of pancreatic beta cells, and has been associated with neonatal diabetes, fasting blood glucose, type 2 diabetes, and congenital hypothyroidism [53,54,55,56]. GLIS3 plays a role in the generation of pancreatic beta cell viability and susceptibility to immune and metabolic-induced stress, such as proinflammatory cytokines and glucose oxidation [57]. The MUC6 gene located in 11p15.5 is expressed in the stomach and pancreas tissue, is associated with hypertrophic cardiomyopathy, and is involved in enhancing innate immune reactivity [58]. Glycosylation of MUC6 is found to upregulate the IL-17 response, which is found to be related to other immune-mediated diseases [59]. The TYK2 gene located in 19p13.2 is associated with the cytoplasmic domain of type I and type II cytokine receptors, and is a component of both the type I and type III interferon signaling pathways [60,61]. The TYK2 gene has a critical importance in the etiology of autoimmune and inflammatory diseases, specifically type 1 and type 2 diabetes, due to its association with the pancreatic β-cell-specific suppression of cytokine response including IFN [61,62,63,64,65].
The TDG gene located in 12q23.3 is induced by β cells and inflammatory mediators that play a key role in initiating the autoimmune response. Given that TDG enzyme control activates DNA demethylation, TDG gene expression was significantly upregulated in the IFN-α–treated islets and lymphocyte cells [66]. The RNASEH2B gene located in 13q14.3 is involved in the activation of the interferon pathway, leading to the infiltration of lymphocytes and mononuclear cells, and local chronic inflammation [67]. The TYRO3 gene located in 15q15.1 regulates immunoregulation, plays an important role in the inhibition of the Toll-like-receptor-mediated innate immune response, and is an essential regulator of immune homeostasis [68,69]. LGALS9C, located in 17p11.2, is involved in cytoplasmic intracellular functions, and controls AMP-activated protein kinase in response to lysosomal damage, which is caused by diabetes, immune responses, and obesity [70]. This gene is also involved in attenuating T-cell expansion, tumor microenvironment, and chronic infections [71,72]. LAMA3, located in 18q11.2, is associated with the binding to cells via a high-affinity receptor through embryogenesis. While LAMA3 is not associated with T1DM and is mainly expressed in lung tissue, it is associated with autoimmune diseases, immunologic isotypes, immune cytolytic activity, and ovarian cancer [73,74].
The CNN2 gene located in 19p13.3 is involved in the structural organization of actin filaments, playing a role in smooth muscle contraction. While it is not associated with type 1 diabetes, it is expressed in fibroblast cells and linked to the innate immune system, myometrial relaxation, and contraction pathways [75,76]. COLGALT1, located in 19p13.11, encodes collagen β (1-O) galactosyltransferase 1 (ColGalT1), and is associated with musculoskeletal defects, cerebral small vessel disease, and congenital porencephaly [77,78,79]. While the COLGALT1 gene is not associated with T1DM, it is correlated with autoimmune diseases as it could potentially antagonize the innate immune response [80]. The BCR gene in 22q11.23 acts as a GTPase-activating protein that encodes a novel serine/threonine kinase activity, and can be considered as a candidate tumor suppressor gene involved in meningioma pathogenesis and chronic myeloid leukemia [76,81]. The gene is expressed in the brain, endocrine, thyroid, and lymphoid tissue, which may trigger the intracellular signaling pathways leading to the expression of genes required for immune response [82,83].

4. Discussion

This study has identified genetic variants associated with the development of T1DM from 13 Emirati case–parent trios. Of the 13 families sampled for this project, 12 had reported de novo variants, with Family 7 having the highest number (nine) of variants linked to T1DM/autoimmune pathways, whilst Family 4 did not have any variants past the filtering steps. Two of the families (3 and 8) were consanguineously related by the third degree, and were the only families that were associated with the CNN2 gene, linked to the innate immune system and autoimmune diseases.
Gene MUC6 had the highest number of variants (seven) associated with nine families, with two variants present in four families, two variants present in three families, and three variants present in two families. Gene MUC6 has been identified as a loss of function (LoF) variant in the UAE population that is common in the local population (AF > 5%), and rare in global population (AF < 1%, as per the gnomAD exon catalog) [84]. MUC6 has been associated with prostate carcinoma in the UAE population, hence promoting tumorigenesis [84,85]. The gene ontology annotations related to this gene include extracellular matrix structural organization, which may lead to glucose-induced endothelial damage and metabolic disturbances [86]. Sanger sequencing is recommended to identify the association of gene MUC6 with T1DM in the Emirati population.
In comparison to global population, there were ten variants of seven genes reportedly associated with T1DM (MST1, rs201139286; TDG, rs760400700, rs764159587; TYRO3, rs746533465, rs750893216, rs757748573; IFIHI, rs141469634; GLIS3, rs113076411; VEGFA, rs750060813; TYK2, 19-10475177-TA-T). In addition to T1DM, the following genes were associated with other conditions from global GWAS analysis: the MST1 gene was associated with multiple chronic inflammatory disease, including inflammatory bowel disease, Crohn’s disease, psoriasis, and ulcerative colitis [28,87]; the TYRO3 gene was associated with the blood pressure, metabolic syndrome, and microalbuminuria [88,89,90]; the IFIH1 gene was associated with psoriasis, systemic lupus erythematosus, chronic inflammatory diseases, autoimmune thyroid disease, and congenital hypothyroidism [23,26,27,28,29,30,31,32,33,34,35,36]; the GLIS3 gene was associated with type 2 diabetes, chronic obstructive pulmonary disease, and asthma [56,91,92]; and the TYK2 gene was associated with systemic lupus erythematosus, psoriasis, inflammatory bowel disease, rheumatoid arthritis, and type 2 diabetes [30,36,93,94].
While not associated with T1DM, there were twenty variants of thirteen genes that were linked to other forms of diabetes (type 2 diabetes (T2D), gestational diabetes, diabetes retinopathy, and nephrogenic diabetes insipidus) or other autoimmune diseases (BCR, rs372013175; CACNA1B, 9-140773611- G-GACGACACGGAGCCCTATTTCATCGGGATCTT; CNN2, 19-1036442-C-A; COLGALT1, 19-17666649-G-A; LAMA3, 18- 21338476-T-G; LGALS9C, rs376412531; MBD4, 3- 129155546-CT-C; MST1L, rs11260920; MUC6, rs368342230, rs376177791, rs754249101, rs761220536, rs766751467, rs766833662; PABC1, rs140822921; RNASEH2B, rs200320729; ZNF596, rs756701581), when compared to the global population. The following genes were associated with other conditions from GWAS analysis from global databases: the CACNA1B gene was associated with acute myeloid leukemia [51]; the LAMA3 gene was associated with ovarian cancer [74]; the MBD4 gene was associated with systemic sclerosis [38]; the MUC6 gene was associated with hypertrophic cardiomyopathy and peptic ulcer disease [58,95]; and the RNASEH2B gene was associated with rheumatoid arthritis and prostatic hyperplasia [96,97]. Given that T1DM is a multifactorial autoimmune endocrine disease, several susceptible genes may be shared across different conditions [23,98,99,100]. These findings increase our understanding of the genetic contribution and biology underlying T1DM development, and suggest overlapping genetic origins with autoimmune disease and other forms of diabetes.
Several limitations affect the power of trio-based genetic analysis, including locus heterogeneity, sample size, pedigree structure, and genotype accuracy. Given the low sample size, widespread mutational recurrence and heterogeneity may be present, hence a larger sample size is required to confirm the de novo genetic findings. Family 5 had weak coverage and mapping quality in the father’s sample; therefore, variants identified in Family 5 should be considered with caution, which includes rs372013175 of BCR, 19-1036442-C-A of COLGALT1, rs200320729 of RNASEH2B, and rs746533465 of TYRO3. Additional functional studies in cell culture and animal models, as well as re-sequencing the identified genes in larger cohorts, are needed to assess the pathogenicity of the genes to T1DM. Since geographic isolation and consanguinity-driven genomic homozygosity may lead to the enrichment of founder mutations in specific ethnic groups, the presence of such rare mutations in our cohort from consanguineous families is not surprising. Therefore, further studies must be conducted in different ethnic groups to further understand the genetic landscape of T1DM.

5. Conclusions

We have identified de novo genes that may play a role in the pathogenesis of T1DM by conducting whole-exome sequencing on 13 trios from singleton Emirati families with T1DM. As per published GWAS, all the identified genes were either associated with T1DM or were associated with an autoimmune disease. The susceptibility loci for T1DM are heterogeneous, hence further work around identifying such genes may play a robust role in the pathogenesis of T1DM. Future studies can also assess the association of genes with different modes of inheritance, such as autosomal recessive, autosomal dominant, and X-linked. This can aid in understanding the linkage to disease pathogenesis on a wider spectrum, and may offer more relationships between non-HLA genes and T1DM among the citizens of the UAE. Finally, focused PCR-RFLP studies of particularly reported genes from GWAS studies on Emirati subjects would be beneficial to understand the impact of these genes in this region, and would allow the confirmation of variants reported in this study. With the findings presented in this paper, the next step ought to be broadening the spectrum of the search for more disease-causing variants that may be indicative of prediction of type 1 diabetes progression and etiology.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/biology12030413/s1, Figure S1: Pedigree of the 13 Emirati case–parent trios. Table S1: Sequencing coverage, number of reads, and additional metrics for exomes for the proband, mother, father, and sibling individuals. Table S2: All detected SNPs and indels before applying filtering pipeline in each.

Author Contributions

This study was conceived by H.A., S.A. and S.E.H.C. developed and carried out the laboratory assays used in the study. A.A.D. was responsible for the recruitment of the patients’ contributed data to the study by collecting demographic data. S.A. cleaned the data; G.D.E. and M.M. analyzed and constructed the figures. S.A. wrote the first draft of the manuscript. H.A., G.T., S.E.H.C., M.A., Z.F., S.S. and M.M. provided critical review during manuscript preparation. All authors have read and agreed to the published version of the manuscript.

Funding

This research received an internal fund from Khalifa University (KUIRFL1-2015).

Institutional Review Board Statement

This study has been approved by the Institutional Review Board (IRB) committee of Mafraq Hospital (MAF-REC_07/2016_04), in concordance with the Mafraq Hospital Abu Dhabi research ethics committee, the site of sample acquisition. All participating subjects, or the guardians of underage subjects, provided written informed consent after ensuring their full understanding of the project and their contribution and privacy, in accordance with the Declaration of Helsinki. All data were de-identified prior to use.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

All materials will be submitted in the Supplementary Material.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Todd, J.A. Etiology of type 1 diabetes. Immunity 2010, 32, 457–467. [Google Scholar] [CrossRef] [Green Version]
  2. Atkinson, M.A.; Eisenbarth, G.S.; Michels, A.W. Type 1 diabetes. Lancet 2014, 383, 69–82. [Google Scholar] [CrossRef] [Green Version]
  3. Katsarou, A.; Gudbjornsdottir, S.; Rawshani, A.; Dabelea, D.; Bonifacio, E.; Anderson, B.J.; Jacobsen, L.M.; Schatz, D.A.; Lernmark, A. Type 1 diabetes mellitus. Nat. Rev. Dis. Prim. 2017, 3, 17016. [Google Scholar] [CrossRef]
  4. Haller, M.J.; Atkinson, M.A.; Schatz, D. Type 1 diabetes mellitus: Etiology, presentation, and management. Pediatric. Clin. 2005, 52, 1553–1578. [Google Scholar] [CrossRef]
  5. Craig, M.E.; Kim, K.W.; Isaacs, S.R.; Penno, M.A.; Hamilton-Williams, E.E.; Couper, J.J.; Rawlinson, W.D. Early-life factors contributing to type 1 diabetes. Diabetologia 2019, 62, 1823–1834. [Google Scholar] [CrossRef] [Green Version]
  6. Olmos, P.; A’hern, R.; Heaton, D.; Millward, B.; Risley, D.; Pyke, D.; Leslie, R. The significance of the concordance rate for type 1 (insulin-dependent) diabetes in identical twins. Diabetologia 1988, 31, 747–750. [Google Scholar] [CrossRef] [Green Version]
  7. Levitsky, L.L.; Misra, M. Epidemiology, Presentation, and Diagnosis of Type 1 Diabetes Mellitus in Children and Adolescents; UpToDate: Waltham, MA, USA, 2007. [Google Scholar]
  8. Redondo, M.J.; Steck, A.K.; Pugliese, A. Genetics of type 1 diabetes. Pediatric. Diabetes 2018, 19, 346–353. [Google Scholar] [CrossRef]
  9. Nyaga, D.M.; Vickers, M.H.; Jefferies, C.; Perry, J.K.; O’Sullivan, J.M. The genetic architecture of type 1 diabetes mellitus. Mol. Cell. Endocrinol. 2018, 477, 70–80. [Google Scholar] [CrossRef] [PubMed]
  10. Lee, H.S.; Hwang, J.S. Genetic aspects of type 1 diabetes. Ann. Pediatr. Endocrinol. Metab. 2019, 24, 143. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. DePristo, M.A.; Banks, E.; Poplin, R.; Garimella, K.V.; Maguire, J.R.; Hartl, C.; Philippakis, A.A.; Del Angel, G.; Rivas, M.A.; Hanna, M. A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat. Genet. 2011, 43, 491. [Google Scholar] [CrossRef] [PubMed]
  12. Andrews, S. FastQC: A Quality Control Tool for High Throughput Sequence Data; Babraham Bioinformatics, Babraham Institute: Cambridge, UK, 2010. [Google Scholar]
  13. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef] [Green Version]
  14. Fujita, P.A.; Rhead, B.; Zweig, A.S.; Hinrichs, A.S.; Karolchik, D.; Cline, M.S.; Goldman, M.; Barber, G.P.; Clawson, H.; Coelho, A. The UCSC genome browser database: Update 2011. Nucleic Acids Res. 2010, 39 (Suppl. S1), D876–D882. [Google Scholar] [CrossRef] [Green Version]
  15. Lander, E.S.; Linton, L.M.; Birren, B.; Nusbaum, C.; Zody, M.C.; Baldwin, J.; Devon, K.; Dewar, K.; Doyle, M.; FitzHugh, W.; et al. Initial sequencing and analysis of the human genome. Nature 2001, 409, 860–921. [Google Scholar] [CrossRef] [Green Version]
  16. Lee, D.-J.; Kwon, T.; Kim, C.-K.; Seol, Y.-J.; Park, D.-S.; Lee, T.-H.; Ahn, B.-O. NGS_SNPAnalyzer: A desktop software supporting genome projects by identifying and visualizing sequence variations from next-generation sequencing data. Genes Genom. 2020, 42, 1311–1317. [Google Scholar] [CrossRef]
  17. McKenna, A.; Hanna, M.; Banks, E.; Sivachenko, A.; Cibulskis, K.; Kernytsky, A.; Garimella, K.; Altshuler, D.; Gabriel, S.; Daly, M. The Genome Analysis Toolkit: A MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010, 20, 1297–1303. [Google Scholar] [CrossRef] [Green Version]
  18. GATK Team. BaseRecalibrator. 2020. Available online: https://gatk.broadinstitute.org/hc/en-us/articles/360036898312-BaseRecalibrator (accessed on 1 January 2021).
  19. Poplin, R.; Ruano-Rubio, V.; DePristo, M.A.; Fennell, T.J.; Carneiro, M.O.; Van der Auwera, G.A.; Kling, D.E.; Gauthier, L.D.; Levy-Moonshine, A.; Roazen, D. Scaling accurate genetic variant discovery to tens of thousands of samples. bioRxiv, 2017; bioRxiv:201178. [Google Scholar] [CrossRef] [Green Version]
  20. Liu, X.; Jian, X.; Boerwinkle, E. dbNSFP: A lightweight database of human nonsynonymous SNPs and their functional predictions. Hum. Mutat. 2011, 32, 894–899. [Google Scholar] [CrossRef]
  21. Paila, U.; Chapman, B.A.; Kirchner, R.; Quinlan, A.R. GEMINI: Integrative exploration of genetic variation and genome annotations. PLoS Comput. Biol. 2013, 9, e1003153. [Google Scholar] [CrossRef]
  22. Stelzer, G.; Rosen, N.; Plaschkes, I.; Zimmerman, S.; Twik, M.; Fishilevich, S.; Stein, T.I.; Nudel, R.; Lieder, I.; Mazor, Y. The GeneCards suite: From gene data mining to disease genome sequence analyses. Curr. Protoc. Bioinform. 2016, 54, 1–30. [Google Scholar] [CrossRef]
  23. Barrett, J.C.; Clayton, D.G.; Concannon, P.; Akolkar, B.; Cooper, J.D.; Erlich, H.A.; Julier, C.; Morahan, G.; Nerup, J.; Nierras, C.; et al. Genome-wide association study and meta-analysis find that over 40 loci affect risk of type 1 diabetes. Nat. Genet. 2009, 41, 703–707. [Google Scholar] [CrossRef] [Green Version]
  24. Ardestani, A.; Paroni, F.; Azizi, Z.; Kaur, S.; Khobragade, V.; Yuan, T.; Frogne, T.; Tao, W.; Oberholzer, J.; Pattou, F. MST1 is a key regulator of beta cell apoptosis and dysfunction in diabetes. Nat. Med. 2014, 20, 385–397. [Google Scholar] [CrossRef] [Green Version]
  25. Burrack, A.L.; Martinov, T.; Fife, B.T. T cell-mediated beta cell destruction: Autoimmunity and alloimmunity in the context of type 1 diabetes. Front. Endocrinol. 2017, 8, 343. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Saevarsdottir, S.; Olafsdottir, T.A.; Ivarsdottir, E.V.; Halldorsson, G.H.; Gunnarsdottir, K.; Sigurdsson, A.; Johannesson, A.; Sigurdsson, J.K.; Juliusdottir, T.; Lund, S.H. FLT3 stop mutation increases FLT3 ligand level and risk of autoimmune thyroid disease. Nature 2020, 584, 619–623. [Google Scholar] [CrossRef] [PubMed]
  27. Jostins, L.; Ripke, S.; Weersma, R.K.; Duerr, R.H.; McGovern, D.P.; Hui, K.Y.; Lee, J.C.; Philip Schumm, L.; Sharma, Y.; Anderson, C.A. Host–microbe interactions have shaped the genetic architecture of inflammatory bowel disease. Nature 2012, 491, 119–124. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Ellinghaus, D.; Jostins, L.; Spain, S.L.; Cortes, A.; Bethune, J.; Han, B.; Park, Y.R.; Raychaudhuri, S.; Pouget, J.G.; Hübenthal, M. Analysis of five chronic inflammatory diseases identifies 27 new associations and highlights disease-specific patterns at shared loci. Nat. Genet. 2016, 48, 510–518. [Google Scholar] [CrossRef] [PubMed]
  29. Jin, Y.; Andersen, G.; Yorgov, D.; Ferrara, T.M.; Ben, S.; Brownson, K.M.; Holland, P.J.; Birlea, S.A.; Siebert, J.; Hartmann, A. Genome-wide association studies of autoimmune vitiligo identify 23 new risk loci and highlight key pathways and regulatory variants. Nat. Genet. 2016, 48, 1418–1424. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Bentham, J.; Morris, D.L.; Cunninghame Graham, D.S.; Pinder, C.L.; Tombleson, P.; Behrens, T.W.; Martín, J.; Fairfax, B.P.; Knight, J.C.; Chen, L. Genetic association analyses implicate aberrant regulation of innate and adaptive immunity genes in the pathogenesis of systemic lupus erythematosus. Nat. Genet. 2015, 47, 1457–1464. [Google Scholar] [CrossRef]
  31. Langefeld, C.D.; Ainsworth, H.C.; Graham, D.S.C.; Kelly, J.A.; Comeau, M.E.; Marion, M.C.; Howard, T.D.; Ramos, P.S.; Croker, J.A.; Morris, D.L. Transancestral mapping and genetic load in systemic lupus erythematosus. Nat. Commun. 2017, 8, 16021. [Google Scholar] [CrossRef] [Green Version]
  32. Onengut-Gumuscu, S.; Chen, W.-M.; Burren, O.; Cooper, N.J.; Quinlan, A.R.; Mychaleckyj, J.C.; Farber, E.; Bonnie, J.K.; Szpak, M.; Schofield, E. Fine mapping of type 1 diabetes susceptibility loci and evidence for colocalization of causal variants with lymphoid gene enhancers. Nat. Genet. 2015, 47, 381–386. [Google Scholar] [CrossRef]
  33. Tsoi, L.C.; Spain, S.L.; Knight, J.; Ellinghaus, E.; Stuart, P.E.; Capon, F.; Ding, J.; Li, Y.; Tejasvi, T.; Gudjonsson, J.E. Identification of 15 new psoriasis susceptibility loci highlights the role of innate immunity. Nat. Genet. 2012, 44, 1341–1348. [Google Scholar] [CrossRef] [Green Version]
  34. Plagnol, V.; Howson, J.M.; Smyth, D.J.; Walker, N.; Hafler, J.P.; Wallace, C.; Stevens, H.; Jackson, L.; Simmonds, M.J.; Consortium, T.D.G. Genome-wide association analysis of autoantibody positivity in type 1 diabetes cases. PLoS Genet. 2011, 7, e1002216. [Google Scholar] [CrossRef]
  35. Yin, X.; Low, H.Q.; Wang, L.; Li, Y.; Ellinghaus, E.; Han, J.; Estivill, X.; Sun, L.; Zuo, X.; Shen, C. Genome-wide meta-analysis identifies multiple novel associations and ethnic heterogeneity of psoriasis susceptibility. Nat. Commun. 2015, 6, 6916. [Google Scholar] [CrossRef] [Green Version]
  36. Patrick, M.T.; Stuart, P.E.; Zhang, H.; Zhao, Q.; Yin, X.; He, K.; Zhou, X.-J.; Mehta, N.N.; Voorhees, J.J.; Boehnke, M. Causal relationship and shared genetic loci between psoriasis and type 2 diabetes through trans-disease meta-analysis. J. Investig. Dermatol. 2021, 141, 1493–1502. [Google Scholar] [CrossRef]
  37. Christakoudi, S.; Evangelou, E.; Riboli, E.; Tsilidis, K.K. GWAS of allometric body-shape indices in UK Biobank identifies loci suggesting associations with morphogenesis, organogenesis, adrenal cell renewal and cancer. Sci. Rep. 2021, 11, 10688. [Google Scholar] [CrossRef]
  38. Gorlova, O.Y.; Li, Y.; Gorlov, I.; Ying, J.; Chen, W.V.; Assassi, S.; Reveille, J.D.; Arnett, F.C.; Zhou, X.; Bossini-Castillo, L. Gene-level association analysis of systemic sclerosis: A comparison of African-Americans and White populations. PLoS ONE 2018, 13, e0189498. [Google Scholar] [CrossRef] [Green Version]
  39. Bellacosa, A. Role of MED1 (MBD4) gene in DNA repair and human cancer. J. Cell. Physiol. 2001, 187, 137–144. [Google Scholar] [CrossRef]
  40. Liao, W.; Li, M.; Wu, H.; Jia, S.; Zhang, N.; Dai, Y.; Zhao, M.; Lu, Q. Down-regulation of MBD4 contributes to hypomethylation and overexpression of CD70 in CD4+ T cells in systemic lupus erythematosus. Clin. Epigenetics 2017, 9, 104. [Google Scholar] [CrossRef] [Green Version]
  41. Gudjonsson, A.; Gudmundsdottir, V.; Axelsson, G.T.; Gudmundsson, E.F.; Jonsson, B.G.; Launer, L.J.; Lamb, J.R.; Jennings, L.L.; Aspelund, T.; Emilsson, V. A genome-wide association study of serum proteins reveals shared loci with common diseases. Nat. Commun. 2022, 13, 480. [Google Scholar] [CrossRef]
  42. Sivaskandarajah, G.A.; Jeansson, M.; Maezawa, Y.; Eremina, V.; Baelde, H.J.; Quaggin, S.E. Vegfa protects the glomerular microvasculature in diabetes. Diabetes 2012, 61, 2958–2966. [Google Scholar] [CrossRef] [Green Version]
  43. Al-Kateb, H.; Mirea, L.; Xie, X.; Sun, L.; Liu, M.; Chen, H.; Bull, S.B.; Boright, A.P.; Paterson, A.D.; Group, D.E.R. Multiple variants in vascular endothelial growth factor (VEGFA) are risk factors for time to severe retinopathy in type 1 diabetes: The DCCT/EDIC genetics study. Diabetes 2007, 56, 2161–2168. [Google Scholar] [CrossRef] [Green Version]
  44. Staels, W.; Heremans, Y.; Heimberg, H.; De Leu, N. VEGF-A and blood vessels: A beta cell perspective. Diabetologia 2019, 62, 1961–1968. [Google Scholar] [CrossRef] [Green Version]
  45. Stupfler, B.; Birck, C.; Séraphin, B.; Mauxion, F. BTG2 bridges PABPC1 RNA-binding domains and CAF1 deadenylase to control cell proliferation. Nat. Commun. 2016, 7, 480. [Google Scholar] [CrossRef] [Green Version]
  46. Takashima, N.; Ishiguro, H.; Kuwabara, Y.; Kimura, M.; Haruki, N.; Ando, T.; Kurehara, H.; Sugito, N.; Mori, R.; Fujii, Y. Expression and prognostic roles of PABPC1 in esophageal cancer: Correlation with tumor progression and postoperative survival. Oncol. Rep. 2006, 15, 667–671. [Google Scholar] [CrossRef]
  47. Chi, C.; Shao, X.; Rhead, B.; Gonzales, E.; Smith, J.B.; Xiang, A.H.; Graves, J.; Waldman, A.; Lotze, T.; Schreiner, T. Admixture mapping reveals evidence of differential multiple sclerosis risk by genetic ancestry. PLoS Genet. 2019, 15, e1007808. [Google Scholar] [CrossRef]
  48. Grigoriadis, N.; Van Pesch, V. A basic overview of multiple sclerosis immunopathology. Eur. J. Neurol. 2015, 22, 3–13. [Google Scholar] [CrossRef]
  49. Ventura, R.E.; Antezana, A.O.; Bacon, T.; Kister, I. Hispanic Americans and African Americans with multiple sclerosis have more severe disease course than Caucasian Americans. Mult. Scler. J. 2017, 23, 1554–1557. [Google Scholar] [CrossRef]
  50. Yang, S.-N.; Berggren, P.-O. The role of voltage-gated calcium channels in pancreatic β-cell physiology and pathophysiology. Endocr. Rev. 2006, 27, 621–676. [Google Scholar] [CrossRef] [Green Version]
  51. Fernandez, C.A.; Smith, C.; Yang, W.; Mullighan, C.G.; Qu, C.; Larsen, E.; Bowman, W.P.; Liu, C.; Ramsey, L.B.; Chang, T. Genome-wide analysis links NFATC2 with asparaginase hypersensitivity. Blood J. Am. Soc. Hematol. 2015, 126, 69–75. [Google Scholar] [CrossRef] [Green Version]
  52. Lv, H.; Zhang, M.; Shang, Z.; Li, J.; Zhang, S.; Lian, D.; Zhang, R. Genome-wide haplotype association study identify the FGFR2 gene as a risk gene for acute myeloid leukemia. Oncotarget 2017, 8, 7891. [Google Scholar] [CrossRef] [Green Version]
  53. Grant, S.F.; Qu, H.-Q.; Bradfield, J.P.; Marchand, L.; Kim, C.E.; Glessner, J.T.; Grabs, R.; Taback, S.P.; Frackelton, E.C.; Eckert, A.W. Follow-up analysis of genome-wide association data identifies novel loci for type 1 diabetes. Diabetes 2009, 58, 290–295. [Google Scholar] [CrossRef] [Green Version]
  54. Wessel, J.; Chu, A.Y.; Willems, S.M.; Wang, S.; Yaghootkar, H.; Brody, J.A.; Dauriz, M.; Hivert, M.-F.; Raghavan, S.; Lipovich, L. Low-frequency and rare exome chip variants associate with fasting glucose and type 2 diabetes susceptibility. Nat. Commun. 2015, 6, 5897. [Google Scholar] [CrossRef] [Green Version]
  55. Li, H.; Gan, W.; Lu, L.; Dong, X.; Han, X.; Hu, C.; Yang, Z.; Sun, L.; Bao, W.; Li, P. A genome-wide association study identifies GRK5 and RASGRP1 as type 2 diabetes loci in Chinese Hans. Diabetes 2013, 62, 291–298. [Google Scholar] [CrossRef] [Green Version]
  56. Wang, X.; Jia, J.; Huang, T. Shared genetic architecture and casual relationship between leptin levels and type 2 diabetes: Large-scale cross-trait meta-analysis and Mendelian randomization analysis. BMJ Open Diabetes Res. Care 2020, 8, e001140. [Google Scholar] [CrossRef]
  57. Nogueira, T.C.; Paula, F.M.; Villate, O.; Colli, M.L.; Moura, R.F.; Cunha, D.A.; Marselli, L.; Marchetti, P.; Cnop, M.; Julier, C. GLIS3, a susceptibility gene for type 1 and type 2 diabetes, modulates pancreatic beta cell apoptosis via regulation of a splice variant of the BH3-only protein Bim. PLoS Genet. 2013, 9, e1003532. [Google Scholar] [CrossRef]
  58. Harper, A.R.; Goel, A.; Grace, C.; Thomson, K.L.; Petersen, S.E.; Xu, X.; Waring, A.; Ormondroyd, E.; Kramer, C.M.; Ho, C.Y. Common genetic variants and modifiable risk factors underpin hypertrophic cardiomyopathy susceptibility and expressivity. Nat. Genet. 2021, 53, 135–142. [Google Scholar] [CrossRef]
  59. Honkanen, J.; Nieminen, J.K.; Gao, R.; Luopajarvi, K.; Salo, H.M.; Ilonen, J.; Knip, M.; Otonkoski, T.; Vaarala, O. IL-17 immunity in human type 1 diabetes. J. Immunol. 2010, 185, 1959–1967. [Google Scholar] [CrossRef] [Green Version]
  60. Sliz, E.; Kalaoja, M.; Ahola-Olli, A.; Raitakari, O.; Perola, M.; Salomaa, V.; Lehtimäki, T.; Karhu, T.; Viinamäki, H.; Salmi, M. Genome-wide association study identifies seven novel loci associating with circulating cytokines and cell adhesion molecules in Finns. J. Med. Genet. 2019, 56, 607–616. [Google Scholar] [CrossRef] [Green Version]
  61. Tao, J.-H.; Zou, Y.-F.; Feng, X.-L.; Li, J.; Wang, F.; Pan, F.-M.; Ye, D.-Q. Meta-analysis of TYK2 gene polymorphisms association with susceptibility to autoimmune and inflammatory diseases. Mol. Biol. Rep. 2011, 38, 4663–4672. [Google Scholar] [CrossRef]
  62. Marroqui, L.; Dos Santos, R.S.; Fløyel, T.; Grieco, F.A.; Santin, I.; Op de Beeck, A.; Marselli, L.; Marchetti, P.; Pociot, F.; Eizirik, D.L. TYK2, a candidate gene for type 1 diabetes, modulates apoptosis and the innate immune response in human pancreatic β-cells. Diabetes 2015, 64, 3808–3817. [Google Scholar] [CrossRef] [Green Version]
  63. Mori, H.; Takahashi, H.; Mine, K.; Higashimoto, K.; Inoue, K.; Kojima, M.; Kuroki, S.; Eguchi, T.; Ono, Y.; Inuzuka, S. TYK2 Promoter Variant Is Associated with Impaired Insulin Secretion and Lower Insulin Resistance in Japanese Type 2 Diabetes Patients. Genes 2021, 12, 400. [Google Scholar] [CrossRef]
  64. Izumi, K.; Mine, K.; Inoue, Y.; Teshima, M.; Ogawa, S.; Kai, Y.; Kurafuji, T.; Hirakawa, K.; Miyakawa, D.; Ikeda, H. Reduced Tyk2 gene expression in β-cells due to natural mutation determines susceptibility to virus-induced diabetes. Nat. Commun. 2015, 6, 6748. [Google Scholar] [CrossRef] [Green Version]
  65. Størling, J.; Pociot, F. Type 1 diabetes candidate genes linked to pancreatic islet cell inflammation and beta-cell apoptosis. Genes 2017, 8, 72. [Google Scholar] [CrossRef] [Green Version]
  66. Almgren, P.; Lindqvist, A.; Krus, U.; Hakaste, L.; Ottosson-Laakso, E.; Asplund, O.; Sonestedt, E.; Prasad, R.B.; Laurila, E.; Orho-Melander, M. Genetic determinants of circulating GIP and GLP-1 concentrations. JCI Insight 2017, 2, e93306. [Google Scholar] [CrossRef] [Green Version]
  67. He, T.; Xia, Y.; Yang, J. Systemic inflammation and chronic kidney disease in a patient due to the RNASEH2B defect. Pediatr. Rheumatol. 2021, 19, 9. [Google Scholar] [CrossRef]
  68. Ji, R.; Meng, L.; Li, Q.; Lu, Q. TAM receptor deficiency affects adult hippocampal neurogenesis. Metab. Brain Dis. 2015, 30, 633–644. [Google Scholar] [CrossRef] [Green Version]
  69. Paolino, M.; Penninger, J.M. The role of TAM family receptors in immune cell function: Implications for cancer therapy. Cancers 2016, 8, 97. [Google Scholar] [CrossRef] [Green Version]
  70. Jia, J.; Bissa, B.; Brecht, L.; Allers, L.; Choi, S.W.; Gu, Y.; Zbinden, M.; Burge, M.R.; Timmins, G.; Hallows, K. AMPK, a regulator of metabolism and autophagy, is activated by lysosomal damage via a novel galectin-directed ubiquitin signal transduction system. Mol. Cell 2020, 77, 951–969.e959. [Google Scholar] [CrossRef]
  71. Sakuishi, K.; Apetoh, L.; Sullivan, J.M.; Blazar, B.R.; Kuchroo, V.K.; Anderson, A.C. Targeting Tim-3 and PD-1 pathways to reverse T cell exhaustion and restore anti-tumor immunity. J. Exp. Med. 2010, 207, 2187–2194. [Google Scholar] [CrossRef]
  72. Silva, I.G.; Yasinska, I.M.; Sakhnevych, S.S.; Fiedler, W.; Wellbrock, J.; Bardelli, M.; Varani, L.; Hussain, R.; Siligardi, G.; Ceccone, G. The Tim-3-galectin-9 secretory pathway is involved in the immune escape of human acute myeloid leukemia cells. EBioMedicine 2017, 22, 44–57. [Google Scholar] [CrossRef] [Green Version]
  73. Jonsson, S.; Sveinbjornsson, G.; de Lapuente Portilla, A.L.; Swaminathan, B.; Plomp, R.; Dekkers, G.; Ajore, R.; Ali, M.; Bentlage, A.E.; Elmer, E. Identification of sequence variants influencing immunoglobulin levels. Nat. Genet. 2017, 49, 1182–1191. [Google Scholar] [CrossRef]
  74. Phelan, C.M.; Kuchenbaecker, K.B.; Tyrer, J.P.; Kar, S.P.; Lawrenson, K.; Winham, S.J.; Dennis, J.; Pirie, A.; Riggan, M.J.; Chornokur, G. Identification of 12 new susceptibility loci for different histotypes of epithelial ovarian cancer. Nat. Genet. 2017, 49, 680–691. [Google Scholar] [CrossRef] [Green Version]
  75. Astle, W.J.; Elding, H.; Jiang, T.; Allen, D.; Ruklisa, D.; Mann, A.L.; Mead, D.; Bouman, H.; Riveros-Mckay, F.; Kostadima, M.A. The allelic landscape of human blood cell trait variation and links to common complex disease. Cell 2016, 167, 1415–1429.e1419. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Chen, M.-H.; Raffield, L.M.; Mousas, A.; Sakaue, S.; Huffman, J.E.; Moscati, A.; Trivedi, B.; Jiang, T.; Akbari, P.; Vuckovic, D. Trans-ethnic and ancestry-specific blood-cell genetics in 746,667 individuals from 5 global populations. Cell 2020, 182, 1198–1213.e1114. [Google Scholar] [CrossRef] [PubMed]
  77. Miyatake, S.; Schneeberger, S.; Koyama, N.; Yokochi, K.; Ohmura, K.; Shiina, M.; Mori, H.; Koshimizu, E.; Imagawa, E.; Uchiyama, Y. Biallelic COLGALT1 variants are associated with cerebral small vessel disease. Ann. Neurol. 2018, 84, 843–853. [Google Scholar] [CrossRef]
  78. Geister, K.A.; Lopez-Jimenez, A.J.; Houghtaling, S.; Ho, T.-H.; Vanacore, R.; Beier, D.R. Loss of function of Colgalt1 disrupts collagen post-translational modification and causes musculoskeletal defects. Dis. Model. Mech. 2019, 12, dmm037176. [Google Scholar] [CrossRef] [Green Version]
  79. Teunissen, M.W.; Kamsteeg, E.-J.; Sallevelt, S.C.; Pennings, M.; Bauer, N.J.; Vermeulen, R.J.; Nicolai, J. Biallelic variants in the COLGALT1 gene causes severe congenital porencephaly: A case report. Neurol. Genet. 2021, 7, e564. [Google Scholar] [CrossRef] [PubMed]
  80. Han, B.; Liu, S.; Xing, Q.; Yu, Y.; Wang, Y. COLGALT1 is a Prognostic Biomarker in Clear Renal Cell Carcinoma Correlated with Immune Infiltrates: A Study Based on TCGA Data. Res. Sq. 2020, 1–24. [Google Scholar] [CrossRef]
  81. Vuckovic, D.; Bao, E.L.; Akbari, P.; Lareau, C.A.; Mousas, A.; Jiang, T.; Chen, M.-H.; Raffield, L.M.; Tardaguila, M.; Huffman, J.E. The polygenic and monogenic basis of blood traits and diseases. Cell 2020, 182, 1214–1231.e1211. [Google Scholar] [CrossRef]
  82. Cooper, M.D.; Alder, M.N. The evolution of adaptive immune systems. Cell 2006, 124, 815–822. [Google Scholar] [CrossRef] [Green Version]
  83. Repo, P.; Jäntti, J.; Järvinen, R.S.; Täll, M.; Raivio, V.; Kivelä, T.; Turunen, J. Germline variants in the MBD4 gene are rare in patients with uveal melanoma. Acta Ophthalmol. 2019, 97, 343. [Google Scholar] [CrossRef]
  84. Daw Elbait, G.; Henschel, A.; Tay, G.K.; Al Safar, H.S. A population-specific major allele reference genome from the United Arab Emirates population. Front. Genet. 2021, 12, 660428. [Google Scholar] [CrossRef]
  85. Liu, B.-H.; Liu, G.-B.; Zhang, B.-B.; Shen, J.; Xie, L.-L.; Liu, X.-Q.; Yao, W.; Dong, R.; Bi, Y.-L.; Dong, K.-R. Tumor Suppressive Role of MUC6 in Wilms Tumor via Autophagy-Dependent β-Catenin Degradation. Front. Oncol. 2022, 12, 756117. [Google Scholar] [CrossRef]
  86. Biswas, S.; Chakrabarti, S. Increased extracellular matrix protein production in chronic diabetic complications: Implications of non-coding RNAs. Non-Coding RNA 2019, 5, 30. [Google Scholar] [CrossRef] [Green Version]
  87. De Lange, K.M.; Moutsianas, L.; Lee, J.C.; Lamb, C.A.; Luo, Y.; Kennedy, N.A.; Jostins, L.; Rice, D.L.; Gutierrez-Achury, J.; Ji, S.-G. Genome-wide association study implicates immune activation of multiple integrin genes in inflammatory bowel disease. Nat. Genet. 2017, 49, 256–261. [Google Scholar] [CrossRef] [Green Version]
  88. Warren, H.R.; Evangelou, E.; Cabrera, C.P.; Gao, H.; Ren, M.; Mifsud, B.; Ntalla, I.; Surendran, P.; Liu, C.; Cook, J.P. Genome-wide association analysis identifies novel blood pressure loci and offers biological insights into cardiovascular risk. Nat. Genet. 2017, 49, 403–415. [Google Scholar] [CrossRef]
  89. Lind, L. Genome-wide association study of the metabolic syndrome in UK biobank. Metab. Syndr. Relat. Disord. 2019, 17, 505–511. [Google Scholar] [CrossRef]
  90. Teumer, A.; Li, Y.; Ghasemi, S.; Prins, B.P.; Wuttke, M.; Hermle, T.; Giri, A.; Sieber, K.B.; Qiu, C.; Kirsten, H. Genome-wide association meta-analyses and fine-mapping elucidate pathways influencing albuminuria. Nat. Commun. 2019, 10, 4130. [Google Scholar] [CrossRef] [Green Version]
  91. Cho, Y.S.; Chen, C.-H.; Hu, C.; Long, J.; Hee Ong, R.T.; Sim, X.; Takeuchi, F.; Wu, Y.; Go, M.J.; Yamauchi, T. Meta-analysis of genome-wide association studies identifies eight new loci for type 2 diabetes in east Asians. Nat. Genet. 2012, 44, 67–72. [Google Scholar] [CrossRef]
  92. Kim, W.; Prokopenko, D.; Sakornsakolpat, P.; Hobbs, B.D.; Lutz, S.M.; Hokanson, J.E.; Wain, L.V.; Melbourne, C.A.; Shrine, N.; Tobin, M.D. Genome-wide gene-by-smoking interaction study of chronic obstructive pulmonary disease. Am. J. Epidemiol. 2021, 190, 875–885. [Google Scholar] [CrossRef]
  93. Sakaue, S.; Kanai, M.; Tanigawa, Y.; Karjalainen, J.; Kurki, M.; Koshiba, S.; Narita, A.; Konuma, T.; Yamamoto, K.; Akiyama, M. A cross-population atlas of genetic associations for 220 human phenotypes. Nat. Genet. 2021, 53, 1415–1424. [Google Scholar] [CrossRef]
  94. Acosta-Herrera, M.; Kerick, M.; González-Serna, D.; Wijmenga, C.; Franke, A.; Gregersen, P.K.; Padyukov, L.; Worthington, J.; Vyse, T.J.; Alarcón-Riquelme, M.E. Genome-wide meta-analysis reveals shared new loci in systemic seropositive rheumatic diseases. Ann. Rheum. Dis. 2019, 78, 311–319. [Google Scholar] [CrossRef] [Green Version]
  95. Wu, Y.; Murray, G.K.; Byrne, E.M.; Sidorenko, J.; Visscher, P.M.; Wray, N.R. GWAS of peptic ulcer disease implicates Helicobacter pylori infection, other gastrointestinal disorders and depression. Nat. Commun. 2021, 12, 1146. [Google Scholar] [CrossRef] [PubMed]
  96. Gudmundsson, J.; Sigurdsson, J.K.; Stefansdottir, L.; Agnarsson, B.A.; Isaksson, H.J.; Stefansson, O.A.; Gudjonsson, S.A.; Gudbjartsson, D.F.; Masson, G.; Frigge, M.L. Genome-wide associations for benign prostatic hyperplasia reveal a genetic correlation with serum levels of PSA. Nat. Commun. 2018, 9, 4568. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  97. Bossini-Castillo, L.; De Kovel, C.; Kallberg, H.; van‘t Slot, R.; Italiaander, A.; Coenen, M.; Tak, P.; Posthumus, M.; Wijmenga, C.; Huizinga, T. A genome-wide association study of rheumatoid arthritis without antibodies against citrullinated peptides. Ann. Rheum. Dis. 2015, 74, e15. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  98. Barker, J.M. Type 1 diabetes-associated autoimmunity: Natural history, genetic associations, and screening. J. Clin. Endocrinol. Metab. 2006, 91, 1210–1217. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  99. Fung, E.; Smyth, D.J.; Howson, J.M.; Cooper, J.D.; Walker, N.M.; Stevens, H.; Wicker, L.S.; Todd, J.A. Analysis of 17 autoimmune disease-associated variants in type 1 diabetes identifies 6q23/TNFAIP3 as a susceptibility locus. Genes Immun. 2009, 10, 188–191. [Google Scholar] [CrossRef] [PubMed]
  100. Becker, K.G. Comparative genetics of type 1 diabetes and autoimmune disease: Common loci, common pathways? Diabetes 1999, 48, 1353–1358. [Google Scholar] [CrossRef]
Table 1. Demographic and anthropometric data of probands and healthy family members.
Table 1. Demographic and anthropometric data of probands and healthy family members.
VariableT1DM (N = 14)Healthy Family Members (N = 36)
Age (years)11.86 ± 5.0831.46 ± 14.85
Gender
Male9 (64.29%)19 (52.78%)
Female5 (35.71%)17 (47.22)
BMI (kg/m2)16.94 ± 2.7727.46 ± 5.52
Waist circumference (cm)64.92 ± 8.4187.82 ± 16.48
Systolic blood pressure (mmHg)108.43 ± 10.89119.29 ± 14.15
Diastolic blood pressure (mmHg)69.79 ± 7.8574.14 ± 10.17
Continuous variables were recorded as mean ± standard deviation; categorical variables were recorded as number of events (percentages %); BMI: body mass index.
Table 2. De novo genes with rare variants identified from 13 Emirati case–parent trios with T1DM.
Table 2. De novo genes with rare variants identified from 13 Emirati case–parent trios with T1DM.
ChrGeneImpactrs-IDFamily IDAFFamily GenotypeLink to T1DMLink to Auto-Immune Diseases
1MST1LStop Gainedrs11260920Family 1 + 125.51 × 10−2G/A,G/G,G/G,G/G,G/GNoYes
2IFIH1Intron Variantrs141469634Family 8 + 97.79 × 10−3CA/C,CA/CA,CA/CA,CA/CAYesYes
3MBD4Frameshift Variantrs747480541Family 2 + 7 + 113.11 × 10−3CT/C,CT/CT,CT/CTNoYes
3MST1Splice Acceptor Variantrs201139286Family 3 + 138.15 × 10−3T/G,T/T,T/T,T/TYesNo
6VEGFAIntron Variantrs750060813Family 5 + 77.24 × 10−4T/C,T/T,T/TYesYes
8PABPC1Frameshift Variantrs140822921Family 3 + 91.15 × 10−3C/CA,C/C,C/CNoYes
8ZNF596Stop Gainedrs756701581Family 2 + 73.98 × 10−5C/A,C/C,C/C,C/CNoYes
9CACNA1BSplice Donor VariantNoneFamily 1 + 9 + 109.04 × 10−3G/GACGACACGGAGCCCTATTTCATCGGGATCTT,G/G,G/G,G/GNoYes
9GLIS3Synonymous Variantrs113076411Family 8 + 104.80 × 10−3CT/C,CT/CT,CT/CT,CT/CTYesYes
11MUC6Frameshift Variantrs368342230Family 3 + 7 + 84.11 × 10−6TG/T,TG/TG,TG/TG,TG/TGNoYes
rs376177791Family 3 + 7 + 81.01 × 10−5G/GT,G/G,G/G,G/G
rs754249101Family 6 + 91.87 × 10−4T/TGC,T/T,T/T,T/T
rs761220536Family 1 + 2 + 7 + 134.08 × 10−6G/GTGACGGT,G/G,G/G,G/G
rs766751467Family 6 + 97.30 × 10−5GCA/G,GCA/GCA,GCA/.,GCA/GCA
rs766833662Family 1 + 2 + 7 + 130.00CTGGTGCG/C,CTGGTGCG/CTGGTGCG,CTGGTGCG/CTGGTGCG,CTGGTGCG/CTGGTGCG
rs776466780Family 8 + 121.51 × 10−4G/GTA,G/G,G/G,G/G
11TYK2Protein Structural International LocusNoneFamily 5 + 62.20 × 10−4TA/T,TA/TA,TA/TA,TA/TAYesYes
12TDGSplice Donor Variantrs760400700Family 2 + 103.63 × 10−2T/G,T/T,T/T,T/TYesNo
rs764159587Family 2 + 104.11 × 10−2G/GA,G/G,G/.,G/G
13RNASEH2BFrameshift Variantrs200320729Family 2 + 53.95 × 10−2T/TA,T/T,T/T,T/TNoYes
15TYRO3Splice Donor Variantrs746533465Family 5 + 137.81 × 10−5G/GAGAGTTTGGTTCAGTGCGGGAGGCCCAGC,G/G,G/G,G/GYesYes
rs750893216Family 1 + 61.37 × 10−3G/GGAGA,G/G,G/G,G/G
rs757748573Family 5 + 7 + 134.82 × 10−4G/C,G/G,G/G,G/G
17LGALS9CFrameshift Variantrs376412531Family 8 + 125.89 × 10−2GGA/G,GGA/GGA,GGA/GGA,GGA/GGANoYes
18LAMA3Splice Donor VariantNoneFamily 1 + 7NAT/G,T/T,T/TNoYes
19CNN2Stop GainedNoneFamily 3 + 8NAC/A,C/C,C/C,C/CNoYes
19COLGALT1Stop GainedNoneFamily 5 + 11NAG/A,G/G,G/G,G/GNoYes
22BCRFrameshift Variantrs372013175Family 5 + 122.79 × 10−3T/TCCGG,T/T,T/T,T/T,T/TNoYes
Chr: chromosome; AF: allele frequency; NA: not applicable; T1DM: type 1 diabetes mellitus. Allele frequency was obtained from the global population using gnomAD.broadinstitute.org. Family genotype is noted as follows, in chronological order: affected child, mother, father, sibling. Genotype of the affected child is in bold.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Mousa, M.; Albarguthi, S.; Albreiki, M.; Farooq, Z.; Sajid, S.; El Hajj Chehadeh, S.; ElBait, G.D.; Tay, G.; Deeb, A.A.; Alsafar, H. Whole-Exome Sequencing in Family Trios Reveals De Novo Mutations Associated with Type 1 Diabetes Mellitus. Biology 2023, 12, 413. https://doi.org/10.3390/biology12030413

AMA Style

Mousa M, Albarguthi S, Albreiki M, Farooq Z, Sajid S, El Hajj Chehadeh S, ElBait GD, Tay G, Deeb AA, Alsafar H. Whole-Exome Sequencing in Family Trios Reveals De Novo Mutations Associated with Type 1 Diabetes Mellitus. Biology. 2023; 12(3):413. https://doi.org/10.3390/biology12030413

Chicago/Turabian Style

Mousa, Mira, Sara Albarguthi, Mohammed Albreiki, Zenab Farooq, Sameeha Sajid, Sarah El Hajj Chehadeh, Gihan Daw ElBait, Guan Tay, Asma Al Deeb, and Habiba Alsafar. 2023. "Whole-Exome Sequencing in Family Trios Reveals De Novo Mutations Associated with Type 1 Diabetes Mellitus" Biology 12, no. 3: 413. https://doi.org/10.3390/biology12030413

APA Style

Mousa, M., Albarguthi, S., Albreiki, M., Farooq, Z., Sajid, S., El Hajj Chehadeh, S., ElBait, G. D., Tay, G., Deeb, A. A., & Alsafar, H. (2023). Whole-Exome Sequencing in Family Trios Reveals De Novo Mutations Associated with Type 1 Diabetes Mellitus. Biology, 12(3), 413. https://doi.org/10.3390/biology12030413

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop