Gene family information facilitates variant interpretation and identification of disease-associated genes in neurodevelopmental disorders

Abstract Background Classifying pathogenicity of missense variants represents a major challenge in clinical practice during the diagnoses of rare and genetic heterogeneous neurodevelopmental disorders (NDDs). While orthologous gene conservation is commonly employed in variant annotation, approximately 80% of known disease-associated genes belong to gene families. The use of gene family information for disease gene discovery and variant interpretation has not yet been investigated on a genome-wide scale. We empirically evaluate whether paralog-conserved or non-conserved sites in human gene families are important in NDDs. Methods Gene family information was collected from Ensembl. Paralog-conserved sites were defined based on paralog sequence alignments; 10,068 NDD patients and 2078 controls were statistically evaluated for de novo variant burden in gene families. Results We demonstrate that disease-associated missense variants are enriched at paralog-conserved sites across all disease groups and inheritance models tested. We developed a gene family de novo enrichment framework that identified 43 exome-wide enriched gene families including 98 de novo variant carrying genes in NDD patients of which 28 represent novel candidate genes for NDD which are brain expressed and under evolutionary constraint. Conclusion This study represents the first method to incorporate gene family information into a statistical framework to interpret variant data for NDDs and to discover new NDD-associated genes

Targeted sequencing of 351 candidate genes for epileptic encephalopathy in a large cohort of patients

Peer reviewe

Genetic variation in CFH predicts phenytoin-induced maculopapular exanthema in European-descent patients

Objective To characterize, among European and Han Chinese populations, the genetic predictors of maculopapular exanthema (MPE), a cutaneous adverse drug reaction common to antiepileptic drugs. Methods We conducted a case-control genome-wide association study of autosomal genotypes, including Class I and II human leukocyte antigen (HLA) alleles, in 323 cases and 1,321 drug-tolerant controls from epilepsy cohorts of northern European and Han Chinese descent. Results from each cohort were meta-analyzed. Results We report an association between a rare variant in the complement factor H–related 4 (CFHR4) gene and phenytoin-induced MPE in Europeans (p = 4.5 × 10–11; odds ratio [95% confidence interval] 7 [3.2–16]). This variant is in complete linkage disequilibrium with a missense variant (N1050Y) in the complement factor H (CFH) gene. In addition, our results reinforce the association between HLA-A*31:01 and carbamazepine hypersensitivity. We did not identify significant genetic associations with MPE among Han Chinese patients. Conclusions The identification of genetic predictors of MPE in CFHR4 and CFH, members of the complement factor H–related protein family, suggest a new link between regulation of the complement system alternative pathway and phenytoin-induced hypersensitivity in European-ancestral patients

Genome-wide association analysis of genetic generalized epilepsies implicates susceptibility loci at 1q43, 2p16.1, 2q22.3 and 17q21.32

Genetic generalized epilepsies (GGEs) have a lifetime prevalence of 0.3% and account for 20-30% of all epilepsies. Despite their high heritability of 80%, the genetic factors predisposing to GGEs remain elusive. To identify susceptibility variants shared across common GGE syndromes, we carried out a two-stage genome-wide association study (GWAS) including 3020 patients with GGEs and 3954 controls of European ancestry. To dissect out syndrome-related variants, we also explored two distinct GGE subgroups comprising 1434 patients with genetic absence epilepsies (GAEs) and 1134 patients with juvenile myoclonic epilepsy (JME). Joint Stage-1 and 2 analyses revealed genome-wide significant associations for GGEs at 2p16.1 (rs13026414, Pmeta = 2.5 × 10−9, OR[T] = 0.81) and 17q21.32 (rs72823592, Pmeta = 9.3 × 10−9, OR[A] = 0.77). The search for syndrome-related susceptibility alleles identified significant associations for GAEs at 2q22.3 (rs10496964, Pmeta = 9.1 × 10−9, OR[T] = 0.68) and at 1q43 for JME (rs12059546, Pmeta = 4.1 × 10−8, OR[G] = 1.42). Suggestive evidence for an association with GGEs was found in the region 2q24.3 (rs11890028, Pmeta = 4.0 × 10−6) nearby the SCN1A gene, which is currently the gene with the largest number of known epilepsy-related mutations. The associated regions harbor high-ranking candidate genes: CHRM3 at 1q43, VRK2 at 2p16.1, ZEB2 at 2q22.3, SCN1A at 2q24.3 and PNPO at 17q21.32. Further replication efforts are necessary to elucidate whether these positional candidate genes contribute to the heritability of the common GGE syndrome

16p11.2 600 kb Duplications confer risk for typical and atypical Rolandic epilepsy

Rolandic epilepsy (RE) is the most common idiopathic focal childhood epilepsy. Its molecular basis is largely unknown and a complex genetic etiology is assumed in the majority of affected individuals. The present study tested whether six large recurrent copy number variants at 1q21, 15q11.2, 15q13.3, 16p11.2, 16p13.11 and 22q11.2 previously associated with neurodevelopmental disorders also increase risk of RE. Our association analyses revealed a significant excess of the 600 kb genomic duplication at the 16p11.2 locus (chr16: 29.5-30.1 Mb) in 393 unrelated patients with typical (n = 339) and atypical (ARE; n = 54) RE compared with the prevalence in 65 046 European population controls (5/393 cases versus 32/65 046 controls; Fisher's exact test P = 2.83 × 10−6, odds ratio = 26.2, 95% confidence interval: 7.9-68.2). In contrast, the 16p11.2 duplication was not detected in 1738 European epilepsy patients with either temporal lobe epilepsy (n = 330) and genetic generalized epilepsies (n = 1408), suggesting a selective enrichment of the 16p11.2 duplication in idiopathic focal childhood epilepsies (Fisher's exact test P = 2.1 × 10−4). In a subsequent screen among children carrying the 16p11.2 600 kb rearrangement we identified three patients with RE-spectrum epilepsies in 117 duplication carriers (2.6%) but none in 202 carriers of the reciprocal deletion. Our results suggest that the 16p11.2 duplication represents a significant genetic risk factor for typical and atypical R

A genome-wide association study identifies risk alleles in plasminogen and P4HA2 associated with giant cell arteritis

Giant cell arteritis (GCA) is the most common form of vasculitis in individuals older than 50 years in Western countries. To shed light onto the genetic background influencing susceptibility for GCA, we performed a genome-wide association screening in a well-powered study cohort. After imputation, 1,844,133 genetic variants were analysed in 2,134 cases and 9,125 unaffected controls from ten independent populations of European ancestry. Our data confirmed HLA class II as the strongest associated region (independent signals: rs9268905, P = 1.94E-54, per-allele OR = 1.79; and rs9275592, P = 1.14E-40, OR = 2.08). Additionally, PLG and P4HA2 were identified as GCA risk genes at the genome-wide level of significance (rs4252134, P = 1.23E-10, OR = 1.28; and rs128738, P = 4.60E-09, OR = 1.32, respectively). Interestingly, we observed that the association peaks overlapped with different regulatory elements related to cell types and tissues involved in the pathophysiology of GCA. PLG and P4HA2 are involved in vascular remodelling and angiogenesis, suggesting a high relevance of these processes for the pathogenic mechanisms underlying this type of vasculitis

Individualised prediction of drug resistance and seizure recurrence after medication withdrawal in people with juvenile myoclonic epilepsy: A systematic review and individual participant data meta-analysis

Summary Background A third of people with juvenile myoclonic epilepsy (JME) are drug-resistant. Three-quarters have a seizure relapse when attempting to withdraw anti-seizure medication (ASM) after achieving seizure-freedom. It is currently impossible to predict who is likely to become drug-resistant and safely withdraw treatment. We aimed to identify predictors of drug resistance and seizure recurrence to allow for individualised prediction of treatment outcomes in people with JME. Methods We performed an individual participant data (IPD) meta-analysis based on a systematic search in EMBASE and PubMed – last updated on March 11, 2021 – including prospective and retrospective observational studies reporting on treatment outcomes of people diagnosed with JME and available seizure outcome data after a minimum one-year follow-up. We invited authors to share standardised IPD to identify predictors of drug resistance using multivariable logistic regression. We excluded pseudo-resistant individuals. A subset who attempted to withdraw ASM was included in a multivariable proportional hazards analysis on seizure recurrence after ASM withdrawal. The study was registered at the Open Science Framework (OSF; https://osf.io/b9zjc/). Findings  368) was predicted by an earlier age at the start of withdrawal, shorter seizure-free interval and more currently used ASMs, resulting in an average internal-external cross-validation concordance-statistic of 0·70 (95%CI 0·68–0·73). Interpretation We were able to predict and validate clinically relevant personalised treatment outcomes for people with JME. Individualised predictions are accessible as nomograms and web-based tools. Funding MING fonds

Identification of Srp9 as a febrile seizure susceptibility gene

Objective: Febrile seizures (FS) are the most common seizure type in young children. Complex FS are a risk factor for mesial temporal lobe epilepsy (mTLE). To identify new FS susceptibility genes we used a forward genetic strategy in mice and subsequently analyzed candidate genes in humans. Methods: We mapped a quantitative trait locus (QTL1) for hyperthermia-induced FS on mouse chromosome 1, containing the signal recognition particle 9 (Srp9) gene. Effects of differential Srp9 expression were assessed in vivo and in vitro. Hippocampal SRP9 expression and genetic association were analyzed in FS and mTLE patients. Results: Srp9 was differentially expressed between parental strains C57BL/6J and A/J. Chromosome substitution strain 1 (CSS1) mice exhibited lower FS susceptibility and Srp9 expression than C57BL/6J mice. In vivo knockdown of brain Srp9 reduced FS susceptibility. Mice with reduced Srp9 expression and FS susceptibility, exhibited reduced hippocampal AMPA and NMDA currents. Downregulation of neuronal Srp9 reduced surface expression of AMPA receptor subunit GluA1. mTLE patients with antecedent FS had higher SRP9 expression than patients without. SRP9 promoter SNP rs12403575(G/A) was genetically associated with FS and mTLE. Interpretation: Our findings identify SRP9 as a novel FS susceptibility gene and indicate that SRP9 conveys its effects through endoplasmic reticulum (ER)-dependent synthesis and trafficking of membrane proteins, such as glutamate receptors. Discovery of this new FS gene and mechanism may provide new leads for early diagnosis and treatment of children with complex FS at risk for mTLE

A Large-Scale Genetic Analysis Reveals a Strong Contribution of the HLA Class II Region to Giant Cell Arteritis Susceptibility

We conducted a large-scale genetic analysis on giant cell arteritis (GCA), a polygenic immune-mediated vasculitis. A case-control cohort, comprising 1,651 case subjects with GCA and 15,306 unrelated control subjects from six different countries of European ancestry, was genotyped by the Immunochip array. We also imputed HLA data with a previously validated imputation method to perform a more comprehensive analysis of this genomic region. The strongest association signals were observed in the HLA region, with rs477515 representing the highest peak (p = 4.05 × 10−40, OR = 1.73). A multivariate model including class II amino acids of HLA-DRβ1 and HLA-DQα1 and one class I amino acid of HLA-B explained most of the HLA association with GCA, consistent with previously reported associations of classical HLA alleles like HLA-DRB1∗04. An omnibus test on polymorphic amino acid positions highlighted DRβ1 13 (p = 4.08 × 10−43) and HLA-DQα1 47 (p = 4.02 × 10−46), 56, and 76 (both p = 1.84 × 10−45) as relevant positions for disease susceptibility. Outside the HLA region, the most significant loci included PTPN22 (rs2476601, p = 1.73 × 10−6, OR = 1.38), LRRC32 (rs10160518, p = 4.39 × 10−6, OR = 1.20), and REL (rs115674477, p = 1.10 × 10−5, OR = 1.63). Our study provides evidence of a strong contribution of HLA class I and II molecules to susceptibility to GCA. In the non-HLA region, we confirmed a key role for the functional PTPN22 rs2476601 variant and proposed other putative risk loci for GCA involved in Th1, Th17, and Treg cell function