Linear mixed effect models are powerful tools used to account for population structure in genome-wide association studies (GWASs) and estimate the genetic architecture of complex traits. However, fully-specified models are computationally demanding and common simplifications often lead to reduced power or biased inference. We describe Grid-LMM (https://github.com/deruncie/GridLMM), an extendable algorithm for repeatedly fitting complex linear models that account for multiple sources of heterogeneity, such as additive and non-additive genetic variance, spatial heterogeneity, and genotype-environment interactions. Grid-LMM can compute approximate (yet highly accurate) frequentist test statistics or Bayesian posterior summaries at a genome-wide scale in a fraction of the time compared to existing general-purpose methods. We apply Grid-LMM to two types of quantitative genetic analyses. The first is focused on accounting for spatial variability and non-additive genetic variance while scanning for QTL; and the second aims to identify gene expression traits affected by non-additive genetic variation. In both cases, modeling multiple sources of heterogeneity leads to new discoveries.

Simultaneous analysis of genetic associations with multiple phenotypes may reveal shared genetic susceptibility across traits (pleiotropy). For a locus exhibiting overall pleiotropy, it is important to identify which specific traits underlie this association. We propose a Bayesian meta-analysis approach (termed CPBayes) that uses summary-level data across multiple phenotypes to simultaneously measure the evidence of aggregate-level pleiotropic association and estimate an optimal subset of traits associated with the risk locus. This method uses a unified Bayesian statistical framework based on a spike and slab prior. CPBayes performs a fully Bayesian analysis by employing the Markov Chain Monte Carlo (MCMC) technique Gibbs sampling. It takes into account heterogeneity in the size and direction of the genetic effects across traits. It can be applied to both cohort data and separate studies of multiple traits having overlapping or non-overlapping subjects. Simulations show that CPBayes can produce higher accuracy in the selection of associated traits underlying a pleiotropic signal than the subset-based meta-analysis ASSET. We used CPBayes to undertake a genome-wide pleiotropic association study of 22 traits in the large Kaiser GERA cohort and detected six independent pleiotropic loci associated with at least two phenotypes. This includes a locus at chromosomal region 1q24.2 which exhibits an association simultaneously with the risk of five different diseases: Dermatophytosis, Hemorrhoids, Iron Deficiency, Osteoporosis and Peripheral Vascular Disease. We provide an R-package 'CPBayes' implementing the proposed method.

For many complex traits, gene regulation is likely to play a crucial mechanistic role. How the genetic architectures of complex traits vary between populations and subsequent effects on genetic prediction are not well understood, in part due to the historical paucity of GWAS in populations of non-European ancestry. We used data from the MESA (Multi-Ethnic Study of Atherosclerosis) cohort to characterize the genetic architecture of gene expression within and between diverse populations. Genotype and monocyte gene expression were available in individuals with African American (AFA, n = 233), Hispanic (HIS, n = 352), and European (CAU, n = 578) ancestry. We performed expression quantitative trait loci (eQTL) mapping in each population and show genetic correlation of gene expression depends on shared ancestry proportions. Using elastic net modeling with cross validation to optimize genotypic predictors of gene expression in each population, we show the genetic architecture of gene expression for most predictable genes is sparse. We found the best predicted gene in each population, TACSTD2 in AFA and CHURC1 in CAU and HIS, had similar prediction performance across populations with R2 > 0.8 in each population. However, we identified a subset of genes that are well-predicted in one population, but poorly predicted in another. We show these differences in predictive performance are due to allele frequency differences between populations. Using genotype weights trained in MESA to predict gene expression in independent populations showed that a training set with ancestry similar to the test set is better at predicting gene expression in test populations, demonstrating an urgent need for diverse population sampling in genomics. Our predictive models and performance statistics in diverse cohorts are made publicly available for use in transcriptome mapping methods at https://github.com/WheelerLab/DivPop.

Infantile spasms (IS) and Lennox-Gastaut syndrome (LGS) are epileptic encephalopathies characterized by early onset, intractable seizures, and poor developmental outcomes. De novo sequence mutations and copy number variants (CNVs) are causative in a subset of cases. We used exome sequence data in 349 trios with IS or LGS to identify putative de novo CNVs. We confirm 18 de novo CNVs in 17 patients (4.8%), 10 of which are likely pathogenic, giving a firm genetic diagnosis for 2.9% of patients. Confirmation of exome-predicted CNVs by array-based methods is still required due to false-positive rates of prediction algorithms. Our exome-based results are consistent with recent array-based studies in similar cohorts and highlight novel candidate genes for IS and LGS.

Emerging evidence indicates that epileptic encephalopathies are genetically highly heterogeneous, underscoring the need for large cohorts of well-characterized individuals to further define the genetic landscape. Through a collaboration between two consortia (EuroEPINOMICS and Epi4K/EPGP), we analyzed exome-sequencing data of 356 trios with the "classical" epileptic encephalopathies, infantile spasms and Lennox Gastaut syndrome, including 264 trios previously analyzed by the Epi4K/EPGP consortium. In this expanded cohort, we find 429 de novo mutations, including de novo mutations in DNM1 in five individuals and de novo mutations in GABBR2, FASN, and RYR3 in two individuals each. Unlike previous studies, this cohort is sufficiently large to show a significant excess of de novo mutations in epileptic encephalopathy probands compared to the general population using a likelihood analysis (p = 8.2 × 10(-4)), supporting a prominent role for de novo mutations in epileptic encephalopathies. We bring statistical evidence that mutations in DNM1 cause epileptic encephalopathy, find suggestive evidence for a role of three additional genes, and show that at least 12% of analyzed individuals have an identifiable causal de novo mutation. Strikingly, 75% of mutations in these probands are predicted to disrupt a protein involved in regulating synaptic transmission, and there is a significant enrichment of de novo mutations in genes in this pathway in the entire cohort as well. These findings emphasize an important role for synaptic dysregulation in epileptic encephalopathies, above and beyond that caused by ion channel dysfunction.

Objective

Classification of epilepsy into types and subtypes is important for both clinical care and research into underlying disease mechanisms. A quantitative, data-driven approach may augment traditional electroclinical classification and shed new light on existing classification frameworks.

Methods

We used latent class analysis, a statistical method that assigns subjects into groups called latent classes based on phenotypic elements, to classify individuals with common familial epilepsies from the Epi4K Multiplex Families study. Phenotypic elements included seizure types, seizure symptoms, and other elements of the medical history. We compared class assignments to traditional electroclinical classifications and assessed familial aggregation of latent classes.

Results

A total of 1120 subjects with epilepsy were assigned to five latent classes. Classes 1 and 2 contained subjects with generalized epilepsy, largely reflecting the distinction between absence epilepsies and younger onset (class 1) versus myoclonic epilepsies and older onset (class 2). Classes 3 and 4 contained subjects with focal epilepsies, and in contrast to classes 1 and 2, these did not adhere as closely to clinically defined focal epilepsy subtypes. Class 5 contained nearly all subjects with febrile seizures plus or unknown epilepsy type, as well as a few subjects with generalized epilepsy and a few with focal epilepsy. Family concordance of latent classes was similar to or greater than concordance of clinically defined epilepsy types.

Significance

Quantitative classification of epilepsy has the potential to augment traditional electroclinical classification by (1) combining some syndromes into a single class, (2) splitting some syndromes into different classes, (3) helping to classify subjects who could not be classified clinically, and (4) defining the boundaries of clinically defined classifications. This approach can guide future research, including molecular genetic studies, by identifying homogeneous sets of individuals that may share underlying disease mechanisms.

Epileptic encephalopathies (EEs) are the most clinically important group of severe early-onset epilepsies. Next-generation sequencing has highlighted the crucial contribution of de novo mutations to the genetic architecture of EEs as well as to their underlying genetic heterogeneity. Our previous whole-exome sequencing study of 264 parent-child trios revealed more than 290 candidate genes in which only a single individual had a de novo variant. We sought to identify additional pathogenic variants in a subset (n = 27) of these genes via targeted sequencing in an unsolved cohort of 531 individuals with a diverse range of EEs. We report 17 individuals with pathogenic variants in seven of the 27 genes, defining a genetic etiology in 3.2% of this unsolved cohort. Our results provide definitive evidence that de novo mutations in SLC1A2 and CACNA1A cause specific EEs and expand the compendium of clinically relevant genotypes for GABRB3. We also identified EEs caused by genetic variants in ALG13, DNM1, and GNAO1 and report a mutation in IQSEC2. Notably, recurrent mutations accounted for 7/17 of the pathogenic variants identified. As a result of high-depth coverage, parental mosaicism was identified in two out of 14 cases tested with mutant allelic fractions of 5%-6% in the unaffected parents, carrying significant reproductive counseling implications. These results confirm that dysregulation in diverse cellular neuronal pathways causes EEs, and they will inform the diagnosis and management of individuals with these devastating disorders.

Gene identification in epilepsy has mainly been limited to large families segregating genes of major effect and de novo mutations in epileptic encephalopathies. Many families that present with common non-acquired focal epilepsies and genetic generalized epilepsies remain unexplained. We assembled a cohort of 'genetically enriched' common epilepsies by collecting and phenotyping families containing multiple individuals with unprovoked seizures. We aimed to determine if specific clinical epilepsy features aggregate within families, and whether this segregation of phenotypes may constitute distinct 'familial syndromes' that could inform genomic analyses. Families with three or more individuals with unprovoked seizures were studied across multiple international centres. Affected individuals were phenotyped and classified according to specific electroclinical syndromes. Families were categorized based on syndromic groupings of affected family members, examined for pedigree structure and phenotypic patterns and, where possible, assigned specific familial epilepsy syndromes. A total of 303 families were assembled and analysed, comprising 1120 affected phenotyped individuals. Of the 303 families, 117 exclusively segregated generalized epilepsy, 62 focal epilepsy, and 22 were classified as genetic epilepsy with febrile seizures plus. Over one-third (102 families) were observed to have mixed epilepsy phenotypes: 78 had both generalized and focal epilepsy features within the same individual (n = 39), or within first or second degree relatives (n = 39). Among the genetic generalized epilepsy families, absence epilepsies were found to cluster within families independently of juvenile myoclonic epilepsy, and significantly more females were affected than males. Of the 62 familial focal epilepsy families, two previously undescribed familial focal syndrome patterns were evident: 15 families had posterior quadrant epilepsies, including seven with occipito-temporal localization and seven with temporo-parietal foci, and four families displayed familial focal epilepsy of childhood with multiple affected siblings that was suggestive of recessive inheritance. The findings suggest (i) specific patterns of syndromic familial aggregation occur, including newly recognized forms of familial focal epilepsy; (ii) although syndrome-specificity usually occurs in multiplex families, the one-third of families with features of both focal and generalized epilepsy is suggestive of shared genetic determinants; and (iii) patterns of features observed across families including pedigree structure, sex, and age of onset may hold clues for future gene identification. Such detailed phenotypic information will be invaluable in the conditioning and interpretation of forthcoming sequencing data to understand the genetic architecture and inter-relationships of the common epilepsy syndromes.

The 22q11.2 deletion syndrome (22q11.2DS) results from non-allelic homologous recombination between low-copy repeats termed LCR22. About 60%-70% of individuals with the typical 3 megabase (Mb) deletion from LCR22A-D have congenital heart disease, mostly of the conotruncal type (CTD), whereas others have normal cardiac anatomy. In this study, we tested whether variants in the hemizygous LCR22A-D region are associated with risk for CTDs on the basis of the sequence of the 22q11.2 region from 1,053 22q11.2DS individuals. We found a significant association (FDR p < 0.05) of the CTD subset with 62 common variants in a single linkage disequilibrium (LD) block in a 350 kb interval harboring CRKL. A total of 45 of the 62 variants were associated with increased risk for CTDs (odds ratio [OR) ranges: 1.64-4.75). Associations of four variants were replicated in a meta-analysis of three genome-wide association studies of CTDs in affected individuals without 22q11.2DS. One of the replicated variants, rs178252, is located in an open chromatin region and resides in the double-elite enhancer, GH22J020947, that is predicted to regulate CRKL (CRK-like proto-oncogene, cytoplasmic adaptor) expression. Approximately 23% of patients with nested LCR22C-D deletions have CTDs, and inactivation of Crkl in mice causes CTDs, thus implicating this gene as a modifier. Rs178252 and rs6004160 are expression quantitative trait loci (eQTLs) of CRKL. Furthermore, set-based tests identified an enhancer that is predicted to target CRKL and is significantly associated with CTD risk (GH22J020946, sequence kernal association test (SKAT) p = 7.21 × 10^-5) in the 22q11.2DS cohort. These findings suggest that variance in CTD penetrance in the 22q11.2DS population can be explained in part by variants affecting CRKL expression.

Schizophrenia occurs in about one in four individuals with 22q11.2 deletion syndrome (22q11.2DS). The aim of this International Brain and Behavior 22q11.2DS Consortium (IBBC) study was to identify genetic factors that contribute to schizophrenia, in addition to the ~20-fold increased risk conveyed by the 22q11.2 deletion. Using whole-genome sequencing data from 519 unrelated individuals with 22q11.2DS, we conducted genome-wide comparisons of common and rare variants between those with schizophrenia and those with no psychotic disorder at age ≥25 years. Available microarray data enabled direct comparison of polygenic risk for schizophrenia between 22q11.2DS and independent population samples with no 22q11.2 deletion, with and without schizophrenia (total n = 35,182). Polygenic risk for schizophrenia within 22q11.2DS was significantly greater for those with schizophrenia (p_adj = 6.73 × 10^-6). Novel reciprocal case-control comparisons between the 22q11.2DS and population-based cohorts showed that polygenic risk score was significantly greater in individuals with psychotic illness, regardless of the presence of the 22q11.2 deletion. Within the 22q11.2DS cohort, results of gene-set analyses showed some support for rare variants affecting synaptic genes. No common or rare variants within the 22q11.2 deletion region were significantly associated with schizophrenia. These findings suggest that in addition to the deletion conferring a greatly increased risk to schizophrenia, the risk is higher when the 22q11.2 deletion and common polygenic risk factors that contribute to schizophrenia in the general population are both present.