Recent studies implicate chromatin modifiers in autism spectrum disorder (ASD) through the identification of recurrent de novo loss of function mutations in affected individuals. ASD risk genes are ...co-expressed in human midfetal cortex, suggesting that ASD risk genes converge in specific regulatory networks during neurodevelopment. To elucidate such networks, we identify genes targeted by CHD8, a chromodomain helicase strongly associated with ASD, in human midfetal brain, human neural stem cells (hNSCs) and embryonic mouse cortex. CHD8 targets are strongly enriched for other ASD risk genes in both human and mouse neurodevelopment, and converge in ASD-associated co-expression networks in human midfetal cortex. CHD8 knockdown in hNSCs results in dysregulation of ASD risk genes directly targeted by CHD8. Integration of CHD8-binding data into ASD risk models improves detection of risk genes. These results suggest loss of CHD8 contributes to ASD by perturbing an ancient gene regulatory network during human brain development.
Gene expression levels vary across developmental stage, cell type, and region in the brain. Genomic variants also contribute to the variation in expression, and some neuropsychiatric disorder loci ...may exert their effects through this mechanism. To investigate these relationships, we present BrainVar, a unique resource of paired whole-genome and bulk tissue RNA sequencing from the dorsolateral prefrontal cortex of 176 individuals across prenatal and postnatal development. Here we identify common variants that alter gene expression (expression quantitative trait loci eQTLs) constantly across development or predominantly during prenatal or postnatal stages. Both “constant” and “temporal-predominant” eQTLs are enriched for loci associated with neuropsychiatric traits and disorders and colocalize with specific variants. Expression levels of more than 12,000 genes rise or fall in a concerted late-fetal transition, with the transitional genes enriched for cell-type-specific genes and neuropsychiatric risk loci, underscoring the importance of cataloging developmental trajectories in understanding cortical physiology and pathology.
Display omitted
•Whole-genome and RNA sequencing across human prefrontal cortex development•Gene-specific developmental trajectories characterize the late-fetal transition•Identification of constant, prenatal-predominant, and postnatal-predominant eQTLs•Integrated analysis implicates genes in loci associated with educational attainment
Werling et al. analyze gene expression across the span of human cerebral cortical development and profile the trajectories of individual genes, coordinated groups of genes, and their relationships to disorders. Integration of genetic variation identifies quantitative trait loci that implicate specific genes in loci associated with neuropsychiatric traits and disorders.
Obsessive-compulsive disorder (OCD) is a chronic and severe psychiatric disorder for which effective treatment options are limited. Structural and functional neuroimaging studies have consistently ...implicated the orbitofrontal cortex (OFC) and striatum in the pathophysiology of the disorder. Recent genetic evidence points to involvement of components of the excitatory synapse in the etiology of OCD. However, the transcriptional alterations that could link genetic risk to known structural and functional abnormalities remain mostly unknown. To assess potential transcriptional changes in the OFC and two striatal regions (caudate nucleus and nucleus accumbens) of OCD subjects relative to unaffected comparison subjects, we sequenced messenger RNA transcripts from these brain regions. In a joint analysis of all three regions, 904 transcripts were differentially expressed between 7 OCD versus 8 unaffected comparison subjects. Region-specific analyses highlighted a smaller number of differences, which concentrated in caudate and nucleus accumbens. Pathway analyses of the 904 differentially expressed transcripts showed enrichment for genes involved in synaptic signaling, with these synapse-associated genes displaying lower expression in OCD subjects relative to unaffected comparison subjects. Finally, we estimated that cell type fractions of medium spiny neurons were lower whereas vascular cells and astrocyte fractions were higher in tissue of OCD subjects. Together, these data provide the first unbiased examination of differentially expressed transcripts in both OFC and striatum of OCD subjects. These transcripts encoded synaptic proteins more often than expected by chance, and thus implicate the synapse as a vulnerable molecular compartment for OCD.
Abstract Background Phenotypic heterogeneity in autism has long been conjectured to be a major hindrance to the discovery of genetic risk factors, leading to numerous attempts to stratify children ...based on phenotype to increase power of discovery studies. This approach, however, is based on the hypothesis that phenotypic heterogeneity closely maps to genetic variation, which has not been tested. Our study examines the impact of subphenotyping of a well-characterized autism spectrum disorder (ASD) sample on genetic homogeneity and the ability to discover common genetic variants conferring liability to ASD. Methods Genome-wide genotypic data of 2576 families from the Simons Simplex Collection were analyzed in the overall sample and phenotypic subgroups defined on the basis of diagnosis, IQ, and symptom profiles. We conducted a family-based association study, as well as estimating heritability and evaluating allele scores for each phenotypic subgroup. Results Association analyses revealed no genome-wide significant association signal. Subphenotyping did not increase power substantially. Moreover, allele scores built from the most associated single nucleotide polymorphisms, based on the odds ratio in the full sample, predicted case status in subsets of the sample equally well and heritability estimates were very similar for all subgroups. Conclusions In genome-wide association analysis of the Simons Simplex Collection sample, reducing phenotypic heterogeneity had at most a modest impact on genetic homogeneity. Our results are based on a relatively small sample, one with greater homogeneity than the entire population; if they apply more broadly, they imply that analysis of subphenotypes is not a productive path forward for discovering genetic risk variants in ASD.
Copy-number variants (CNVs) are a major contributor to the pathophysiology of autism spectrum disorders (ASDs), but the functional impact of CNVs remains largely unexplored. Because brain tissue is ...not available from most samples, we interrogated gene expression in lymphoblasts from 244 families with discordant siblings in the Simons Simplex Collection in order to identify potentially pathogenic variation. Our results reveal that the overall frequency of significantly misexpressed genes (which we refer to here as outliers) identified in probands and unaffected siblings does not differ. However, in probands, but not their unaffected siblings, the group of outlier genes is significantly enriched in neural-related pathways, including neuropeptide signaling, synaptogenesis, and cell adhesion. We demonstrate that outlier genes cluster within the most pathogenic CNVs (rare de novo CNVs) and can be used for the prioritization of rare CNVs of potentially unknown significance. Several nonrecurrent CNVs with significant gene-expression alterations are identified (these include deletions in chromosomal regions 3q27, 3p13, and 3p26 and duplications at 2p15), suggesting that these are potential candidate ASD loci. In addition, we identify distinct expression changes in 16p11.2 microdeletions, 16p11.2 microduplications, and 7q11.23 duplications, and we show that specific genes within the 16p CNV interval correlate with differences in head circumference, an ASD-relevant phenotype. This study provides evidence that pathogenic structural variants have a functional impact via transcriptome alterations in ASDs at a genome-wide level and demonstrates the utility of integrating gene expression with mutation data for the prioritization of genes disrupted by potentially pathogenic mutations.
Ulcerative colitis is a chronic inflammatory disease of the colon that presents as diarrhea and gastrointestinal bleeding. We performed a genome-wide association study using DNA samples from 1,052 ...individuals with ulcerative colitis and preexisting data from 2,571 controls, all of European ancestry. In an analysis that controlled for gender and population structure, ulcerative colitis loci attaining genome-wide significance and subsequent replication in two independent populations were identified on chromosomes 1p36 (rs6426833, combined P = 5.1 x 10(-13), combined odds ratio OR = 0.73) and 12q15 (rs1558744, combined P = 2.5 x 10(-12), combined OR = 1.35). In addition, combined genome-wide significant evidence for association was found in a region spanning BTNL2 to HLA-DQB1 on chromosome 6p21 (rs2395185, combined P = 1.0 x 10(-16), combined OR = 0.66) and at the IL23R locus on chromosome 1p31 (rs11209026, combined P = 1.3 x 10(-8), combined OR = 0.56; rs10889677, combined P = 1.3 x 10(-8), combined OR = 1.29).
DNA methylation (DNAm), the addition of a methyl group to a cytosine in DNA, plays an important role in the regulation of gene expression. Single-nucleotide polymorphisms (SNPs) associated with ...schizophrenia (SZ) by genome-wide association studies (GWAS) often influence local DNAm levels. Thus, DNAm alterations, acting through effects on gene expression, represent one potential mechanism by which SZ-associated SNPs confer risk. In this study, we investigated genome-wide DNAm in postmortem superior temporal gyrus from 44 subjects with SZ and 44 non-psychiatric comparison subjects using Illumina Infinium MethylationEPIC BeadChip microarrays, and extracted cell-type-specific methylation signals by applying tensor composition analysis. We identified SZ-associated differential methylation at 242 sites, and 44 regions containing two or more sites (FDR cutoff of q = 0.1) and determined a subset of these were cell-type specific. We found mitotic arrest deficient 1-like 1 (MAD1L1), a gene within an established GWAS risk locus, harbored robust SZ-associated differential methylation. We investigated the potential role of MAD1L1 DNAm in conferring SZ risk by assessing for colocalization among quantitative trait loci for methylation and gene transcripts (mQTLs and tQTLs) in brain tissue and GWAS signal at the locus using multiple-trait-colocalization analysis. We found that mQTLs and tQTLs colocalized with the GWAS signal (posterior probability >0.8). Our findings suggest that alterations in MAD1L1 methylation and transcription may mediate risk for SZ at the MAD1L1-containing locus. Future studies to identify how SZ-associated differential methylation affects MAD1L1 biological function are indicated.
Polygenic scores (PGSs) are quantitative metrics for predicting phenotypic values, such as human height or disease status. Some PGS methods require only summary statistics of a relevant genome-wide ...association study (GWAS) for their score. One such method is Lassosum, which inherits the model selection advantages of Lasso to select a meaningful subset of the GWAS single-nucleotide polymorphisms as predictors from their association statistics. However, even efficient scores like Lassosum, when derived from European-based GWASs, are poor predictors of phenotype for subjects of non-European ancestry; that is, they have limited portability to other ancestries. To increase the portability of Lassosum, when GWAS information and estimates of linkage disequilibrium are available for both ancestries, we propose Joint-Lassosum (JLS). In the simulation settings we explore, JLS provides more accurate PGSs compared to other methods, especially when measured in terms of fairness. In analyses of UK Biobank data, JLS was computationally more efficient but slightly less accurate than a Bayesian comparator, SDPRX. Like all PGS methods, JLS requires selection of predictors, which are determined by data-driven tuning parameters. We describe a new approach to selecting tuning parameters and note its relevance for model selection for any PGS. We also draw connections to the literature on algorithmic fairness and discuss how JLS can help mitigate fairness-related harms that might result from the use of PGSs in clinical settings. While no PGS method is likely to be universally portable, due to the diversity of human populations and unequal information content of GWASs for different ancestries, JLS is an effective approach for enhancing portability and reducing predictive bias.
Genetic scores of individuals can be used to predict their traits, such as eventual disease status. However, these scores work best in the ancestry from which they were developed. Lack of portability across ancestries violates fairness principles and could generate clinical harm. We propose methods to expand portability.
The Autism Sequencing Consortium identified 102 high-confidence autism spectrum disorder (ASD) genes, showing that individuals with ASD and with potentially damaging single nucleotide variation ...(pdSNV) in these genes had lower cognitive levels and delayed age at walking, when compared to ASD participants without pdSNV. Here, we made use of a Swedish sample of individuals with ASD (called PAGES, for Population-Based Autism Genetics & Environment Study) to evaluate the frequency of pdSNV and their impact on medical and psychiatric phenotypes, using an epidemiological frame and universal health reporting. We then combine findings with those for potentially damaging copy number variation (pdCNV).
SNV and CNV calls were generated from whole-exome sequencing and chromosome microarray data, respectively. Birth and medical register data were used to collect phenotypes.
Of 808 individuals assessed by sequencing, 69 (9%) had pdSNV in the 102 ASC genes, and 144 (18%) had pdSNV in the 102 ASC genes or in a larger set of curated neurodevelopmental genes (from the Deciphering Developmental Disorders study, the gene2phenotype database, and the Radboud University gene lists). Three or more individuals had pdSNV in GRIN2B, POGZ, SATB1, DYNC1H1, SCN8A, or CREBBP. In comparison, out of the 996 individuals from whom CNV were called, 105 (11%) carried one or more pdCNV, including four or more individuals with CNV in the recurrent 15q11q13, 22q11.2, and 16p11.2 loci. Carriers of pdSNV were more likely to have intellectual disability (ID) and epilepsy, while carriers of pdCNV showed increased rates of congenital anomalies and scholastic skill disorders. Carriers of either pdSNV or pdCNV were more likely to have ID, scholastic skill disorders, and epilepsy.
The cohort only included individuals with autistic disorder, the more severe form of ASD, and phenotypes are defined from medical registers. Not all genes studied are definitively ASD genes, and we did not have de novo information to aid in classification.
In this epidemiological sample, rare pdSNV were more common than pdCNV and the combined yield of potentially damaging variation was substantial at 27%. The results provide compelling rationale for the use of high-throughout sequencing as part of routine clinical workup for ASD and support the development of precision medicine in ASD.