Clinical and genetic correlates of islet-autoimmune signatures in juvenile-onset type 1 diabetes.

AIMS/HYPOTHESIS: Heterogeneity in individuals with type 1 diabetes has become more generally appreciated, but has not yet been extensively and systematically characterised. Here, we aimed to characterise type 1 diabetes heterogeneity by creating immunological, genetic and clinical profiles for individuals with juvenile-onset type 1 diabetes in a cross-sectional study. METHODS: Participants were HLA-genotyped to determine HLA-DR-DQ risk, and SNP-genotyped to generate a non-HLA genetic risk score (GRS) based on 93 type 1 diabetes-associated SNP variants outside the MHC region. Islet autoimmunity was assessed as T cell proliferation upon stimulation with the beta cell antigens GAD65, islet antigen-2 (IA-2), preproinsulin (PPI) and defective ribosomal product of the insulin gene (INS-DRIP). Clinical parameters were collected retrospectively. RESULTS: Of 80 individuals, 67 had proliferation responses to one or more islet antigens, with vast differences in the extent of proliferation. Based on the multitude and amplitude of the proliferation responses, individuals were clustered into non-, intermediate and high responders. High responders could not be characterised entirely by enrichment for the highest risk HLA-DR3-DQ2/DR4-DQ8 genotype. However, high responders did have a significantly higher non-HLA GRS. Clinically, high T cell responses to beta cell antigens did not reflect in worsened glycaemic control, increased complications, development of associated autoimmunity or younger age at disease onset. The number of beta cell antigens that an individual responded to increased with disease duration, pointing to chronic islet autoimmunity and epitope spreading. CONCLUSIONS/INTERPRETATION: Collectively, these data provide new insights into type 1 diabetes disease heterogeneity and highlight the importance of stratifying patients on the basis of their genetic and autoimmune signatures for immunotherapy and personalised disease management.

Assessment of parental mosaicism in SCN1A-related epilepsy by single-molecule molecular inversion probes and next-generation sequencing.

BACKGROUND: Dravet syndrome is a severe genetic encephalopathy, caused by pathogenic variants in SCN1A. Low-grade parental mosaicism occurs in a substantial proportion of families (7%-13%) and has important implications for recurrence risks. However, parental mosaicism can remain undetected by methods regularly used in diagnostics. In this study, we use single-molecule molecular inversion probes (smMIP), a technique with high sensitivity for detecting low-grade mosaic variants and high cost-effectiveness, to investigate the incidence of parental mosaicism of SCN1A variants in a cohort of 90 families and assess the feasibility of this technique. METHODS: Deep sequencing of SCN1A was performed using smMIPs. False positive rates for each of the proband's pathogenic variants were determined in 145 unrelated samples. If parents showed corresponding variant alleles at a significantly higher rate than the established noise ratio, mosaicism was confirmed by droplet digital PCR (ddPCR). RESULTS: Sequence coverage of at least 100× at the location of the corresponding pathogenic variant was reached for 80 parent couples. The variant ratio was significantly higher than the established noise ratio in eight parent couples, of which four (5%) were regarded as true mosaics, based on ddPCR results. The false positive rate of smMIP analysis without ddPCR was therefore 50%. Three of these variants had previously been considered de novo in the proband by Sanger sequencing. CONCLUSION: smMIP technology combined withnext generation sequencing (NGS) performs better than Sanger sequencing in the detection of parental mosaicism. Because parental mosaicism has important implications for genetic counselling and recurrence risks, we stress the importance of implementing high-sensitivity NGS-based assays in standard diagnostics.

Mosaicism of de novo pathogenic SCN1A variants in epilepsy is a frequent phenomenon that correlates with variable phenotypes.

OBJECTIVE: Phenotypes caused by de novo SCN1A pathogenic variants are very variable, ranging from severely affected patients with Dravet syndrome to much milder genetic epilepsy febrile seizures plus cases. The most important determinant of disease severity is the type of variant, with variants that cause a complete loss of function of the SCN1A protein (α-subunit of the neuronal sodium channel Nav1.1) being detected almost exclusively in Dravet syndrome patients. However, even within Dravet syndrome disease severity ranges greatly, and consequently other disease modifiers must exist. A better prediction of disease severity is very much needed in daily practice to improve counseling, stressing the importance of identifying modifying factors in this patient group. We evaluated 128 participants with de novo, pathogenic SCN1A variants to investigate whether mosaicism, caused by postzygotic mutation, is a major modifier in SCN1A-related epilepsy. METHODS: Mosaicism was investigated by reanalysis of the pathogenic SCN1A variants using single molecule molecular inversion probes and next generation sequencing with high coverage. Allelic ratios of pathogenic variants were used to determine whether mosaicism was likely. Selected mosaic variants were confirmed by droplet digital polymerase chain reaction and sequencing of different tissues. Developmental outcome was classified based on available data on intelligence quotient and school functioning/education. RESULTS: Mosaicism was present for 7.5% of de novo pathogenic SCN1A variants in symptomatic patients. Mosaic participants were less severely affected than nonmosaic participants if only participants with truncating variants are considered (distribution of developmental outcome scores, Mann-Whitney U, P = .023). SIGNIFICANCE: Postzygotic mutation is a common phenomenon in SCN1A-related epilepsies. Participants with mosaicism have on average milder phenotypes, suggesting that mosaicism can be a major modifier of SCN1A-related diseases. Detection of mosaicism has important implications for genetic counseling and can be achieved by deep sequencing of unique reads.

A genome-wide association study identifies a functional ERAP2 haplotype associated with birdshot chorioretinopathy.

Birdshot chorioretinopathy (BSCR) is a rare form of autoimmune uveitis that can lead to severe visual impairment. Intriguingly, >95% of cases carry the HLA-A29 allele, which defines the strongest documented HLA association for a human disease. We have conducted a genome-wide association study in 96 Dutch and 27 Spanish cases, and 398 unrelated Dutch and 380 Spanish controls. Fine-mapping the primary MHC association through high-resolution imputation at classical HLA loci, identified HLA-A*29:02 as the principal MHC association (odds ratio (OR) = 157.5, 95% CI 91.6-272.6, P = 6.6 × 10(-74)). We also identified two novel susceptibility loci at 5q15 near ERAP2 (rs7705093; OR = 2.3, 95% CI 1.7-3.1, for the T allele, P = 8.6 × 10(-8)) and at 14q32.31 in the TECPR2 gene (rs150571175; OR = 6.1, 95% CI 3.2-11.7, for the A allele, P = 3.2 × 10(-8)). The association near ERAP2 was confirmed in an independent British case-control samples (combined meta-analysis P = 1.7 × 10(-9)). Functional analyses revealed that the risk allele of the polymorphism near ERAP2 is strongly associated with high mRNA and protein expression of ERAP2 in B cells. This study further defined an extremely strong MHC risk component in BSCR, and detected evidence for a novel disease mechanism that affects peptide processing in the endoplasmic reticulum.

Increased paternal age and the influence on burden of genomic copy number variation in the general population.

Genomic copy number variations (CNVs) and increased parental age are both associated with the risk to develop a variety of clinical neuropsychiatric disorders such as autism, schizophrenia and bipolar disorder. At the same time, it has been shown that the rate of transmitted de novo single nucleotide mutations is increased with paternal age. To address whether paternal age also affects the burden of structural genomic deletions and duplications, we examined various types of CNV burden in a large population sample from the Netherlands. Healthy participants with parental age information (n = 6,773) were collected at different University Medical Centers. CNVs were called with the PennCNV algorithm using Illumina genome-wide SNP array data. We observed no evidence in support of a paternal age effect on CNV load in the offspring. Our results were negative for global measures as well as several proxies for de novo CNV events in this unique sample. While recent studies suggest de novo single nucleotide mutation rate to be dominated by the age of the father at conception, our results strongly suggest that at the level of global CNV burden there is no influence of increased paternal age. While it remains possible that local genomic effects may exist for specific phenotypes, this study indicates that global CNV burden and increased father's age may be independent disease risk factors.

Long RNA Sequencing and Ribosome Profiling of Inflamed ß-Cells Reveal an Extensive Translatome Landscape.

Type 1 diabetes (T1D) is an autoimmune disease characterized by autoreactive T cell-mediated destruction of the insulin-producing pancreatic ß-cells. Increasing evidence suggest that the ß-cells themselves contribute to their own destruction by generating neoantigens through the production of aberrant or modified proteins that escape central tolerance. We recently demonstrated that ribosomal infidelity amplified by stress could lead to the generation of neoantigens in human ß-cells, emphasizing the participation of nonconventional translation events in autoimmunity, as occurring in cancer or virus-infected tissues. Using a transcriptome-wide profiling approach to map translation initiation start sites in human ß-cells under standard and inflammatory conditions, we identify a completely new set of polypeptides derived from noncanonical start sites and translation initiation within long noncoding RNA. Our data underline the extreme diversity of the ß-cell translatome and may reveal new functional biomarkers for ß-cell distress, disease prediction and progression, and therapeutic intervention in T1D.

Using genome-wide pathway analysis to unravel the etiology of complex diseases.

Several genome-wide association studies (GWAS) have been published on various complex diseases. Although, new loci are found to be associated with these diseases, still only very little of the genetic risk for these diseases can be explained. As GWAS are still underpowered to find small main effects, and gene-gene interactions are likely to play a role, the data might currently not be analyzed to its full potential. In this study, we evaluated alternative methods to study GWAS data. Instead of focusing on the single nucleotide polymorphisms (SNPs) with the highest statistical significance, we took advantage of prior biological information and tried to detect overrepresented pathways in the GWAS data. We evaluated whether pathway classification analysis can help prioritize the biological pathways most likely to be involved in the disease etiology. In this study, we present the various benefits and limitations of pathway-classification tools in analyzing GWAS data. We show multiple differences in outcome between pathway tools analyzing the same dataset. Furthermore, analyzing randomly selected SNPs always results in significantly overrepresented pathways, large pathways have a higher chance of becoming statistically significant and the bioinformatics tools used in this study are biased toward detecting well-defined pathways. As an example, we analyzed data from two GWAS on type 2 diabetes (T2D): the Diabetes Genetics Initiative (DGI) and the Wellcome Trust Case Control Consortium (WTCCC). Occasionally the results from the DGI and the WTCCC GWAS showed concordance in overrepresented pathways, but discordance in the corresponding genes. Thus, incorporating gene networks and pathway classification tools into the analysis can point toward significantly overrepresented molecular pathways, which cannot be picked up using traditional single-locus analyses. However, the limitations discussed in this study, need to be addressed before these methods can be widely used.

Modifier genes in SCN1A-related epilepsy syndromes.

BACKGROUND: SCN1A is one of the most important epilepsy-related genes, with pathogenic variants leading to a range of phenotypes with varying disease severity. Different modifying factors have been hypothesized to influence SCN1A-related phenotypes. We investigate the presence of rare and more common variants in epilepsy-related genes as potential modifiers of SCN1A-related disease severity. METHODS: 87 patients with SCN1A-related epilepsy were investigated. Whole-exome sequencing was performed by the Beijing Genomics Institute (BGI). Functional variants in 422 genes associated with epilepsy and/or neuronal excitability were investigated. Differences in proportions of variants between the epilepsy genes and four control gene sets were calculated, and compared to the proportions of variants in the same genes in the ExAC database. RESULTS: Statistically significant excesses of variants in epilepsy genes were observed in the complete cohort and in the combined group of mildly and severely affected patients, particularly for variants with minor allele frequencies of <0.05. Patients with extreme phenotypes showed much greater excesses of epilepsy gene variants than patients with intermediate phenotypes. CONCLUSION: Our results indicate that relatively common variants in epilepsy genes, which would not necessarily be classified as pathogenic, may play a large role in modulating SCN1A phenotypes. They may modify the phenotypes of both severely and mildly affected patients. Our results may be a first step toward meaningful testing of modifier gene variants in regular diagnostics for individual patients, to provide a better estimation of disease severity for newly diagnosed patients.

Influence of common SCN1A promoter variants on the severity of SCN1A-related phenotypes.

BACKGROUND: Pathogenic variants in SCN1A cause variable epilepsy disorders with different disease severities. We here investigate whether common variation in the promoter region of the unaffected SCN1A allele could reduce normal expression, leading to a decreased residual function of Nav1.1, and therefore to more severe clinical outcomes in patients affected by pathogenic SCN1A variants. METHODS: Five different SCN1A promoter-haplotypes were functionally assessed in SH-SY5Y cells using Firefly and Renilla luciferase assays. The SCN1A promoter region was analyzed in a cohort of 143 participants with SCN1A pathogenic variants. Differences in clinical features and outcomes between participants with and without common variants in the SCN1A promoter-region of their unaffected allele were investigated. RESULTS: All non-wildtype haplotypes showed a significant reduction in luciferase expression, compared to the wildtype promoter-region (65%-80%, p = 0.039-0.0023). No statistically significant differences in clinical outcomes were observed between patients with and without common promoter variants. However, patients with a wildtype promoter-haplotype on their unaffected SCN1A allele showed a nonsignificant trend for milder phenotypes. CONCLUSION: The nonsignificant observed trends in our study warrant replication studies in larger cohorts to explore the potential modifying role of these common SCN1A promoter-haplotypes.

Exome-Wide Association Analysis of Coronary Artery Disease in the Kingdom of Saudi Arabia Population.

Coronary Artery Disease (CAD) remains the leading cause of mortality worldwide. Mortality rates associated with CAD have shown an exceptional increase particularly in fast developing economies like the Kingdom of Saudi Arabia (KSA). Over the past twenty years, CAD has become the leading cause of death in KSA and has reached epidemic proportions. This rise is undoubtedly caused by fast urbanization that is associated with a life-style that promotes CAD. However, the question remains whether genetics play a significant role and whether genetic susceptibility is increased in KSA compared to the well-studied Western European populations. Therefore, we performed an Exome-wide association study (EWAS) in 832 patients and 1,076 controls of Saudi Arabian origin to test whether population specific, strong genetic risk factors for CAD exist, or whether the polygenic risk score for known genetic risk factors for CAD, lipids, and Type 2 Diabetes show evidence for an enriched genetic burden. Our results do not show significant associations for a single genetic locus. However, the heritability estimate for CAD for this population was high (h(2) = 0.53, S.E. = 0.1, p = 4e(-12)) and we observed a significant association of the polygenic risk score for CAD that demonstrates that the population of KSA, at least in part, shares the genetic risk associated to CAD in Western populations.

A gene co-expression network in whole blood of schizophrenia patients is independent of antipsychotic-use and enriched for brain-expressed genes.

Despite large-scale genome-wide association studies (GWAS), the underlying genes for schizophrenia are largely unknown. Additional approaches are therefore required to identify the genetic background of this disorder. Here we report findings from a large gene expression study in peripheral blood of schizophrenia patients and controls. We applied a systems biology approach to genome-wide expression data from whole blood of 92 medicated and 29 antipsychotic-free schizophrenia patients and 118 healthy controls. We show that gene expression profiling in whole blood can identify twelve large gene co-expression modules associated with schizophrenia. Several of these disease related modules are likely to reflect expression changes due to antipsychotic medication. However, two of the disease modules could be replicated in an independent second data set involving antipsychotic-free patients and controls. One of these robustly defined disease modules is significantly enriched with brain-expressed genes and with genetic variants that were implicated in a GWAS study, which could imply a causal role in schizophrenia etiology. The most highly connected intramodular hub gene in this module (ABCF1), is located in, and regulated by the major histocompatibility (MHC) complex, which is intriguing in light of the fact that common allelic variants from the MHC region have been implicated in schizophrenia. This suggests that the MHC increases schizophrenia susceptibility via altered gene expression of regulatory genes in this network.

ATG16L1 and IL23R are associated with inflammatory bowel diseases but not with celiac disease in the Netherlands.

BACKGROUND: Inflammatory bowel disease (IBD)--Crohn's disease (CD) and ulcerative colitis (UC)--and celiac disease are intestinal inflammatory disorders with a complex genetic background. Recently, two novel genes were found to be associated with IBD susceptibility. One, an uncommon coding variant (rs11209026) in the gene encoding for the interleukin-23 receptor (IL23R), conferred strong protection against CD. The other, rs2241880 in the autophagy-related 16-like 1 gene (ATG16L1), was associated with CD. We performed a case-control study for the association of IBD with IL23R and ATG16L1 in a Dutch cohort. We also looked at the association of IL23R and ATG16L1 with celiac disease. METHODS: Five hundred eighteen Dutch white IBD patients (311 CD and 207 UC, including 176 trios of patients with both parents), 508 celiac disease patients, and 893 healthy controls were studied for association with the rs11209026 (IL23R) and rs2241880 (ATG16L1) single nucleotide polymorphisms (SNP). RESULTS: The rs11209026 SNP in IL23R had a protective effect for IBD in the case-control analysis (odds ratio [OR] 0.19, 95% confidence interval [CI] 0.10-0.37, P= 6.6E-09). Both CD (OR 0.14, CI 0.06-0.37, P= 3.9E-07) and UC (OR 0.33, CI 0.15-0.73, P= 1.4E-03) were associated with IL23R. For ATG16L1, the rs2241880 SNP was associated with CD susceptibility (OR 1.36, CI 1.12-1.66, P= 0.0017). The population-attributable risk of carrying allele G is 0.24 and is 0.19 for homozygosity for allele G in CD. No association was found between IL23R or ATG16L1 and celiac disease. CONCLUSIONS: We confirmed the association of IL23R and ATG16L1 with CD susceptibility and also the association of IL23R with UC. We found IL23R and ATG16L1 were not associated with celiac disease susceptibility.