Base-resolution analyses of sequence and parent-of-origin dependent DNA methylation in the mouse genome.

Differential methylation of the two parental genomes in placental mammals is essential for genomic imprinting and embryogenesis. To systematically study this epigenetic process, we have generated a base-resolution, allele-specific DNA methylation (ASM) map in the mouse genome. We find parent-of-origin dependent (imprinted) ASM at 1,952 CG dinucleotides. These imprinted CGs form 55 discrete clusters including virtually all known germline differentially methylated regions (DMRs) and 23 previously unknown DMRs, with some occurring at microRNA genes. We also identify sequence-dependent ASM at 131,765 CGs. Interestingly, methylation at these sites exhibits a strong dependence on the immediate adjacent bases, allowing us to define a conserved sequence preference for the mammalian DNA methylation machinery. Finally, we report a surprising presence of non-CG methylation in the adult mouse brain, with some showing evidence of imprinting. Our results provide a resource for understanding the mechanisms of imprinting and allele-specific gene expression in mammalian cells.

Integrated analysis of single-cell chromatin state and transcriptome identified common vulnerability despite glioblastoma heterogeneity.

In 2021, the World Health Organization reclassified glioblastoma, the most common form of adult brain cancer, into isocitrate dehydrogenase (IDH)-wild-type glioblastomas and grade IV IDH mutant (G4 IDHm) astrocytomas. For both tumor types, intratumoral heterogeneity is a key contributor to therapeutic failure. To better define this heterogeneity, genome-wide chromatin accessibility and transcription profiles of clinical samples of glioblastomas and G4 IDHm astrocytomas were analyzed at single-cell resolution. These profiles afforded resolution of intratumoral genetic heterogeneity, including delineation of cell-to-cell variations in distinct cell states, focal gene amplifications, as well as extrachromosomal circular DNAs. Despite differences in IDH mutation status and significant intratumoral heterogeneity, the profiled tumor cells shared a common chromatin structure defined by open regions enriched for nuclear factor 1 transcription factors (NFIA and NFIB). Silencing of NFIA or NFIB suppressed in vitro and in vivo growths of patient-derived glioblastomas and G4 IDHm astrocytoma models. These findings suggest that despite distinct genotypes and cell states, glioblastoma/G4 astrocytoma cells share dependency on core transcriptional programs, yielding an attractive platform for addressing therapeutic challenges associated with intratumoral heterogeneity.

Identification of brain cell types underlying genetic association with word reading and correlated traits.

Neuroimaging studies implicate multiple cortical regions in reading ability/disability. However, the neural cell types integral to the reading process are unknown. To contribute to this gap in knowledge, we integrated genetic results from genome-wide association studies for word reading (n = 5054) with gene expression datasets from adult/fetal human brain. Linkage disequilibrium score regression (LDSC) suggested that variants associated with word reading were enriched in genes expressed in adult excitatory neurons, specifically layer 5 and 6 FEZF2 expressing neurons and intratelencephalic (IT) neurons, which express the marker genes LINC00507, THEMIS, or RORB. Inhibitory neurons (VIP, SST, and PVALB) were also found. This finding was interesting as neurometabolite studies previously implicated excitatory-inhibitory imbalances in the etiology of reading disabilities (RD). We also tested traits that shared genetic etiology with word reading (previously determined by polygenic risk scores): attention-deficit/hyperactivity disorder (ADHD), educational attainment, and cognitive ability. For ADHD, we identified enrichment in L4 IT adult excitatory neurons. For educational attainment and cognitive ability, we confirmed previous studies identifying multiple subclasses of adult cortical excitatory and inhibitory neurons, as well as astrocytes and oligodendrocytes. For educational attainment and cognitive ability, we also identified enrichment in multiple fetal cortical excitatory and inhibitory neurons, intermediate progenitor cells, and radial glial cells. In summary, this study supports a role of excitatory and inhibitory neurons in reading and excitatory neurons in ADHD and contributes new information on fetal cell types enriched in educational attainment and cognitive ability, thereby improving our understanding of the neurobiological basis of reading/correlated traits.

EGFR Mutation Promotes Glioblastoma through Epigenome and Transcription Factor Network Remodeling.

Epidermal growth factor receptor (EGFR) gene amplification and mutations are the most common oncogenic events in glioblastoma (GBM), but the mechanisms by which they promote aggressive tumor growth are not well understood. Here, through integrated epigenome and transcriptome analyses of cell lines, genotyped clinical samples, and TCGA data, we show that EGFR mutations remodel the activated enhancer landscape of GBM, promoting tumorigenesis through a SOX9 and FOXG1-dependent transcriptional regulatory network in vitro and in vivo. The most common EGFR mutation, EGFRvIII, sensitizes GBM cells to the BET-bromodomain inhibitor JQ1 in a SOX9, FOXG1-dependent manner. These results identify the role of transcriptional/epigenetic remodeling in EGFR-dependent pathogenesis and suggest a mechanistic basis for epigenetic therapy.

Language and reading impairments are associated with increased prevalence of non-right-handedness.

Handedness has been studied for association with language-related disorders because of its link with language hemispheric dominance. No clear pattern has emerged, possibly because of small samples, publication bias, and heterogeneous criteria across studies. Non-right-handedness (NRH) frequency was assessed in N = 2503 cases with reading and/or language impairment and N = 4316 sex-matched controls identified from 10 distinct cohorts (age range 6-19 years old; European ethnicity) using a priori set criteria. A meta-analysis (Ncases  = 1994) showed elevated NRH % in individuals with language/reading impairment compared with controls (OR = 1.21, CI = 1.06-1.39, p = .01). The association between reading/language impairments and NRH could result from shared pathways underlying brain lateralization, handedness, and cognitive functions.

Partitioning the heritability of Tourette syndrome and obsessive compulsive disorder reveals differences in genetic architecture.

The direct estimation of heritability from genome-wide common variant data as implemented in the program Genome-wide Complex Trait Analysis (GCTA) has provided a means to quantify heritability attributable to all interrogated variants. We have quantified the variance in liability to disease explained by all SNPs for two phenotypically-related neurobehavioral disorders, obsessive-compulsive disorder (OCD) and Tourette Syndrome (TS), using GCTA. Our analysis yielded a heritability point estimate of 0.58 (se = 0.09, p = 5.64e-12) for TS, and 0.37 (se = 0.07, p = 1.5e-07) for OCD. In addition, we conducted multiple genomic partitioning analyses to identify genomic elements that concentrate this heritability. We examined genomic architectures of TS and OCD by chromosome, MAF bin, and functional annotations. In addition, we assessed heritability for early onset and adult onset OCD. Among other notable results, we found that SNPs with a minor allele frequency of less than 5% accounted for 21% of the TS heritability and 0% of the OCD heritability. Additionally, we identified a significant contribution to TS and OCD heritability by variants significantly associated with gene expression in two regions of the brain (parietal cortex and cerebellum) for which we had available expression quantitative trait loci (eQTLs). Finally we analyzed the genetic correlation between TS and OCD, revealing a genetic correlation of 0.41 (se = 0.15, p = 0.002). These results are very close to previous heritability estimates for TS and OCD based on twin and family studies, suggesting that very little, if any, heritability is truly missing (i.e., unassayed) from TS and OCD GWAS studies of common variation. The results also indicate that there is some genetic overlap between these two phenotypically-related neuropsychiatric disorders, but suggest that the two disorders have distinct genetic architectures.

The Dyslexia-associated gene KIAA0319L is involved in neuronal migration in the developing chick visual system.

The gene KIAA0319-Like (KIAA0319L) is thought to confer susceptibility for developmental dyslexia. Dyslexia may be caused by alterations in neuronal migration, and in utero knockdown of KIAA0319L in rats indicated migration errors. However, studies carried out with KIAA0319L knockout mice did not reveal an altered neuronal migration phenotype. Gene knockout may activate compensatory mechanisms to buffer against genetic mutations during development. Here we assessed the role of KIAA0319L on migrating neurons in the chick developing tectum. Whole mount in situ hybridization was performed for KIAA0319L on embryonic day (E)3 - E5 chick embryos and in situ hybridization on sections was performed at later stages. The specificity and efficiency of engineered microRNA (miRNA) constructs targeting KIAA0319L for knocking down KIAA0319L were verified. miRNAs were electroporated into E5 chick optic tecta. Our studies demonstrate that KIAA0319L is expressed in the developing chick visual system, as well as in the otic vesicles. Knockdown of KIAA0319L in the optic tectum results in abnormal neuronal migration, strengthening the argument that KIAA0319L is involved in this developmental process.

Hypothesis-driven genome-wide association studies provide novel insights into genetics of reading disabilities.

Reading Disability (RD) is often characterized by difficulties in the phonology of the language. While the molecular mechanisms underlying it are largely undetermined, loci are being revealed by genome-wide association studies (GWAS). In a previous GWAS for word reading (Price, 2020), we observed that top single-nucleotide polymorphisms (SNPs) were located near to or in genes involved in neuronal migration/axon guidance (NM/AG) or loci implicated in autism spectrum disorder (ASD). A prominent theory of RD etiology posits that it involves disturbed neuronal migration, while potential links between RD-ASD have not been extensively investigated. To improve power to identify associated loci, we up-weighted variants involved in NM/AG or ASD, separately, and performed a new Hypothesis-Driven (HD)-GWAS. The approach was applied to a Toronto RD sample and a meta-analysis of the GenLang Consortium. For the Toronto sample (n = 624), no SNPs reached significance; however, by gene-set analysis, the joint contribution of ASD-related genes passed the threshold (p~1.45 × 10-2, threshold = 2.5 × 10-2). For the GenLang Cohort (n = 26,558), SNPs in DOCK7 and CDH4 showed significant association for the NM/AG hypothesis (sFDR q = 1.02 × 10-2). To make the GenLang dataset more similar to Toronto, we repeated the analysis restricting to samples selected for reading/language deficits (n = 4152). In this GenLang selected subset, we found significant association for a locus intergenic between BTG3-C21orf91 for both hypotheses (sFDR q < 9.00 × 10-4). This study contributes candidate loci to the genetics of word reading. Data also suggest that, although different variants may be involved, alleles implicated in ASD risk may be found in the same genes as those implicated in word reading. This finding is limited to the Toronto sample suggesting that ascertainment influences genetic associations.

Language Difficulties in School-Age Children With Developmental Dyslexia.

Developmental dyslexia (DD) is a common reading disability, affecting 5% to 11% of children in North America. Children classified as having DD often have a history of early language delay (ELD) or language impairments. Nevertheless, studies have reported conflicting results as to the association between DD-ELD and the extent of current language difficulties in children with DD. To examine these relationships, we queried the parents of school-age children with reading difficulties on their child's early and current language ability. Siblings were also examined. Children were directly assessed using quantitative tests of language and reading skills. To compare this study with the literature, we divided the sample (N = 674) into three groups: DD, intermediate readers (IR), and skilled readers (SR). We found a significant association between DD and ELD, with parents of children in the DD/IR groups reporting their children put words together later than the SR group. We also found a significant association between DD and language difficulties, with children with low reading skills having low expressive/receptive language abilities. Finally, we identified early language predicted current language, which predicted reading skills. These data contribute to research indicating that children with DD experience language difficulties, suggesting early recognition may help identify reading problems.

Genetic variation in the KIAA0319 5' region as a possible contributor to dyslexia.

Reading disabilities (RD) have been linked and associated with markers on chromosome 6p with results from multiple independent samples pointing to KIAA0319 as a risk gene and specifically, the 5' region of this gene. Here we focus genetic studies on a 2.3 kb region spanning the predicted promoter, the first untranslated exon, and part of the first intron, a region we identified as a region of open chromatin. Using DNA from probands with RD, we screened for genetic variants and tested select variants for association. We identified 17 DNA variants in this sample of probands, 16 of which were previously reported in public databases and one previously identified in a screen of this region. Based on the allele frequencies in the probands compared to public databases, and on possible functional consequences of the variation, we selected seven variants to test for association in a sample of families with RD, in addition to four variants which had been tested previously. We also tested two markers 5' of this region that were previously reported as associated. The strongest evidence for association was observed with alleles of the microsatellite marker located in the first untranslated exon and haplotypes of that marker. These results support previous studies indicating the 5' region of the KIAA0319 gene as the location of risk alleles contributing to RD.

Synaptic processes and immune-related pathways implicated in Tourette syndrome.

Tourette syndrome (TS) is a neuropsychiatric disorder of complex genetic architecture involving multiple interacting genes. Here, we sought to elucidate the pathways that underlie the neurobiology of the disorder through genome-wide analysis. We analyzed genome-wide genotypic data of 3581 individuals with TS and 7682 ancestry-matched controls and investigated associations of TS with sets of genes that are expressed in particular cell types and operate in specific neuronal and glial functions. We employed a self-contained, set-based association method (SBA) as well as a competitive gene set method (MAGMA) using individual-level genotype data to perform a comprehensive investigation of the biological background of TS. Our SBA analysis identified three significant gene sets after Bonferroni correction, implicating ligand-gated ion channel signaling, lymphocytic, and cell adhesion and transsynaptic signaling processes. MAGMA analysis further supported the involvement of the cell adhesion and trans-synaptic signaling gene set. The lymphocytic gene set was driven by variants in FLT3, raising an intriguing hypothesis for the involvement of a neuroinflammatory element in TS pathogenesis. The indications of involvement of ligand-gated ion channel signaling reinforce the role of GABA in TS, while the association of cell adhesion and trans-synaptic signaling gene set provides additional support for the role of adhesion molecules in neuropsychiatric disorders. This study reinforces previous findings but also provides new insights into the neurobiology of TS.

No evidence of an association between two genes, EDN1 and ACE, and childhood-onset mood disorders.

Recent evidence supports a pathological link between heart disease and depressive symptoms, suggesting that depression is both etiologic and prognostic to heart disease. Thus, biological molecules which are at the interface between heart and mind are plausible candidate genes for depressive disorders. To investigate this line of enquiry we have investigated two genes, Endothelin 1 (EDN1) and Angiotensin-converting enzyme (ACE) in a family-based sample with childhood-onset mood disorders (COMDs). EDN1 is highly expressed in endothelium where it acts as a potent vasoconstrictor, and is also expressed in the brain where it exhibits neurotransmitter characteristics. ACE acts as a potent vasopressor, and interacts with the hypothalamic-pituitary-adrenocortical (HPA) system, which is often dysregulated in mood disorders. Furthermore, ACE has recently been found to be associated with major depression. Polymorphisms were selected to best capture the genetic variation at the two loci, and to replicate previous associations. The markers were genotyped across EDN1 and ACE in a sample comprised of 382 Hungarian nuclear families ascertained through affected probands diagnosed with a mood disorders before the age of 15. We found no evidence of association between either of these genes and COMD. Consequently, we were unable to support our hypothesis that these two genes, which are involved in both vascular and brain functions are contributing to the susceptibility to mood disorders of children/adolescents.

Association of reading disabilities with regions marked by acetylated H3 histones in KIAA0319.

Reading disabilities (RDs) have been associated with chromosome 6p with recent studies pointing to two genes, DCDC2 and KIAA0319. In this study, markers across the 6p region were tested for association with RD. Our strongest findings were for association with markers in KIAA0319, although with the opposite alleles compared with a previous study. We also found association with markers in VMP, but not with DCDC2. Current evidence indicates that differential regulation of KIAA0319 and DCDC2 contributes to RD, thus we used chromatin immunoprecipitation coupled with genomic tiling arrays (ChIP-chip) to map acetylated histones, a molecular marker for regulatory elements, across a 500 kb genomic region covering the RD locus on 6p. This approach identified several regions marked by acetylated histones that mapped near associated markers, including intron 7 of DCDC2 and the 5' region of KIAA0319. The latter is located within the 70 kb region previously associated with differential expression of KIAA0319. Interestingly, five markers associated with RD in independent studies were also located within the 2.7 kb acetylated region, and six additional associated markers, including the most significant one in this study, were located within a 22 kb haplotype block that encompassed this region. Our data indicates that this putative regulatory region is a likely site of genetic variation contributing to RD in our sample, further narrowing the candidate region.

DNA Variant in the RPGRIP1L Gene Influences Alternative Splicing.

The retinitis pigmentosa GTPase regulator interacting protein 1-like (RPGRIP1L) gene encodes a ciliary protein that is critical for processes related to brain development, including development of left-right asymmetry, sonic hedgehog signaling, and neural tube formation. RPGRIP1L is a risk factor for retinal degeneration, and rare, deleterious variants in the RPGRIP1L gene cause Joubert syndrome and Meckel syndrome, both autosomal recessive disorders. These syndromes are characterized by dysfunctional primary cilia that result in abnormal development - and even lethality in the case of Meckel syndrome. Genetic studies have also implicated RPGRIP1L in psychiatric disorders by suggestive findings from genome-wide association studies and findings from rare-variant exome analyses for bipolar disorder and de novo mutations in autism. In this study we identify a common variant in RPGRIP1L, rs7203525, that influences alternative splicing, increasing the inclusion of exon 20 of RPGRIP1L. We detected this alternative splicing association in human postmortem brain tissue samples and, using a minigene assay combined with in vitro mutagenesis, confirmed that the alternative splicing is attributable to the alleles of this variant. The predominate RPGRIP1L isoform expressed in adult brains does not contain exon 20; thus, a shift to include this exon may impact brain function.

Genome-wide association study of word reading: Overlap with risk genes for neurodevelopmental disorders.

Reading disabilities (RD) are the most common neurocognitive disorder, affecting 5% to 17% of children in North America. These children often have comorbid neurodevelopmental/psychiatric disorders, such as attention deficit/hyperactivity disorder (ADHD). The genetics of RD and their overlap with other disorders is incompletely understood. To contribute to this, we performed a genome-wide association study (GWAS) for word reading. Then, using summary statistics from neurodevelopmental/psychiatric disorders, we computed polygenic risk scores (PRS) and used them to predict reading ability in our samples. This enabled us to test the shared aetiology between RD and other disorders. The GWAS consisted of 5.3 million single nucleotide polymorphisms (SNPs) and two samples; a family-based sample recruited for reading difficulties in Toronto (n = 624) and a population-based sample recruited in Philadelphia [Philadelphia Neurodevelopmental Cohort (PNC)] (n = 4430). The Toronto sample SNP-based analysis identified suggestive SNPs (P ~ 5 × 10-7 ) in the ARHGAP23 gene, which is implicated in neuronal migration/axon pathfinding. The PNC gene-based analysis identified significant associations (P < 2.72 × 10-6 ) for LINC00935 and CCNT1, located in the region of the KANSL2/CCNT1/LINC00935/SNORA2B/SNORA34/MIR4701/ADCY6 genes on chromosome 12q, with near significant SNP-based analysis. PRS identified significant overlap between word reading and intelligence (R2 = 0.18, P = 7.25 × 10-181 ), word reading and educational attainment (R2 = 0.07, P = 4.91 × 10-48 ) and word reading and ADHD (R2 = 0.02, P = 8.70 × 10-6 ; threshold for significance = 7.14 × 10-3 ). Overlap was also found between RD and autism spectrum disorder (ASD) as top-ranked genes were previously implicated in autism by rare and copy number variant analyses. These findings support shared risk between word reading, cognitive measures, educational outcomes and neurodevelopmental disorders, including ASD.

Further genetic evidence implicates the vasopressin system in childhood-onset mood disorders.

Studies in both animals and humans advocate a role for the vasopressin (AVP) system in the aetiology of depressive symptoms. Attention has particularly focused on the role of AVP in the overactivation of the hypothalamic-pituitary-adrenal (HPA)-axis in mood disorders. Elevated AVP plasma levels have been found in mood disorder patients, which are often positively correlated with the severity of symptoms. We recently reported an association between childhood-onset mood disorders (COMD) and polymorphisms in the receptor responsible for the AVP-mediated activation of the HPA-axis (AVPR1B). As genetic variation in the vasopressinergic system could provide a mechanism to explain the endocrine alterations observed in mood disorders, we investigated other genes in this system. The gene encoding AVP is the strongest candidate, particularly as genetic variation in this gene in rodents is associated with anxiety-related behaviours. Six single-nucleotide polymorphisms (SNPs) were genotyped across the AVP gene in a sample comprised of 586 Hungarian nuclear families ascertained through affected probands with a diagnosis of COMD. In addition, AVP coding and putative regulatory regions were screened for mutations using denaturing high-performance liquid chromatography. One SNP, 3' to the AVP, gene reached significance (P = 0.03), as did the overtransmission of a five-marker haplotype with a frequency of 22% (P = 0.0001). The subsequent mutation screen failed to identify any putative functional polymorphisms. The outcome of this study, combined with our previous association between COMD and AVPR1B, implicates genetic variation in vasopressinergic genes in mediating vulnerability to COMD.

Chromosome 5 and Gilles de la Tourette syndrome: Linkage in a large pedigree and association study of six candidates in the region.

Gilles de la Tourette Syndrome (TS) is a neuropsychiatric disorder characterized by both motor and vocal tics. In our previous genome scan for TS we identified evidence for linkage to the centromeric region of chromosome 5 in a single large family of 32 individuals with 10 family members with TS or chronic multiple tics (CMT). In this paper we report further analyses of the 5p-centromeric region in this pedigree. An additional 11 family members were identified and screened for TS. Using a set of 14 microsatellite markers we refined the linked region to a approximately 28 Mb interval between the markers D5S1506 and D5S76. A set of six candidate genes located in this region were selected to be tested for genetic association with TS. These genes were GDNF, ITGA1, ISL1, FGF10, HCN1 and SLC1A3. The TDT statistic was used for the association tests in a sample of 171 independent nuclear families with 241 affected children with TS. We found no evidence for an association between TS and markers in these genes in this sample of families. This study represents the first efforts to narrow the linkage region in the extended pedigree and the first tests of candidate genes in the chromosome 5 region linked to TS.

Association of SLITRK1 to Gilles de la Tourette Syndrome.

Previously the Slit and Trk-like family member 1 (SLITRK1) gene was identified as a candidate gene for Gilles de la Tourette Syndrome (GTS) based on a patient that carried a chromosomal inversion on 13q, as well as the identification of two rare DNA variants in the SLITRK1 gene. Since that report, studies have tested for the two rare variants in GTS and either did not find them, or when found, they did not segregate with the disorder in families, casting doubt on the relationship of this gene to GTS. We tested for these two rare variants and genotyped three polymorphisms that tag the currently identified major haplotypes of this gene in a sample of 154 nuclear families with GTS. In addition, the entire coding region was screened for novel DNA variants. We did not find the two reported rare variants in any of the probands or siblings in these families. We did however find significant evidence for association of a single polymorphism and of haplotypes of the three tagging polymorphisms. These findings provide the first support for the original finding indicating SLITRK1 as a susceptibility gene for GTS and indicate that further study of this gene in GTS is warranted.

Tagging SNP association study of the IL-1beta gene (IL1B) and childhood-onset mood disorders.

Given substantial evidence for IL-1beta involvement in the etiology of depression, the IL1B gene is a strong candidate for involvement in susceptibility to depressive disorders. However, association studies investigating this, to date, have been limited to just two polymorphisms (rs1143627[-31T/C] and rs16944[-511C/T]) that constitute only a fraction of the genetic variation that is actually present across this gene in the population. Here, in a family-based association study of childhood-onset mood disorders (COMD), characterized by onset of depression before the age of 15, we have used a gene-wide approach, employing a panel of five tagging SNPs spanning the entire gene. Based on TDT analyses of both individual alleles and haplotypes, in a study sample of 646 families (with 782 affected children), none of the SNPs, including those implicated in transcriptional regulation of the gene, showed evidence for association with COMD. This is the largest and most comprehensive study of IL1B in relation to mood disorders that has been carried out, to date. The results do not support the involvement of IL1B as a major factor in genetic risk for early-onset mood disorders.

G72/G30 (DAOA) and juvenile-onset mood disorders.

The chromosome 13q region has been linked to bipolar disorder in a number of genome scans as well as focused linkage studies. Previously we identified linkage to the 13q32 region in a genome scan of 146 affected sibling pair families from Hungary with juvenile-onset mood disorders. Within this region are the overlapping genes G72/G30, with G72 now officially named as D-amino-acid oxidase activator (DAOA). This locus has been associated with panic disorder, schizophrenia, and bipolar disorder. In this study, we tested for association to 11 markers in these genes and mood disorders in a sample of 646 nuclear families identified with a proband with onset of a mood disorder before 14.9 years of age. We identified evidence for association to three markers within the gene (rs2391191, rs3918341, rs1935062), two of which had been associated with bipolar disorder in previous studies. When corrected for the number of markers tested, the results were no longer significant, however the prior evidence for association of this gene in multiple studies points to this gene as a potential contributor to juvenile-onset mood disorders.