Next generation sequencing reveals the association of DRB3*02:02 with type 1 diabetes

Fifty years of HLA-associated type 1 diabetes risk: history, current knowledge, and future directions

More than 50 years have elapsed since the association of human leukocyte antigens (HLA) with type 1 diabetes (T1D) was first reported. Since then, methods for identification of HLA have progressed from cell based to DNA based, and the number of recognized HLA variants has grown from a few to tens of thousands. Current genotyping methodology allows for exact identification of all HLA-encoding genes in an individual's genome, with statistical analysis methods evolving to digest the enormous amount of data that can be produced at an astonishing rate. The HLA region of the genome has been repeatedly shown to be the most important genetic risk factor for T1D, and the original reported associations have been replicated, refined, and expanded. Even with the remarkable progress through 50 years and over 5,000 reports, a comprehensive understanding of all effects of HLA on T1D remains elusive. This report represents a summary of the field as it evolved and as it stands now, enumerating many past and present challenges, and suggests possible paradigm shifts for moving forward with future studies in hopes of finally understanding all the ways in which HLA influences the pathophysiology of T1D.

Shared and distinct genetics of pure type 1 diabetes and type 1 diabetes with celiac disease, homology in their auto-antigens and immune dysregulation states: a study from North India

AIM

This study was undertaken to explicate the shared and distinctive genetic susceptibility and immune dysfunction in patients with T1D alone and T1D with CD (T1D + CD).

METHODS

A total of 100 T1D, 50 T1D + CD and 150 healthy controls were recruited. HLA-DRB1/DQB1 alleles were determined by PCR-sequence-specific primer method, SNP genotyping for CTLA-4 and PTPN22 was done by simple probe-based SNP-array and genotyping for INS-23 Hph1 A/T was done by RFLP. Autoantibodies and cytokine estimation was done by ELISA. Immune-regulation was analysed by flow-cytometry. Clustering of autoantigen epitopes was done by epitope cluster analytical tool.

RESULTS

Both T1D alone and T1D + CD had a shared association of DRB1*03:01, DRB1*04, DRB3*01:07/15 and DQB1*02. DRB3*01:07/15 confers the highest risk for T1D with relative risk of 11.32 (5.74-22.31). Non-HLA gene polymorphisms PTPN22 and INS could discriminate between T1D and T1D + CD. T1D + CD have significantly higher titers of autoantibodies, expression of costimulatory molecules on CD4 and CD8 cells, and cytokine IL-17A and TGF-β1 levels compared to T1D patients. Epitopes from immunodominant regions of autoantigens of T1D and CD clustered together with 40% homology.

CONCLUSION

Same HLA genes provide susceptibility for both T1D and CD. Non-HLA genes CTLA4, PTPN22 and INS provide further susceptibility while different polymorphisms in PTPN22 and INS can discriminate between T1D and T1D + CD. Epitope homology between autoantigens of two diseases further encourages the two diseases to occur together. The T1D + CD being more common in females along with co-existence of thyroid autoimmunity, and have more immune dysregulated state than T1D alone.

Evaluation of HLA typing content of next-generation sequencing datasets from family trios and individuals of arab ethnicity

Introduction: HLA typing is a critical tool in both clinical and research applications at the individual and population levels. Benchmarking studies have indicated HLA-HD as the preferred tool for accurate and comprehensive HLA allele calling. The advent of next-generation sequencing (NGS) has revolutionized genetic analysis by providing high-throughput sequencing data. This study aims to evaluate, using the HLA-HD tool, the HLA typing content of whole exome, whole genome, and HLA-targeted panel sequence data from the consanguineous population of Arab ethnicity, which has been underrepresented in prior benchmarking studies. Methods: We utilized sequence data from family trios and individuals, sequenced on one or more of the whole exome, whole genome, and HLA-targeted panel sequencing technologies. The performance and resolution across various HLA genes were evaluated. We incorporated a comparative quality control analysis, assessing the results obtained from HLA-HD by comparing them with those from the HLA-Twin tool to authenticate the accuracy of the findings. Results: Our analysis found that alleles across 29 HLA loci can be successfully and consistently typed from NGS datasets. Clinical-grade whole exome sequencing datasets achieved the highest consistency rate at three-field resolution, followed by targeted HLA panel, research-grade whole exome, and whole genome datasets. Discussion: The study catalogues HLA typing consistency across NGS datasets for a large array of HLA genes and highlights assessments regarding the feasibility of utilizing available NGS datasets in HLA allele studies. These findings underscore the reliability of HLA-HD for HLA typing in underrepresented populations and demonstrate the utility of various NGS technologies in achieving accurate HLA allele calling.

HLA Class II (DR, DQ, DP) Genes Were Separately Associated With the Progression From Seroconversion to Onset of Type 1 Diabetes Among Participants in Two Diabetes Prevention Trials (DPT-1 and TN07)

OBJECTIVE

To explore associations of HLA class II genes (HLAII) with the progression of islet autoimmunity from asymptomatic to symptomatic type 1 diabetes (T1D).

RESEARCH DESIGN AND METHODS

Next-generation targeted sequencing was used to genotype eight HLAII genes (DQA1, DQB1, DRB1, DRB3, DRB4, DRB5, DPA1, DPB1) in 1,216 participants from the Diabetes Prevention Trial-1 and Randomized Diabetes Prevention Trial with Oral Insulin sponsored by TrialNet. By the linkage disequilibrium, DQA1 and DQB1 are haplotyped to form DQ haplotypes; DP and DR haplotypes are similarly constructed. Together with available clinical covariables, we applied the Cox regression model to assess HLAII immunogenic associations with the disease progression.

RESULTS

First, the current investigation updated the previously reported genetic associations of DQA1*03:01-DQB1*03:02 (hazard ratio [HR] = 1.25, P = 3.50*10-3) and DQA1*03:03-DQB1*03:01 (HR = 0.56, P = 1.16*10-3), and also uncovered a risk association with DQA1*05:01-DQB1*02:01 (HR = 1.19, P = 0.041). Second, after adjusting for DQ, DPA1*02:01-DPB1*11:01 and DPA1*01:03-DPB1*03:01 were found to have opposite associations with progression (HR = 1.98 and 0.70, P = 0.021 and 6.16*10-3, respectively). Third, DRB1*03:01-DRB3*01:01 and DRB1*03:01-DRB3*02:02, sharing the DRB1*03:01, had opposite associations (HR = 0.73 and 1.44, P = 0.04 and 0.019, respectively), indicating a role of DRB3. Meanwhile, DRB1*12:01-DRB3*02:02 and DRB1*01:03 alone were found to associate with progression (HR = 2.6 and 2.32, P = 0.018 and 0.039, respectively). Fourth, through enumerating all heterodimers, it was found that both DQ and DP could exhibit associations with disease progression.

CONCLUSIONS

These results suggest that HLAII polymorphisms influence progression from islet autoimmunity to T1D among at-risk subjects with islet autoantibodies.

Genetic and Epigenetic Aspects of Type 1 Diabetes Mellitus: Modern View on the Problem

Omics technologies accumulated an enormous amount of data that advanced knowledge about the molecular pathogenesis of type 1 diabetes mellitus and identified a number of fundamental problems focused on the transition to personalized diabetology in the future. Among them, the most significant are the following: (1) clinical and genetic heterogeneity of type 1 diabetes mellitus; (2) the prognostic significance of DNA markers beyond the HLA genes; (3) assessment of the contribution of a large number of DNA markers to the polygenic risk of disease progress; (4) the existence of ethnic population differences in the distribution of frequencies of risk alleles and genotypes; (5) the infancy of epigenetic research into type 1 diabetes mellitus. Disclosure of these issues is one of the priorities of fundamental diabetology and practical healthcare. The purpose of this review is the systemization of the results of modern molecular genetic, transcriptomic, and epigenetic investigations of type 1 diabetes mellitus in general, as well as its individual forms. The paper summarizes data on the role of risk HLA haplotypes and a number of other candidate genes and loci, identified through genome-wide association studies, in the development of this disease and in alterations in T cell signaling. In addition, this review assesses the contribution of differential DNA methylation and the role of microRNAs in the formation of the molecular pathogenesis of type 1 diabetes mellitus, as well as discusses the most currently central trends in the context of early diagnosis of type 1 diabetes mellitus.

Identification of gene mutations associated with type 1 diabetes by next-generation sequencing in affected Palestinian families

Introduction: Diabetes Mellitus is a group of metabolic disorders characterized by hyperglycemia secondary to insulin resistance or deficiency. It is considered a major health problem worldwide. T1DM is a result of a combination of genetics, epigenetics, and environmental factors. Several genes have been associated with T1DM, including HLA, INS, CTLA4, and PTPN22. However, none of these findings have been based on linkage analysis because it is rare to find families with several diabetic individuals. Two Palestinian families with several afflicted members with variations in the mode of inheritance were identified and selected for this study. This study aimed to identify the putative causative gene(s) responsible for T1DM development in these families to improve our understanding of the molecular genetics of the disease. Methods: One afflicted member from each family was selected for Whole-Exome Sequencing. Data were mapped to the reference of the human genome, and the resulting VCF file data were filtered. The variants with the highest phenotype correlation score were checked by Sanger sequencing for all family members. The confirmed variants were analyzed in silico by bioinformatics tools. Results: In one family, the IGF1R p.V579F variant, which follows autosomal dominant inheritance, was confirmed and segregated in the family. In another family, the NEUROD1 p.P197H variant, which follows autosomal recessive inheritance, was positively confirmed and segregated. Conclusion: IGF1R p.V579F and NEUROD1 p.P197H variants were associated with T1DM development in the two inflicted families. Further analysis and functional assays will be performed, including the generation of mutant model cell systems, to unravel their specific molecular mechanism in the disease development.

A vascularized 3D model of the human pancreatic islet forex vivostudy of immune cell-islet interaction

Insulin is an essential regulator of blood glucose homeostasis that is produced exclusively byβcells within the pancreatic islets of healthy individuals. In those affected by diabetes, immune inflammation, damage, and destruction of isletβcells leads to insulin deficiency and hyperglycemia. Current efforts to understand the mechanisms underlyingβcell damage in diabetes rely onin vitro-cultured cadaveric islets. However, isolation of these islets involves removal of crucial matrix and vasculature that supports islets in the intact pancreas. Unsurprisingly, these islets demonstrate reduced functionality over time in standard culture conditions, thereby limiting their value for understanding native islet biology. Leveraging a novel, vascularized micro-organ (VMO) approach, we have recapitulated elements of the native pancreas by incorporating isolated human islets within a three-dimensional matrix nourished by living, perfusable blood vessels. Importantly, these islets show long-term viability and maintain robust glucose-stimulated insulin responses. Furthermore, vessel-mediated delivery of immune cells to these tissues provides a model to assess islet-immune cell interactions and subsequent islet killing-key steps in type 1 diabetes pathogenesis. Together, these results establish the islet-VMO as a novel,ex vivoplatform for studying human islet biology in both health and disease.

Update on the ZNT8 epitope and its role in the pathogenesis of type 1 diabetes

Type 1 diabetes (T1D) is an organ-specific chronic autoimmune disease mediated by autoreactive T cells. ZnT8 is a pancreatic islet-specific zinc transporter that is mainly located in β cells. It not only participates in the synthesis, storage and secretion of insulin but also maintains the structural integrity of insulin. ZnT8 is the main autoantigen recognized by autoreactive CD8+ T cells in children and adults with T1D. This article summarizes the latest research results on the T lymphocyte epitope and B lymphocyte epitope of ZnT8 in the current literature. The structure and expression of ZnT8, the role of ZnT8 in insulin synthesis and its role in autoimmunity are reviewed. ZnT8 is the primary autoantigen of T1D and is specifically expressed in pancreatic islets. Thus, it is one of biomarkers for the diagnosis of T1D. It has broad prospects for further research on immunomodulators for the treatment of T1D.

Investigation of molecular interaction and conformational stability of disease concomitant to HLA-DRβ3

Human leucocyte antigen DRβ3 is associated with specific autoimmune thyroid disease and plays a vital role in the progression of Grave's disease. The available crystallographic structure of the HLA DRA, DRβ3*0101, was selected and used to generate mutation at position 57 from valine amino acid to Aspartic acid (D), Glutamic acid (E), Alanine (A), and Serine (S) amino acids by computational modeling approach. Mutant models were minimized, and stable conformation was chosen based on the lowest root mean square deviation value. Molecular docking assessed the best binding affinity of ligands C1, C2, C3, and C4 with wild-type and mutant HLA-DRβ3 models. Molecular dynamics simulation studies were executed to evaluate the stability of selected hits with wild-type and mutant dock complexes. The C3 has shown good binding affinity with wild-type and selected mutants; V57A, V57E, and V57D. Structural and molecular dynamics insights reveal the differences between wild-type and mutant-type HLA-DRβ3, which could help design novel antagonist molecules against autoimmune thyroid disorder.Communicated by Ramaswamy H. Sarma.

The β Cell in Diabetes: Integrating Biomarkers With Functional Measures

The pathogenesis of hyperglycemia observed in most forms of diabetes is intimately tied to the islet β cell. Impairments in propeptide processing and secretory function, along with the loss of these vital cells, is demonstrable not only in those in whom the diagnosis is established but typically also in individuals who are at increased risk of developing the disease. Biomarkers are used to inform on the state of a biological process, pathological condition, or response to an intervention and are increasingly being used for predicting, diagnosing, and prognosticating disease. They are also proving to be of use in the different forms of diabetes in both research and clinical settings. This review focuses on the β cell, addressing the potential utility of genetic markers, circulating molecules, immune cell phenotyping, and imaging approaches as biomarkers of cellular function and loss of this critical cell. Further, we consider how these biomarkers complement the more long-established, dynamic, and often complex measurements of β-cell secretory function that themselves could be considered biomarkers.

Identification of the novel HLA-DRB3*02:02:19 allele

The HLA-DRB3*02:02:19 allele differs from DRB3*02:02:01:02 by a single nucleotide change in exon 2.

High-resolution genotyping indicates that children with type 1 diabetes and celiac disease share three HLA class II loci in DRB3, DRB4 and DRB5 genes

Type 1 diabetes (T1D) and celiac disease (CD) share common genetic loci, mainly within the human leukocyte antigen (HLA) class II complex. Extended genotyping of HLA class II alleles and their potential risk for developing both diseases remains to be studied. The present study compared extended HLA-class II gene polymorphisms in children with T1D, CD, and a subgroup diagnosed with both diseases (T1D w/CD). Next-generation targeted sequencing (NGTS) of HLA-DRB3, DRB4, DRB5, DRB1, DQA1, DQB1, DPA1, and DPB1 alleles from DNA collected from 68 T1D, 219 CD, and seven T1D w/CD patients were compared with 636 HLA-genotyped Swedish children from the general population selected as controls. In comparison to controls, the DRB4*01:03:01 allele occurred more frequently in T1D w/CD (odds ratio (OR) = 7.84; 95% confidence interval (95% CI) = (2.24, 34.5), P = 0.0002) and T1D (OR = 3.86; 95% CI, (2.69, 5.55), P = 1.07 × 10^-14 ), respectively. The DRB3*01:01:02 allele occurred more frequently in CD as compared to controls (OR = 7.87; 95% CI, (6.17, 10.03), P = 4.24 × 10^-71 ), but less frequently in T1D (OR = 2.59; 95% CI, (1.76, 3.81), P = 7.29 × 10^-07 ) and T1D w/CD (OR = 0.87; 95% CI, (0.09, 3.96), P ≤ 0.999). The frequency of the DRB4*01:03:01-DRB1*04:01:01-DQA1*03:01:01-DQB1*03:02:01 (DR4-DQ8) haplotype was higher in T1D w/CD (OR = 12.88; 95% CI (4.35, 38.14) P = 3.75 × 10^-9 ), and moderately higher in T1D (OR = 2.13; 95% CI (1.18, 3.83) P = 0.01) compared with controls, but comparable in CD (OR = 1.45; 95% CI (0.94, 2.21), P = 0.08) and controls. Children with T1D and CD are associated with DRB4*01:03:01, DRB3*01:01:02, and DRB3*02:02:01 of which DRB4*01:03:01 confers the strongest risk allele for developing T1D w/CD.

Genome wide association studies for treatment-related adverse effects of pediatric acute lymphoblastic leukemia

Acute lymphoblastic leukemia (ALL) is the most common pediatric hematological malignancy; notwithstanding the success of ALL therapy, severe adverse drugs effects represent a serious issue in pediatric oncology, because they could be both an additional life threatening condition for ALL patients per se and a reason to therapy delay or discontinuation with important fallouts on final outcome. Cancer treatment-related toxicities have generated a significant need of finding predictive pharmacogenomic markers for the a priori identification of at risk patients. In the era of precision medicine, high throughput genomic screening such as genome wide association studies (GWAS) might provide useful markers to tailor therapy intensity on patients' genetic profile. Furthermore, these findings could be useful in basic research for better understanding the mechanistic and regulatory pathways of the biological functions associated with ALL treatment toxicities. The purpose of this review is to give an overview of high throughput genomic screening of the last 10 years that had investigated the landscape of ALL treatment-associated toxicities. This article is categorized under: Cancer > Genetics/Genomics/Epigenetics.

Genome-wide analyses disclose the distinctive HLA architecture and the pharmacogenetic landscape of the Somali population

African populations are underrepresented in medical genomics studies. For the Somali population, there is virtually no information on genomic markers with significance to precision medicine. Here, we analyzed nearly 900,000 genomic markers in samples collected from 95 unrelated individuals in the North Eastern Somalia. ADMIXTURE program for estimation of individual ancestries revealed a homogenous Somali population. Principal component analysis with PLINK software showed approximately 60% East African and 40% West Eurasian genes in the Somali population, with a close relation to the Cushitic and Semitic speaking Ethiopian populations. We report the unique features of human leukocyte antigens (HLA) in the Somali population, which seem to differentiate from all other neighboring regions compared. Current study identified high prevalence of the diabetes type 1 (T1D) predisposing HLA DR-DQ haplotypes in Somalia. This finding may explain the increased T1D risk observed among Somali children. In addition, ethnic Somalis were found to host the highest frequencies observed thus far for several pharmacogenetic variants, including UGT1A4*2. In conclusion, we report that the Somali population displays genetic traits of significance to health and disease. The Somali dataset is publicly available and will add more information to the few genomic datasets available for African populations.

Molecular mimicry between Anoctamin 2 and Epstein-Barr virus nuclear antigen 1 associates with multiple sclerosis risk

Multiple sclerosis (MS) is a chronic inflammatory, likely autoimmune disease of the central nervous system with a combination of genetic and environmental risk factors, among which Epstein-Barr virus (EBV) infection is a strong suspect. We have previously identified increased autoantibody levels toward the chloride-channel protein Anoctamin 2 (ANO2) in MS. Here, IgG antibody reactivity toward ANO2 and EBV nuclear antigen 1 (EBNA1) was measured using bead-based multiplex serology in plasma samples from 8,746 MS cases and 7,228 controls. We detected increased anti-ANO2 antibody levels in MS (P = 3.5 × 10^-36) with 14.6% of cases and 7.8% of controls being ANO2 seropositive (odds ratio [OR] = 1.6; 95% confidence intervals [95%CI]: 1.5 to 1.8). The MS risk increase in ANO2-seropositive individuals was dramatic when also exposed to 3 known risk factors for MS: HLA-DRB1*15:01 carriage, absence of HLA-A*02:01, and high anti-EBNA1 antibody levels (OR = 24.9; 95%CI: 17.9 to 34.8). Reciprocal blocking experiments with ANO2 and EBNA1 peptides demonstrated antibody cross-reactivity, mapping to ANO2 [aa 140 to 149] and EBNA1 [aa 431 to 440]. HLA gene region was associated with anti-ANO2 antibody levels and HLA-DRB1*04:01 haplotype was negatively associated with ANO2 seropositivity (OR = 0.6; 95%CI: 0.5 to 0.7). Anti-ANO2 antibody levels were not increased in patients from 3 other inflammatory disease cohorts. The HLA influence and the fact that specific IgG production usually needs T cell help provides indirect evidence for a T cell ANO2 autoreactivity in MS. We propose a hypothesis where immune reactivity toward EBNA1 through molecular mimicry with ANO2 contributes to the etiopathogenesis of MS.

High resolution HLA analysis reveals independent class I haplotypes and amino-acid motifs protective for multiple sclerosis

We investigated association between HLA class I and class II alleles and haplotypes, and KIR loci and their HLA class I ligands, with multiple sclerosis (MS) in 412 European American MS patients and 419 ethnically matched controls, using next-generation sequencing. The DRB1*15:01~DQB1*06:02 haplotype was highly predisposing (odds ratio (OR) = 3.98; 95% confidence interval (CI) = 3-5.31; p-value (p) = 2.22E-16), as was DRB1*03:01~DQB1*02:01 (OR = 1.63; CI = 1.19-2.24; p = 1.41E-03). Hardy-Weinberg (HW) analysis in MS patients revealed a significant DRB1*03:01~DQB1*02:01 homozyote excess (15 observed; 8.6 expected; p = 0.016). The OR for this genotype (5.27; CI = 1.47-28.52; p = 0.0036) suggests a recessive MS risk model. Controls displayed no HW deviations. The C*03:04~B*40:01 haplotype (OR = 0.27; CI = 0.14-0.51; p = 6.76E-06) was highly protective for MS, especially in haplotypes with A*02:01 (OR = 0.15; CI = 0.04-0.45; p = 6.51E-05). By itself, A*02:01 is moderately protective, (OR = 0.69; CI = 0.54-0.87; p = 1.46E-03), and haplotypes of A*02:01 with the HLA-B Thr80 Bw4 variant (Bw4T) more so (OR = 0.53; CI = 0.35-0.78; p = 7.55E-04). Protective associations with the Bw4 KIR ligand resulted from linkage disequilibrium (LD) with DRB1*15:01, but the Bw4T variant was protective (OR = 0.64; CI = 0.49-0.82; p = 3.37-04) independent of LD with DRB1*15:01. The Bw4I variant was not associated with MS. Overall, we find specific class I HLA polymorphisms to be protective for MS, independent of the strong predisposition conferred by DRB1*15:01.

Association of protein tyrosine phosphatase non-receptor type 22 gene functional variant C1858T, HLA-DQ/DR genotypes and autoantibodies with susceptibility to type-1 diabetes mellitus in Kuwaiti Arabs

The incidence of type-1 Diabetes Mellitus (T1DM) has increased steadily in Kuwait during recent years and it is now considered amongst the high-incidence countries. An interaction between susceptibility genes, immune system mediators and environmental factors predispose susceptible individuals to T1DM. We have determined the prevalence of protein tyrosine phosphatase non-receptor type 22 (PTPN22) gene functional variant (C1858T; R620W, rs2476601), HLA-DQ and DR alleles and three autoantibodies in Kuwaiti children with T1DM to evaluate their impact on genetic predisposition of the disease. This study included 253 Kuwaiti children with T1DM and 214 ethnically matched controls. The genotypes of PTPN22 gene functional variant C1858T (R620W; rs2476601) were detected by PCR-RFLP method and confirmed by DNA sequencing. HLA-DQ and DR alleles were determined by sequence-specific PCR. Three autoantibodies were detected in the T1DM patients using radio-immunoassays. A significant association was detected between the variant genotype of the PTPN22 gene (C1858T, rs2476601) and T1DM in Kuwaiti Arabs. HLA-DQ2 and DQ8 alleles showed a strong association with T1DM. In T1DM patients which carried the variant TT-genotype of the PTPN22 gene, 93% had at least one DQ2 allele and 60% carried either a DQ2 or a DQ8 allele. Amongst the DR alleles, the DR3-DRB5, DR3-3, DR3-4 and DR4-4 showed a strong association with T1DM. Majority of T1DM patients who carried homozygous variant (TT) genotype of the PTPN22 gene had either DR3-DRB5 or DRB3-DRB4 genotypes. In T1DM patients who co-inherited the high risk HLA DQ, DR alleles with the variant genotype of PTPN22 gene, the majority were positive for three autoantibodies. Our data demonstrate that the variant T-allele of the PTPN22 gene along with HLA-DQ2 and DQ8 alleles constitute significant determinants of genetic predisposition of T1DM in Kuwaiti children.

Genetics of type 1 diabetes

Type 1 diabetes (T1D) results from immune-mediated loss of pancreatic beta cells leading to insulin deficiency. It is the most common form of diabetes in children, and its incidence is on the rise. This article reviews the current knowledge on the genetics of T1D. In particular, we discuss the influence of HLA and non-HLA genes on T1D risk and disease progression through the preclinical stages of the disease, and the development of genetic scores that can be applied to disease prediction. Racial/ethnic differences, challenges and future directions in the genetics of T1D are also discussed.

Early prediction of autoimmune (type 1) diabetes

Underlying type 1 diabetes is a genetic aetiology dominated by the influence of specific HLA haplotypes involving primarily the class II DR-DQ region. In genetically predisposed children with the DR4-DQ8 haplotype, exogenous factors, yet to be identified, are thought to trigger an autoimmune reaction against insulin, signalled by insulin autoantibodies as the first autoantibody to appear. In children with the DR3-DQ2 haplotype, the triggering reaction is primarily against GAD signalled by GAD autoantibodies (GADA) as the first-appearing autoantibody. The incidence rate of insulin autoantibodies as the first-appearing autoantibody peaks during the first years of life and declines thereafter. The incidence rate of GADA as the first-appearing autoantibody peaks later but does not decline. The first autoantibody may variably be followed, in an apparently non-HLA-associated pathogenesis, by a second, third or fourth autoantibody. Although not all persons with a single type of autoantibody progress to diabetes, the presence of multiple autoantibodies seems invariably to be followed by loss of functional beta cell mass and eventually by dysglycaemia and symptoms. Infiltration of mononuclear cells in and around the islets appears to be a late phenomenon appearing in the multiple-autoantibody-positive with dysglycaemia. As our understanding of the aetiology and pathogenesis of type 1 diabetes advances, the improved capability for early prediction should guide new strategies for the prevention of type 1 diabetes.

Impact of Human Leukocyte Antigen Loci and Haplotypes on Intestinal Acute Graft-versus-host Disease after Human Leukocyte Antigen-matched Sibling Peripheral Blood Stem Cell Transplantation

Background:

Acute graft-versus-host disease (aGVHD) is a common and severe complication of allogeneic hematopoietic stem cell transplantation (allo-HSCT). Some studies have found that the presence of certain specific human leukocyte antigen (HLA) loci could affect the occurrence of aGVHD. Meanwhile, the impact of HLA haplotypes on aGVHD has been rarely studied. This study aimed to investigate the effects of HLA loci and haplotypes on intestinal aGVHD.

Methods:

Totally, 345 consecutive patients undergoing first HLA-matched sibling peripheral blood stem cell transplantation (PBSCT) from February 2004 to June 2013 at Institute of Hematology and Blood Diseases Hospital, Chinese Academy of Medical Sciences and Peking Union Medical College, were enrolled in this study. HLA loci and haplotypes of recipients with frequency over 5% were searched and their effects on intestinal aGVHD were investigated. Other important factors including donor age, recipient age, donor-recipient sex combinations, and conditioning regimens were also evaluated using logistic regression. Pure upper gastrointestinal tract aGVHD without diarrhea was excluded because the histological proof was unavailable. The follow-up end-point was 6 months after HSCT.

Results:

The cumulative incidence of intestinal aGVHD was 19.4%, with 18.0% of the patients classified as classic aGVHD and 1.4% as persistent, recurrent, or late aGVHD. Multivariate analysis showed that HLA-A31 locus (odds ratio [OR] 2.893, 95% confidence interval [CI] [1.054, 7.935], P = 0.039), HLA B40-DR15 (OR 3.133, 95% CI [1.250, 7.857], P = 0.015), and HLA B46-DR9 haplotypes (OR 2.580, 95% CI [1.070, 6.220], P = 0.035), female donor for male recipient (OR 2.434, 95% CI [1.319, 4.493], P = 0.004) were risk factors for intestinal aGVHD.

Conclusion:

The presence of certain HLA loci and haplotypes may influence the occurrence of intestinal aGVHD in PBSCT with HLA-identical sibling donors.

Different DRB1*03:01-DQB1*02:01 haplotypes confer different risk for celiac disease

Celiac disease is associated with the HLA-DR3-DQA1*05:01-DQB1*02:01 and DR4-DQA1*03:01-DQB1*03:02 haplotypes. In addition, there are currently over 40 non-HLA loci associated with celiac disease. This study extends previous analyses on different HLA haplotypes in celiac disease using next generation targeted sequencing. Included were 143 patients with celiac disease and 135 non-celiac disease controls investigated at median 9.8 years (1.4-18.3 years). PCR-based amplification of HLA and sequencing with Illumina MiSeq technology were used for extended sequencing of the HLA class II haplotypes HLA-DRB1, DRB3, DRB4, DRB5, DQA1 and DQB1, respectively. Odds ratios were computed marginally for every allele and haplotype as the ratio of allelic frequency in patients and controls as ratio of exposure rates (RR), when comparing a null reference with equal exposure rates in cases and controls. Among the extended HLA haplotypes, the strongest risk haplotype for celiac disease was shown for DRB3*01:01:02 in linkage with DQA1*05:01-DQB1*02:01 (RR = 6.34; P-value < .0001). In a subpopulation analysis, DRB3*01:01:02-DQA1*05:01-DQB1*02:01 remained the most significant in patients with Scandinavian ethnicity (RR = 4.63; P < .0001) whereas DRB1*07:01:01-DRB4*01:03:01-DQA1*02:01-DQB1*02:02:01 presented the highest risk of celiac disease among non-Scandinavians (RR = 7.94; P = .011). The data also revealed 2 distinct celiac disease risk DR3-DQA1*05:01-DQB*02:01 haplotypes distinguished by either the DRB3*01:01:02 or DRB3*02:02:01 alleles, indicating that different DRB1*03:01-DQB1*02:01 haplotypes confer different risk for celiac disease. The associated risk of celiac disease for DR3-DRB3*01:01:02-DQA1*05:01-DQB1*02:01 is predominant among patients of Scandinavian ethnicity.

Type 1 diabetes mellitus

Type 1 diabetes mellitus (T1DM), also known as autoimmune diabetes, is a chronic disease characterized by insulin deficiency due to pancreatic β-cell loss and leads to hyperglycaemia. Although the age of symptomatic onset is usually during childhood or adolescence, symptoms can sometimes develop much later. Although the aetiology of T1DM is not completely understood, the pathogenesis of the disease is thought to involve T cell-mediated destruction of β-cells. Islet-targeting autoantibodies that target insulin, 65 kDa glutamic acid decarboxylase, insulinoma-associated protein 2 and zinc transporter 8 - all of which are proteins associated with secretory granules in β-cells - are biomarkers of T1DM-associated autoimmunity that are found months to years before symptom onset, and can be used to identify and study individuals who are at risk of developing T1DM. The type of autoantibody that appears first depends on the environmental trigger and on genetic factors. The pathogenesis of T1DM can be divided into three stages depending on the absence or presence of hyperglycaemia and hyperglycaemia-associated symptoms (such as polyuria and thirst). A cure is not available, and patients depend on lifelong insulin injections; novel approaches to insulin treatment, such as insulin pumps, continuous glucose monitoring and hybrid closed-loop systems, are in development. Although intensive glycaemic control has reduced the incidence of microvascular and macrovascular complications, the majority of patients with T1DM are still developing these complications. Major research efforts are needed to achieve early diagnosis, prevent β-cell loss and develop better treatment options to improve the quality of life and prognosis of those affected.

Predominance of DR3 in Somali children with type 1 diabetes in the twin cities, Minnesota

BACKGROUND

Minnesota is home to the largest Somali population in USA, and pediatric diabetes teams are seeing increasing numbers of Somali children with diabetes.

OBJECTIVE

To assess the immune basis of diabetes in Somali children in the Twin Cities, Minnesota.

METHODS

A total of 31 Somali children ≤19 yr were treated for type 1 diabetes (T1D) at the University of Minnesota Masonic Children's Hospital and Children's Hospitals and Clinics of Minnesota underwent analysis of human leukocyte antigen (HLA) alleles (n = 30) and diabetes autoantibodies [glutamic acid decarboxylase (GAD65), islet antigen 2 (IA-2), zinc transporter 8 (ZnT8); n = 31]. HLA alleles were analyzed in 49 Somalis without diabetes (controls). Anti-transglutaminase autoantibodies (TGA) for celiac disease were also measured.

RESULTS

In Somali children with T1D aged 13.5 ± 5 yr (35% female, disease duration 6.5 ± 3.6 yr), the most common HLA allele was DRB1*03:01 (93%, compared with 45% of Somali controls), followed by DRB1*13:02 (27%). There was a relatively low frequency of DR4 (13%). Controls showed a similar pattern. All 31 participants were positive for at least one diabetes autoantibody. Insulin antibodies were positive in 84% (all were on insulin). Excluding insulin antibodies, 23 (74%) subjects tested positive for at least one other diabetes autoantibody; 32% had 1 autoantibody, 32% had 2 autoantibodies, and 10% had 3 autoantibodies. GAD65 autoantibodies were found in 56% of subjects, IA-2 in 29%, and ZnT8 in 26%. Four (13%) were TGA positive.

CONCLUSION

The autoantibody and HLA profiles of Somali children with diabetes are consistent with autoimmune diabetes. Their HLA profile is unique with an exceptionally high prevalence of DRB1*03:01 allele and relative paucity of DR4 alleles compared with African Americans with T1D.

Integrated analysis of different microarray studies to identify candidate genes in type 1 diabetes

BACKGROUND

Type 1 diabetes (T1D), an autoimmune disease, occurs most commonly in children. Identifying altered gene expression in peripheral blood mononuclear cells (PBMCs) of T1D may lead to new strategies for preserving or improving β-ell function in patients with T1D.

METHODS

The Gene Expression Omnibus database was searched for microarray studies in PBMCs of T1D. Subsequently, gene expression datasets from multiple microarray studies were integrated to obtain differentially expressed genes (DEGs) between T1D and normal controls (NC). Gene function analysis was performed to determine the functions of the DEGs identified.

RESULTS

Four microarray studies were available for analysis, including 199 T1D samples and 74 NC samples. Analysis revealed 695 genes that were significantly differentially expressed in PBMCs from T1D compared with NC samples, with 450 upregulated and 245 downregulated. Signal transduction (gene ontology [GO]: 0007165; false discovery rate [FDR] = 1.54 × 10^-7 ) and protein binding (GO: 0005515; FDR = 2.93 × 10^-24 ) were significantly enriched for the GO categories of biological processes and molecular functions, respectively. The most significant pathway in the Kyoto Encyclopedia of Genes and Genomes analysis was arachidonic acid metabolism (FDR = 1.44 × 10^-3 ). Protein-protein interaction network analysis showed that the significant hub proteins contained immature colon carcinoma transcript 1 (ICT1; degree = 214; clustering coefficient [C] = 4.39 × 10^-5 ), zinc finger and BTB domain containing 16 (ZBTB16; degree = 112; C = 8.04 × 10^-4 ), and SERTA domain containing 1 (SERTAD1; degree = 38; C = 0.0014).

CONCLUSIONS

This integrated analysis will help develop improved therapies and interventions for T1D by identifying novel drug targets.

Genetic risk factors for type 1 diabetes

Type 1 diabetes is diagnosed at the end of a prodrome of β-cell autoimmunity. The disease is most likely triggered at an early age by autoantibodies primarily directed against insulin or glutamic acid decarboxylase, or both, but rarely against islet antigen-2. After the initial appearance of one of these autoantibody biomarkers, a second, third, or fourth autoantibody against either islet antigen-2 or the ZnT8 transporter might also appear. The larger the number of β-cell autoantibody types, the greater the risk of rapid progression to clinical onset of diabetes. This association does not necessarily mean that the β-cell autoantibodies are pathogenic, but rather that they represent reproducible biomarkers of the pathogenesis. The primary risk factor for β-cell autoimmunity is genetic, mainly occurring in individuals with either HLA-DR3-DQ2 or HLA-DR4-DQ8 haplotypes, or both, but a trigger from the environment is generally needed. The pathogenesis can be divided into three stages: 1, appearance of β-cell autoimmunity, normoglycaemia, and no symptoms; 2, β-cell autoimmunity, dysglycaemia, and no symptoms; and 3, β-cell autoimmunity, dysglycaemia, and symptoms of diabetes. The genetic association with each one of the three stages can differ. Type 1 diabetes could serve as a disease model for organ-specific autoimmune disorders such as coeliac disease, thyroiditis, and Addison's disease, which show similar early markers of a prolonged disease process before clinical diagnosis.

Deep sequencing of the MHC region in the Chinese population contributes to studies of complex disease

The human major histocompatibility complex (MHC) region has been shown to be associated with numerous diseases. However, it remains a challenge to pinpoint the causal variants for these associations because of the extreme complexity of the region. We thus sequenced the entire 5-Mb MHC region in 20,635 individuals of Han Chinese ancestry (10,689 controls and 9,946 patients with psoriasis) and constructed a Han-MHC database that includes both variants and HLA gene typing results of high accuracy. We further identified multiple independent new susceptibility loci in HLA-C, HLA-B, HLA-DPB1 and BTNL2 and an intergenic variant, rs118179173, associated with psoriasis and confirmed the well-established risk allele HLA-C*06:02. We anticipate that our Han-MHC reference panel built by deep sequencing of a large number of samples will serve as a useful tool for investigating the role of the MHC region in a variety of diseases and thus advance understanding of the pathogenesis of these disorders.

Next-Generation Sequencing Reveals That HLA-DRB3, -DRB4, and -DRB5 May Be Associated With Islet Autoantibodies and Risk for Childhood Type 1 Diabetes

The possible contribution of HLA-DRB3, -DRB4, and -DRB5 alleles to type 1 diabetes risk and to insulin autoantibody (IAA), GAD65 (GAD autoantibody [GADA]), IA-2 antigen (IA-2A), or ZnT8 against either of the three amino acid variants R, W, or Q at position 325 (ZnT8RA, ZnT8WA, and ZnT8QA, respectively) at clinical diagnosis is unclear. Next-generation sequencing (NGS) was used to determine all DRB alleles in consecutively diagnosed patients ages 1-18 years with islet autoantibody-positive type 1 diabetes (n = 970) and control subjects (n = 448). DRB3, DRB4, or DRB5 alleles were tested for an association with the risk of DRB1 for autoantibodies, type 1 diabetes, or both. The association between type 1 diabetes and DRB1*03:01:01 was affected by DRB3*01:01:02 and DRB3*02:02:01. These DRB3 alleles were associated positively with GADA but negatively with ZnT8WA, IA-2A, and IAA. The negative association between type 1 diabetes and DRB1*13:01:01 was affected by DRB3*01:01:02 to increase the risk and by DRB3*02:02:01 to maintain a negative association. DRB4*01:03:01 was strongly associated with type 1 diabetes (P = 10(-36)), yet its association was extensively affected by DRB1 alleles from protective (DRB1*04:03:01) to high (DRB1*04:01:01) risk, but its association with DRB1*04:05:01 decreased the risk. HLA-DRB3, -DRB4, and -DRB5 affect type 1 diabetes risk and islet autoantibodies. HLA typing with NGS should prove useful to select participants for prevention or intervention trials.

Population-dependent contribution of the major histocompatibility complex region to schizophrenia susceptibility

There is consistent data from European cohorts suggesting a genetic contribution from the major histocompatibility complex (MHC) to the pathogenesis of schizophrenia. However, the genomic complexity and ethnicity-specific diversity found in the MHC cause difficulties in identifying causal variants or genes, and there is a need for studies encompassing the entire MHC region in multiple ethnic populations. Here, we report on association signals in the MHC region, with schizophrenia in the Japanese population. We genotyped and imputed a total of 10,131 single nucleotide polymorphisms (SNPs), spanning the entire MHC interval. The analysis included 3302 participants (1518 schizophrenics and 1784 healthy controls) from the Japanese population. In this study, we present evidence for association at rs494620, located in the SLC44A4 gene. The association survived after correction for multiple testing (unadjusted P=7.78×10(-5), empirical P=0.0357). The imputation results detected the highest association at rs707937 in the MSH5-SAPCD1 gene (imputed P=8.40×10(-5)). In expression analysis using postmortem brains from schizophrenia and control samples, MSH5-SAPCD1 showed marginally significant expression differences in Brodmann's area 46 (P=0.044 by unpaired t test with Welch's correction, P=0.099 by Mann-Whitney U test). Our study further strengthens evidence for the involvement of the MHC in schizophrenia across populations, and provides insight into population-specific mechanisms for the MHC region in schizophrenia susceptibility.

DNA storage under high temperature conditions does not affect performance in human leukocyte antigen genotyping via next-generation sequencing (DNA integrity maintained in extreme conditions)

BACKGROUND

Stable dry-state storage of DNA is desirable to minimize required storage space and to reduce electrical and shipping costs. DNA purified from various commercially available dry-state stabilization matrices has been used successfully in downstream molecular applications (e.g., quantitative polymerase chain reaction [qPCR], microarray, and sequence-based genotyping). However, standard DNA storage conditions still include freezing of DNA eluted in aqueous buffers or nuclease-free water. Broad implementation of dry-state, long-term DNA storage requires enhancement of such dry-state DNA stabilization products to control for temperature fluctuations at specimen collection, transit, and storage. This study tested the integrity of genomic DNA subjected to long-term storage on GenTegra(™) DNA stabilization matrices (GenTegra LLC, Pleasanton, CA) at extreme conditions, as defined by a 4-year storage period at ambient temperature with an initial incubation for 7 months at 37°C, 56°C, or ambient temperature. Subsequently, purified DNA performance and integrity were measured by qPCR and next-generation sequencing (NGS)-based human leokocyte antigen (HLA) genotyping.

RESULTS

High molecular weight genomic DNA samples were recovered from the GenTegra product matrix and exhibited integrity comparable to a highly characterized commercial standard under assessment by qPCR. Samples were genotyped for classical HLA loci using next generation sequencing-based methodolgy on the Roche 454 GS Junior instrument. Amplification efficiency, sequence coverage, and sequence quality were all comparable with those produced from a cell line DNA sequenced as a control. No significant differences were observed in the mean, median, or mode quality scores between samples and controls (p≥0.4).

CONCLUSIONS

Next generation HLA genotyping was chosen to test the integrity of GenTegra-treated genomic DNA due to the requirment for long sequence reads to genotype the highly polymorphic classical HLA genes. Experimental results demonstrate the efficacy of the GenTegra product as a suitable genomic DNA preservation tool for collection and long-term biobanking of DNA at fluctuating and high temperatures.

HLA class II genotypes are not associated with age related macular degeneration in a case-control, population-based study

Multiple lines of evidence support an immunologic basis and genetic disposition for the development of age-related macular degeneration (AMD). Comprehensive human leukocyte antigens (HLA) class II typing at four loci (DRB1, DQA1, DQB1, and DPB1) was assessed using next generation sequencing methods and tested for association with age-related macular degeneration (AMD) in a case-control study of 456 AMD cases and 499 controls from the population-based Study of Osteoporotic Fractures (SOF) cohort. No statistically significant associations were identified for any of the class II loci and a previously identified association between DRB1*13:01 was not replicated in this dataset. These results reported here suggest that common HLA class II genetic variation does not contribute to AMD disease risk.

Type I diabetes mellitus: genetic factors and presumptive enteroviral etiology or protection

We review type 1 diabetes and host genetic components, as well as epigenetics and viruses associated with type 1 diabetes, with added emphasis on the enteroviruses, which are often associated with triggering the disease. Genus Enterovirus is classified into twelve species of which seven (Enterovirus A, Enterovirus B, Enterovirus C, and Enterovirus D and Rhinovirus A, Rhinovirus B, and Rhinovirus C) are human pathogens. These viruses are transmitted mainly by the fecal-oral route; they may also spread via the nasopharyngeal route. Enterovirus infections are highly prevalent, but these infections are usually subclinical or cause a mild flu-like illness. However, infections caused by enteroviruses can sometimes be serious, with manifestations of meningoencephalitis, paralysis, myocarditis, and in neonates a fulminant sepsis-like syndrome. These viruses are often implicated in chronic (inflammatory) diseases as chronic myocarditis, chronic pancreatitis, and type 1 diabetes. In this review we discuss the currently suggested mechanisms involved in the viral induction of type 1 diabetes. We recapitulate current basic knowledge and definitions.

HLA typing by next-generation sequencing - getting closer to reality

Next generation sequencing (NGS) denotes novel sequencing technologies that enable the generation of a large number of clonal sequences in a single sequencing run. NGS was initially introduced for whole genome sequencing and for quantitation of viral variants or genetic mutations in tumor tissues; more recently, the potential for high resolution HLA typing and high throughput analyses has been explored. It became clear that the complexity of the HLA system implicates new challenges, especially for bioinformatics. From an economical point of view, NGS is becoming increasingly attractive for HLA typing laboratories currently relying on Sanger based sequencing. Realizing the full potential of NGS will require the development of specifically adapted typing strategies and software algorithms. In the present review, three laboratories that were among the first to perform HLA-typing using different NGS platforms, the Roche 454, the Illumina Miseq and the Ion Torrent system, respectively, give an overview of these applications and point out advantages and limitations.

One year remission of type 1 diabetes mellitus in a patient treated with sitagliptin

Type 1 diabetes mellitus (T1DM) is a chronic disease characterized by the autoimmune destruction of pancreatic β-cells. This paper describes the case of a 19-year-old male patient who presented with glutamic acid decarboxylase (GAD) antibody positive and diabetic ketoacidosis, which mandated intensive insulin treatment. Once the ketoacidosis was controlled, an oral dose of 100 mg of sitagliptin was administered once a day. Ketoacidosis was managed by insulin and insulin daily requirement began to dwindle after one month, until its complete withdrawal at 8 weeks, when partial remission was reached. The patient has now remained on sitagliptin treatment alone for a year, without requiring insulin. The benefit observed with this medication is possibly associated with its immunological effects. Inhibition of dipeptidyl peptidase 4 in animal models deregulates the Th1 immune response, increases secretion of Th2 cytokines, activates CD4+CD25+FoxP3+ regulatory T-cells, and prevents IL17 production.

HLA-DRB1*07:01 is associated with a higher risk of asparaginase allergies

Asparaginase is a therapeutic enzyme used to treat leukemia and lymphoma, with immune responses resulting in suboptimal drug exposure and a greater risk of relapse. To elucidate whether there is a genetic component to the mechanism of asparaginase-induced immune responses, we imputed human leukocyte antigen (HLA) alleles in patients of European ancestry enrolled on leukemia trials at St. Jude Children's Research Hospital (n = 541) and the Children's Oncology Group (n = 1329). We identified a higher incidence of hypersensitivity and anti-asparaginase antibodies in patients with HLA-DRB1*07:01 alleles (P = 7.5 × 10(-5), odds ratio [OR] = 1.64; P = 1.4 × 10(-5), OR = 2.92, respectively). Structural analysis revealed that high-risk amino acids were located within the binding pocket of the HLA protein, possibly affecting the interaction between asparaginase epitopes and the HLA-DRB1 protein. Using a sequence-based consensus approach, we predicted the binding affinity of HLA-DRB1 alleles for asparaginase epitopes, and patients whose HLA genetics predicted high-affinity binding had more allergy (P = 3.3 × 10(-4), OR = 1.38). Our results suggest a mechanism of allergy whereby HLA-DRB1 alleles that confer high-affinity binding to asparaginase epitopes lead to a higher frequency of reactions. These trials were registered at www.clinicaltrials.gov as NCT00137111, NCT00549848, NCT00005603, and NCT00075725.

Full-length novel MHC class I allele discovery by next-generation sequencing: two platforms are better than one

Deep sequencing has revolutionized major histocompatibility complex (MHC) class I analysis of nonhuman primates by enabling high-throughput, economical, and comprehensive genotyping. Full-length MHC class I cDNA sequences, which are required to generate reagents such as MHC-peptide tetramers, cannot be directly obtained by short read deep sequencing. We combined data from two next-generation sequencing platforms to discover novel full-length MHC class I mRNA/cDNA transcripts in Chinese rhesus macaques. We first genotyped macaques by Roche/454 pyrosequencing using a 530-bp amplicon spanning the densely polymorphic exons 2 through 4 of the MHC class I loci that encode the peptide-binding region. We then mapped short paired-end 250 bp Illumina sequence reads spanning the full-length transcript to each 530-bp amplicon at high stringency and used paired-end information to reconstruct full-length allele sequences. We characterized 65 full-length sequences from six Chinese rhesus macaques. Overall, approximately 70 % of the alleles distinguished in these six animals contained new sequence information, including 29 novel transcripts. The flexibility of this approach should make full-length MHC class I allele genotyping accessible for any nonhuman primate population of interest. We are currently optimizing this method for full-length characterization of other highly polymorphic, duplicated loci such as the MHC class II DRB and killer immunoglobulin-like receptors. We anticipate that this method will facilitate rapid expansion and near completion of sequence libraries of polymorphic loci, such as MHC class I, within a few years.

The changing landscape of type 1 diabetes: recent developments and future frontiers

Type 1 diabetes (T1D) research has made great strides over the past decade with advances in understanding the pathogenesis, natural history, candidate environmental exposures, exposure triggering time, disease prediction, and diagnosis. Major monitoring efforts have provided baseline historical measures, leading to better epidemiological studies incorporating longitudinal biosamples (ie, biobanks), which have allowed for new technologies ('omics') to further expose the etiological agents responsible for the initiation, progression, and eventual clinical onset of T1D. These new frontiers have brought forth high-dimensionality data, which have furthered the evidence of the heterogeneous nature of T1D pathogenesis and allowed for a more mechanistic approach in understanding the etiology of T1D. This review will expand on the most recent advances in the quest for T1D determinants, drawing upon novel research tools that epidemiology, genetics, microbiology, and immunology have provided, linking them to the major hypotheses associated with T1D etiology, and discussing the future frontiers.

SeqSIMLA: a sequence and phenotype simulation tool for complex disease studies

Background

Association studies based on next-generation sequencing (NGS) technology have become popular, and statistical association tests for NGS data have been developed rapidly. A flexible tool for simulating sequence data in either unrelated case–control or family samples with different disease and quantitative trait models would be useful for evaluating the statistical power for planning a study design and for comparing power among statistical methods based on NGS data.

Results

We developed a simulation tool, SeqSIMLA, which can simulate sequence data with user-specified disease and quantitative trait models. We implemented two disease models, in which the user can flexibly specify the number of disease loci, effect sizes or population attributable risk, disease prevalence, and risk or protective loci. We also implemented a quantitative trait model, in which the user can specify the number of quantitative trait loci (QTL), proportions of variance explained by the QTL, and genetic models. We compiled recombination rates from the HapMap project so that genomic structures similar to the real data can be simulated.

Conclusions

SeqSIMLA can efficiently simulate sequence data with disease or quantitative trait models specified by the user. SeqSIMLA will be very useful for evaluating statistical properties for new study designs and new statistical methods using NGS. SeqSIMLA can be downloaded for free at http://seqsimla.sourceforge.net.