Loci nominally associated with autism from genome-wide analysis show enrichment of brain expression quantitative trait loci but not lymphoblastoid cell line expression quantitative trait loci

Apically localized PANX1 impacts neuroepithelial expansion in human cerebral organoids

Dysfunctional paracrine signaling through Pannexin 1 (PANX1) channels is linked to several adult neurological pathologies and emerging evidence suggests that PANX1 plays an important role in human brain development. It remains unclear how early PANX1 influences brain development, or how loss of PANX1 alters the developing human brain. Using a cerebral organoid model of early human brain development, we find that PANX1 is expressed at all stages of organoid development from neural induction through to neuroepithelial expansion and maturation. Interestingly, PANX1 cellular distribution and subcellular localization changes dramatically throughout cerebral organoid development. During neural induction, PANX1 becomes concentrated at the apical membrane domain of neural rosettes where it co-localizes with several apical membrane adhesion molecules. During neuroepithelial expansion, PANX1-/- organoids are significantly smaller than control and exhibit significant gene expression changes related to cell adhesion, WNT signaling and non-coding RNAs. As cerebral organoids mature, PANX1 expression is significantly upregulated and is primarily localized to neuronal populations outside of the ventricular-like zones. Ultimately, PANX1 protein can be detected in all layers of a 21-22 post conception week human fetal cerebral cortex. Together, these results show that PANX1 is dynamically expressed by numerous cell types throughout embryonic and early fetal stages of human corticogenesis and loss of PANX1 compromises neuroepithelial expansion due to dysregulation of cell-cell and cell-matrix adhesion, perturbed intracellular signaling, and changes to gene regulation.

Overlap in synaptic neurological condition susceptibility pathways and the neural pannexin 1 interactome revealed by bioinformatics analyses

Many neurological conditions exhibit synaptic impairments, suggesting mechanistic convergence. Additionally, the pannexin 1 (PANX1) channel and signaling scaffold is linked to several of these neurological conditions and is an emerging regulator of synaptic development and plasticity; however, its synaptic pathogenic contributions are relatively unexplored. To this end, we explored connections between synaptic neurodevelopmental disorder and neurodegenerative disease susceptibility genes discovered by genome-wide association studies (GWASs), and the neural PANX1 interactome (483 proteins) identified from mouse Neuro2a (N2a) cells. To identify shared susceptibility genes, we compared synaptic suggestive GWAS candidate genes amongst autism spectrum disorders, schizophrenia, Parkinson's disease, and Alzheimer's disease. To further probe PANX1 signaling pathways at the synapse, we used bioinformatics tools to identify PANX1 interactome signaling pathways and protein-protein interaction clusters. To shed light on synaptic disease mechanisms potentially linking PANX1 and these four neurological conditions, we performed additional cross-analyses between gene ontologies enriched for the PANX1 synaptic and disease-susceptibility gene sets. Finally, to explore the regional specificity of synaptic PANX1-neurological condition connections, we identified brain region-specific elevations of synaptic PANX1 interactome and GWAS candidate gene set transcripts. Our results confirm considerable overlap in risk genes for autism spectrum disorders and schizophrenia and identify potential commonalities in genetic susceptibility for neurodevelopmental disorders and neurodegenerative diseases. Our findings also pinpointed novel putative PANX1 links to synaptic disease-associated pathways, such as regulation of vesicular trafficking and proteostasis, warranting further validation.

Transcriptomic analysis reveals mitochondrial pathways associated with distinct adolescent behavioral phenotypes and stress response

Adolescent individuals exhibit great variability in cortical dynamics and behavioral outcomes. The developing adolescent brain is highly sensitive to social experiences and environmental insults, influencing how personality traits emerge. A distinct pattern of mitochondrial gene expression in the prefrontal cortex (PFC) during adolescence underscores the essential role of mitochondria in brain maturation and the development of mental illnesses. Mitochondrial features in certain brain regions account for behavioral differences in adulthood. However, it remains unclear whether distinct adolescent behavioral phenotypes and the behavioral consequences of early adolescent stress exposure in rats are accompanied by changes in PFC mitochondria-related genes and mitochondria respiratory chain capacity. We performed a behavioral characterization during late adolescence (postnatal day, PND 47-50), including naïve animals and a group exposed to stress from PND 31-40 (10 days of footshock and 3 restraint sessions) by z-normalized data from three behavioral domains: anxiety (light-dark box tests), sociability (social interaction test) and cognition (novel-object recognition test). Employing principal component analysis, we identified three clusters: naïve with higher-behavioral z-score (HBZ), naïve with lower-behavioral z-score (LBZ), and stressed animals. Genome-wide transcriptional profiling unveiled differences in the expression of mitochondria-related genes in both naïve LBZ and stressed animals compared to naïve HBZ. Genes encoding subunits of oxidative phosphorylation complexes were significantly down-regulated in both naïve LBZ and stressed animals and positively correlated with behavioral z-score of phenotypes. Our network topology analysis of mitochondria-associated genes found Ndufa10 and Cox6a1 genes as central identifiers for naïve LBZ and stressed animals, respectively. Through high-resolution respirometry analysis, we found that both naïve LBZ and stressed animals exhibited a reduced prefrontal phosphorylation capacity and redox dysregulation. Our findings identify an association between mitochondrial features and distinct adolescent behavioral phenotypes while also underscoring the detrimental functional consequences of adolescent stress on the PFC.

The intersection of genome, epigenome and social experience in autism spectrum disorder: Exploring modifiable pathways for intervention

The number of children diagnosed with autism spectrum disorder (ASD) has increased substantially over the past two decades. Current research suggests that both genetic and environmental risk factors are involved in the etiology of ASD. The goal of this paper is to examine how one specific environmental factor, early social experience, may be correlated with DNA methylation (DNAm) changes in genes associated with ASD. We present an innovative model which proposes that polygenic risk and changes in DNAm due to social experience may both contribute to the symptoms of ASD. Previous research on genetic and environmental factors implicated in the etiology of ASD will be reviewed, with an emphasis on the oxytocin receptor gene, which may be epigenetically altered by early social experience, and which plays a crucial role in social and cognitive development. Identifying an environmental risk factor for ASD (e.g., social experience) that could be modified via early intervention and which results in epigenetic (DNAm) changes, could transform our understanding of this condition, facilitate earlier identification of ASD, and guide early intervention efforts.

Multiple Recurrent Copy Number Variations (CNVs) in Chromosome 22 Including 22q11.2 Associated with Autism Spectrum Disorder

Introduction

Autism spectrum disorder (ASD) is a developmental disorder that can cause substantial social, communication, and behavioral challenges. Genetic factors play a significant role in ASD, where the risk of ASD has been increased for unclear reasons. Twin studies have shown important evidence of both genetic and environmental contributions in ASD, where the level of contribution of these factors has not been proven yet. It has been suggested that copy number variation (CNV) duplication and the deletion of many genes in chromosome 22 (Ch22) may have a strong association with ASD. This study screened the CNVs in Ch22 in autistic Saudi children and assessed the candidate gene in the CNVs region of Ch22 that is most associated with ASD.

Methods

This study included 15 autistic Saudi children as well as 4 healthy children as controls; DNA was extracted from samples and analyzed using array comparative genomic hybridization (aCGH) and DNA sequencing.

Results

The aCGH detected (in only 6 autistic samples) deletion and duplication in many regions of Ch22, including some critical genes. Moreover, DNA sequencing determined a genetic mutation in the TBX1 gene sequence in autistic samples. This study, carried out using aCGH, found that six autistic patients had CNVs in Ch22, and DNA sequencing revealed mutations in the TBX1 gene in autistic samples but none in the control.

Conclusion

CNV deletion and the duplication of the TBX1 gene could be related to ASD; therefore, this gene needs more analysis in terms of expression levels.

The Role of Genetics, Epigenetics, and the Environment in ASD: A Mini Review

According to recent findings, variances in autism spectrum disorder (ASD) risk factors might be determined by several factors, including molecular genetic variants. Accumulated evidence has also revealed the important role of biological and chemical pathways in ASD aetiology. In this paper, we assess several reviews with regard to their quality of evidence and provide a brief outline of the presumed mechanisms of the genetic, epigenetic, and environmental risk factors of ASD. We also review some of the critical literature, which supports the basis of each factor in the underlying and specific risk patterns of ASD. Finally, we consider some of the implications of recent research regarding potential molecular targets for future investigations.

Predicting the Risk Genes of Autism Spectrum Disorders

Autism spectrum disorder (ASD) refers to a wide spectrum of neurodevelopmental disorders that emerge during infancy and continue throughout a lifespan. Although substantial efforts have been made to develop therapeutic approaches, core symptoms persist lifelong in ASD patients. Identifying the brain temporospatial regions where the risk genes are expressed in ASD patients may help to improve the therapeutic strategies. Accordingly, this work aims to predict the risk genes of ASD and identify the temporospatial regions of the brain structures at different developmental time points for exploring the specificity of ASD gene expression in the brain that would help in possible ASD detection in the future. A dataset consisting of 13 developmental stages ranging from 8 weeks post-conception to 8 years from 26 brain structures was retrieved from the BrainSpan atlas. This work proposes a support vector machine–based risk gene prediction method ASD-Risk to distinguish the risk genes of ASD and non-ASD genes. ASD-Risk used an optimal feature selection algorithm called inheritable bi-objective combinatorial genetic algorithm to identify the brain temporospatial regions for prediction of the risk genes of ASD. ASD-Risk achieved a 10-fold cross-validation accuracy, sensitivity, specificity, area under a receiver operating characteristic curve, and a test accuracy of 81.83%, 0.84, 0.79, 0.84, and 72.27%, respectively. We prioritized the temporospatial features according to their contribution to the prediction accuracy. The top identified temporospatial regions of the brain for risk gene prediction included the posteroventral parietal cortex at 13 post-conception weeks feature. The identified temporospatial features would help to explore the risk genes that are specifically expressed in different brain regions of ASD patients.

Purinergic signaling in nervous system health and disease: Focus on pannexin 1

Purinergic signaling encompasses the cycle of adenosine 5' triphosphate (ATP) release and its metabolism into nucleotide and nucleoside derivatives, the direct release of nucleosides, and subsequent receptor-triggered downstream intracellular pathways. Since the discovery of nerve terminal and glial ATP release into the neuropil, purinergic signaling has been implicated in the modulation of nervous system development, function, and disease. In this review, we detail our current understanding of the roles of the pannexin 1 (PANX1) ATP-release channel in neuronal development and plasticity, glial signaling, and neuron-glial-immune interactions. We additionally provide an overview of PANX1 structure, activation, and permeability to orientate readers and highlight recent research developments. We identify areas of convergence between PANX1 and purinergic receptor actions. Additional highlights include data on PANX1's participation in the pathophysiology of nervous system developmental, degenerative, and inflammatory disorders. Our aim in combining this knowledge is to facilitate the movement of our current understanding of PANX1 in the context of other nervous system purinergic signaling mechanisms one step closer to clinical translation.

Pannexin 1 mutation found in melanoma tumor reduces phosphorylation, glycosylation, and trafficking of the channel-forming protein

Pannexin 1 (PANX1) is a glycoprotein that forms large pore channels capable of passing ions and metabolites such as ATP for cellular communication. PANX1 has been implicated in many diseases including breast cancer and melanoma, where inhibition or deletion of PANX1 reduced the tumorigenic and metastatic properties of the cancer cells. We interrogated the effect of single amino acid changes in various PANX1 domains using naturally occurring variants reported in cancer patient tumors. We found that a previously reported variant (Q5H) is present in cancer cells, but was not different from the wild type (Q5) in glycosylation, trafficking, or channel function and did not affect cellular properties. We discovered that the Q5H variant is in fact the highly conserved ancestral allele of PANX1 with 89% of humans carrying at least one Q5H allele. Another mutated form Y150F, found in a melanoma patient tumor, prevented phosphorylation at Y150 as well as complex N-glycosylation while increasing intracellular localization. Sarcoma (SRC) is the predicted kinase to phosphorylate the Y150 residue, and its phosphorylation is not likely to be constitutive, but rather dynamically regulated. The Y150 phosphorylation site is the first one reported to play a role in regulating posttranslational modifications and trafficking of PANX1, with potential consequences on its large-pore channel structure and function in melanoma cells.

Pannexin 1 Regulates Dendritic Protrusion Dynamics in Immature Cortical Neurons

The integration of neurons into networks relies on the formation of dendritic spines. These specialized structures arise from dynamic filopodia-like dendritic protrusions. It was recently reported that cortical neurons lacking the channel protein pannexin 1 (PANX1) exhibited higher dendritic spine densities. Here, we expanded on those findings to investigate, at an earlier developmental time point (with more abundant dendritic protrusions), whether differences in the properties of dendritic protrusion dynamics could contribute to this previously discovered phenomenon. Using a fluorescent membrane tag (mCherry-CD9-10) to visualize dendritic protrusions in developing neurons [at 10 d in vitro (DIV10)], we confirmed that lack of PANX1 led to higher protrusion density, while transient transfection of Panx1 led to decreased protrusion density. To quantify the impact of PANX1 expression on protrusion formation, elimination, and motility, we used live cell imaging in DIV10 neurons (one frame every 5 s for 10 min). We discovered that at DIV10, loss of PANX1 stabilized protrusions. Notably, re-expression of PANX1 in Panx1 knock-out (KO) neurons resulted in a significant increase in protrusion motility and turnover. In summary, these new data revealed that PANX1 could regulate the development of dendritic spines, in part, by controlling dendritic protrusion dynamics.

Pannexin 1 Regulates Network Ensembles and Dendritic Spine Development in Cortical Neurons

Dendritic spines are the postsynaptic targets of excitatory synaptic inputs that undergo extensive proliferation and maturation during the first postnatal month in mice. However, our understanding of the molecular mechanisms that regulate spines during this critical period is limited. Previous work has shown that pannexin 1 (Panx1) regulates neurite growth and synaptic plasticity. We therefore investigated the impact of global Panx1 KO on spontaneous cortical neuron activity using Ca²⁺ imaging and in silico network analysis. Panx1 KO increased both the number and size of spontaneous co-active cortical neuron network ensembles. To understand the basis for these findings, we investigated Panx1 expression in postnatal synaptosome preparations from early postnatal mouse cortex. Between 2 and 4 postnatal weeks, we observed a precipitous drop in cortical synaptosome protein levels of Panx1, suggesting it regulates synapse proliferation and/or maturation. At the same time points, we observed significant enrichment of the excitatory postsynaptic density proteins PSD-95, GluA1, and GluN2a in cortical synaptosomes from global Panx1 knock-out mice. Ex vivo analysis of pyramidal neuron structure in somatosensory cortex revealed a consistent increase in dendritic spine densities in both male and female Panx1 KO mice. Similar findings were observed in an excitatory neuron-specific Panx1 KO line (Emx1-Cre driven; Panx1 cKO^E) and in primary Panx1 KO cortical neurons cultured in vitro. Altogether, our study suggests that Panx1 negatively regulates cortical dendritic spine development.

Calculating the statistical significance of rare variants causal for Mendelian and complex disorders

BACKGROUND

With the expanding use of next-gen sequencing (NGS) to diagnose the thousands of rare Mendelian genetic diseases, it is critical to be able to interpret individual DNA variation. To calculate the significance of finding a rare protein-altering variant in a given gene, one must know the frequency of seeing a variant in the general population that is at least as damaging as the variant in question.

METHODS

We developed a general method to better interpret the likelihood that a rare variant is disease causing if observed in a given gene or genic region mapping to a described protein domain, using genome-wide information from a large control sample. Based on data from 2504 individuals in the 1000 Genomes Project dataset, we calculated the number of individuals who have a rare variant in a given gene for numerous filtering threshold scenarios, which may be used for calculating the significance of an observed rare variant being causal for disease. Additionally, we calculated mutational burden data on the number of individuals with rare variants in genic regions mapping to protein domains.

RESULTS

We describe methods to use the mutational burden data for calculating the significance of observing rare variants in a given proportion of sequenced individuals. We present SORVA, an implementation of these methods as a web tool, and we demonstrate application to 20 relevant but diverse next-gen sequencing studies. Specifically, we calculate the statistical significance of findings involving multi-family studies with rare Mendelian disease and a large-scale study of a complex disorder, autism spectrum disorder. If we use the frequency counts to rank genes based on intolerance for variation, the ranking correlates well with pLI scores derived from the Exome Aggregation Consortium (ExAC) dataset (ρ = 0.515), with the benefit that the scores are directly interpretable.

CONCLUSIONS

We have presented a strategy that is useful for vetting candidate genes from NGS studies and allows researchers to calculate the significance of seeing a variant in a given gene or protein domain. This approach is an important step towards developing a quantitative, statistics-based approach for presenting clinical findings.

Revealing the complex genetic architecture of obsessive-compulsive disorder using meta-analysis

Two obsessive-compulsive disorder (OCD) genome-wide association studies (GWASs) have been published by independent OCD consortia, the International Obsessive-Compulsive Disorder Foundation Genetics Collaborative (IOCDF-GC) and the OCD Collaborative Genetics Association Study (OCGAS), but many of the top-ranked signals were supported in only one study. We therefore conducted a meta-analysis from the two consortia, investigating a total of 2688 individuals of European ancestry with OCD and 7037 genomically matched controls. No single-nucleotide polymorphisms (SNPs) reached genome-wide significance. However, in comparison with the two individual GWASs, the distribution of P-values shifted toward significance. The top haplotypic blocks were tagged with rs4733767 (P=7.1 × 10^-7; odds ratio (OR)=1.21; confidence interval (CI): 1.12-1.31, CASC8/CASC11), rs1030757 (P=1.1 × 10^-6; OR=1.18; CI: 1.10-1.26, GRID2) and rs12504244 (P=1.6 × 10^-6; OR=1.18; CI: 1.11-1.27, KIT). Variants located in or near the genes ASB13, RSPO4, DLGAP1, PTPRD, GRIK2, FAIM2 and CDH20, identified in linkage peaks and the original GWASs, were among the top signals. Polygenic risk scores for each individual study predicted case-control status in the other by explaining 0.9% (P=0.003) and 0.3% (P=0.0009) of the phenotypic variance in OCGAS and the European IOCDF-GC target samples, respectively. The common SNP heritability in the combined OCGAS and IOCDF-GC sample was estimated to be 0.28 (s.e.=0.04). Strikingly, ∼65% of the SNP-based heritability in the OCGAS sample was accounted for by SNPs with minor allele frequencies of ⩾40%. This joint analysis constituting the largest single OCD genome-wide study to date represents a major integrative step in elucidating the genetic causes of OCD.

Learning-dependent chromatin remodeling highlights noncoding regulatory regions linked to autism

Autism spectrum disorder (ASD) is a prevalent neurodevelopmental disorder that is associated with genetic risk factors. Most human disease-associated single-nucleotide polymorphisms (SNPs) are not located in genes but rather are in regulatory regions that control gene expression. The function of regulatory regions is determined through epigenetic mechanisms. Parallels between the cellular basis of development and the formation of long-term memory have long been recognized, particularly the role of epigenetic mechanisms in both processes. We analyzed how learning alters chromatin accessibility in the mouse hippocampus using a new high-throughput sequencing bioinformatics strategy we call DEScan (differential enrichment scan). DEScan, which enabled the analysis of data from epigenomic experiments containing multiple replicates, revealed changes in chromatin accessibility at 2365 regulatory regions-most of which were promoters. Learning-regulated promoters were active during forebrain development in mice and were enriched in epigenetic modifications indicative of bivalent promoters. These promoters were disproportionally intronic, showed a complex relationship with gene expression and alternative splicing during memory consolidation and retrieval, and were enriched in the data set relative to known ASD risk genes. Genotyping in a clinical cohort within one of these promoters (SHANK3 promoter 6) revealed that the SNP rs6010065 was associated with ASD. Our data support the idea that learning recapitulates development at the epigenetic level and demonstrate that behaviorally induced epigenetic changes in mice can highlight regulatory regions relevant to brain disorders in patients.

Cross-tissue integration of genetic and epigenetic data offers insight into autism spectrum disorder

Integration of emerging epigenetic information with autism spectrum disorder (ASD) genetic results may elucidate functional insights not possible via either type of information in isolation. Here we use the genotype and DNA methylation (DNAm) data from cord blood and peripheral blood to identify SNPs associated with DNA methylation (meQTL lists). Additionally, we use publicly available fetal brain and lung meQTL lists to assess enrichment of ASD GWAS results for tissue-specific meQTLs. ASD-associated SNPs are enriched for fetal brain (OR = 3.55; P < 0.001) and peripheral blood meQTLs (OR = 1.58; P < 0.001). The CpG targets of ASD meQTLs across cord, blood, and brain tissues are enriched for immune-related pathways, consistent with other expression and DNAm results in ASD, and reveal pathways not implicated by genetic findings. This joint analysis of genotype and DNAm demonstrates the potential of both brain and blood-based DNAm for insights into ASD and psychiatric phenotypes more broadly.

“There have been a number of recent epigenetic studies on autism spectrum disorder. Here, the authors integrate genetic and epigenetic data from cord and peripheral blood and also from brain tissues to show the potential of blood-based epigenetic data to provide insights into psychiatric disorders.”

Integrative genetic analysis suggests that skin color modifies the genetic architecture of melanoma

Melanoma is the deadliest form of skin cancer and presents a significant health care burden in many countries. In addition to ultraviolet radiation in sunlight, the main causal factor for melanoma, genetic factors also play an important role in melanoma susceptibility. Although genome-wide association studies have identified many single nucleotide polymorphisms associated with melanoma, little is known about the proportion of disease risk attributable to these loci and their distribution throughout the genome. Here, we investigated the genetic architecture of melanoma in 1,888 cases and 990 controls of European non-Hispanic ancestry. We estimated the overall narrow-sense heritability of melanoma to be 0.18 (P < 0.03), indicating that genetics contributes significantly to the risk of sporadically-occurring melanoma. We then demonstrated that only a small proportion of this risk is attributable to known risk variants, suggesting that much remains unknown of the role of genetics in melanoma. To investigate further the genetic architecture of melanoma, we partitioned the heritability by chromosome, minor allele frequency, and functional annotations. We showed that common genetic variation contributes significantly to melanoma risk, with a risk model defined by a handful of genomic regions rather than many risk loci distributed throughout the genome. We also demonstrated that variants affecting gene expression in skin account for a significant proportion of the heritability, and are enriched among melanoma risk loci. Finally, by incorporating skin color into our analyses, we observed both a shift in significance for melanoma-associated loci and an enrichment of expression quantitative trait loci among melanoma susceptibility variants. These findings suggest that skin color may be an important modifier of melanoma risk. We speculate that incorporating skin color and other non-genetic factors into genetic studies may allow for an improved understanding of melanoma susceptibility and guide future investigations to identify melanoma risk genes.

Transcriptional and post-translational regulation of pannexins

Pannexins are a 3-membered family of proteins that form large pore ion and metabolite channels in vertebrates. The impact of pannexins on vertebrate biology is intricately tied to where and when they are expressed, and how they are modified, once produced. The purpose of this review is therefore to outline our current understanding of transcriptional and post-translational regulation of pannexins. First, we briefly summarize their discovery and characteristics. Next, we describe several aspects of transcriptional regulation, including cell and tissue-specific expression, dynamic expression over development and disease, as well as new insights into the underlying molecular machinery involved. Following this, we delve into the role of post-translational modifications in the regulation of trafficking and channel properties, highlighting important work on glycosylation, phosphorylation, S-nitrosylation and proteolytic cleavage. Embedded throughout, we also highlight important knowledge gaps and avenues of future research. This article is part of a Special Issue entitled: Gap Junction Proteins edited by Jean Claude Herve.

The role of Pannexin gene variants in schizophrenia: systematic analysis of phenotypes

Pannexins are a group of brain-expressed channel proteins thought to be regulators of schizophrenia-linked pathways including glutamate release, synaptic plasticity and neural stem proliferation. We got evidence for linkage of a catatonic phenotype to the PANX2 locus in a family study. Aim of our study was to evaluate the role of Pannexins in schizophrenia and clinical phenotypes, particularly with regard to periodic catatonia. We genotyped six single-nucleotide polymorphisms at PANX1, five at PANX2 and three at PANX3 in 1173 German cases with schizophrenia according to DSM-5 and 480 controls. Our sample included 338 cases with periodic catatonia corresponding to Leonhard's classification. Association with schizophrenia according to DSM-5 was limited to genotype rs4838858-TT [p = 0.02, odds ratio (OR) 3.1] and haplotype rs4838858T-rs5771206G (p = 0.02, OR 2.7) at PANX2. We found no significant association with clinical phenotypes. Our limited findings do not support a major contribution of PANX1-3 to disease risk of schizophrenia according to DSM-5. We cannot confirm an association of the PANX2 loci at chromosome 22q13 with periodic catatonia.

Variants in Adjacent Oxytocin/Vasopressin Gene Region and Associations with ASD Diagnosis and Other Autism Related Endophenotypes

BACKGROUND

There has been increasing interest in oxytocin (peptide: OT, gene: OXT) as a treatment pathway for neurodevelopmental disorders such as Autism Spectrum Disorder (ASD). Neurodevelopmental disorders affect functional, social, and intellectual abilities. With advances in molecular biology, research has connected multiple gene regions to the clinical presentation of ASD. Studies have also shown that the neuropeptide hormones OT and arginine vasopressin (AVP) influence mammalian social and territorial behaviors and may have treatment potential for neurodevelopmental disorders. Published data examining molecular and phenotypic variation in ASD, such as cognitive abilities, are limited. Since most studies have focused on the receptors in the OT-AVP system, we investigated genetic variation within peptide genes for association with phenotypic ASD features that help identify subgroups within the spectrum.

METHODS

In this study, TDT analysis was carried out utilizing FBAT in 207 probands (156 trios) and a European Ancestry (EA) subsample (108 trios).The evolutionarily related and adjacent genes of OXT and AVP were studied for associations between the tagged single nucleotide polymorphisms and ASD diagnosis, social abilities, restrictive and repetitive behaviors, and IQ for cognitive abilities. Additionally, relationships with whole blood serotonin (WB5HT) were explored because of the developmental relationships connecting plasma levels of OT and WB5HT within ASD.

RESULTS

RESULTS indicate significant association between OXT rs6084258 (p = 0.001) and ASD. Associations with several endophenotypes were also noted: OXT rs6133010 was associated with IQ (full scale IQ, p = 0.008; nonverbal IQ, p = 0.010, verbal IQ, p = 0.006); and OXT rs4813625 and OXT rs877172 were associated with WB5HT levels (EA, p = 0.027 and p = 0.033, respectively). Additionally, we measured plasma OT (pOT) levels in a subsample (N = 54). RESULTS show the three polymorphisms, OXT rs6084258, OXT rs11697250, and OXT rs877172, have significant association with pOT (EA, p = 0.011, p = 0.010, and p = 0.002, respectively).

CONCLUSIONS

These findings suggest that SNPs near OXT and AVP are associated with diagnosis of ASD, social behaviors, restricted and repetitive behaviors, IQ, pOT, and WB5HT. Future studies need to replicate these findings and examine gene-interactions in other neurodevelopmental disorders. Mechanisms of action may influence early social and cognitive development that may or may not be limited to ASD diagnosis.

Top associated SNPs in prostate cancer are significantly enriched in cis-expression quantitative trait loci and at transcription factor binding sites

While genome-wide association studies (GWAS) have revealed thousands of disease risk single nucleotide polymorphisms (SNPs), their functions remain largely unknown. Recent studies have suggested the regulatory roles of GWAS risk variants in several common diseases; however, the complex regulatory structure in prostate cancer is unclear.

We investigated the potential regulatory roles of risk variants in two prostate cancer GWAS datasets by their interactions with expression quantitative trait loci (eQTL) and/or transcription factor binding sites (TFBSs) in three populations.

Our results indicated that the moderately associated GWAS SNPs were significantly enriched with cis-eQTLs and TFBSs in Caucasians (CEU), but not in African Americans (AA) or Japanese (JPT); this was also observed in an independent pan-cancer related SNPs from the GWAS Catalog. We found that the eQTL enrichment in the CEU population was tissue-specific to eQTLs from CEU lymphoblastoid cell lines. Importantly, we pinpointed two SNPs, rs2861405 and rs4766642, by overlapping results from cis-eQTL and TFBS as applied to the CEU data.

These results suggested that prostate cancer associated SNPs and pan-cancer associated SNPs are likely to play regulatory roles in CEU. However, the negative enrichment results in AA or JPT and the potential mechanisms remain to be elucidated in additional samples.

Enrichment of inflammatory bowel disease and colorectal cancer risk variants in colon expression quantitative trait loci

Background

Genome-wide association studies (GWAS) have identified single nucleotide polymorphisms (SNPs) associated with diseases of the colon including inflammatory bowel diseases (IBD) and colorectal cancer (CRC). However, the functional role of many of these SNPs is largely unknown and tissue-specific resources are lacking. Expression quantitative trait loci (eQTL) mapping identifies target genes of disease-associated SNPs. This study provides a comprehensive eQTL map of distal colonic samples obtained from 40 healthy African Americans and demonstrates their relevance for GWAS of colonic diseases.

Results

8.4 million imputed SNPs were tested for their associations with 16,252 expression probes representing 12,363 unique genes. 1,941 significant cis-eQTL, corresponding to 122 independent signals, were identified at a false discovery rate (FDR) of 0.01. Overall, among colon cis-eQTL, there was significant enrichment for GWAS variants for IBD (Crohn’s disease [CD] and ulcerative colitis [UC]) and CRC as well as type 2 diabetes and body mass index. ERAP2, ADCY3, INPP5E, UBA7, SFMBT1, NXPE1 and REXO2 were identified as target genes for IBD-associated variants. The CRC-associated eQTL rs3802842 was associated with the expression of C11orf93 (COLCA2). Enrichment of colon eQTL near transcription start sites and for active histone marks was demonstrated, and eQTL with high population differentiation were identified.

Conclusions

Through the comprehensive study of eQTL in the human colon, this study identified novel target genes for IBD- and CRC-associated genetic variants. Moreover, bioinformatic characterization of colon eQTL provides a tissue-specific tool to improve understanding of biological differences in diseases between different ethnic groups.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1292-z) contains supplementary material, which is available to authorized users.

Structural architecture of SNP effects on complex traits

Despite the discovery of copy-number variation (CNV) across the genome nearly 10 years ago, current SNP-based analysis methodologies continue to collapse the homozygous (i.e., A/A), hemizygous (i.e., A/0), and duplicative (i.e., A/A/A) genotype states, treating the genotype variable as irreducible or unaltered by other colocalizing forms of genetic (e.g., structural) variation. Our understanding of common, genome-wide CNVs suggests that the canonical genotype construct might belie the enormous complexity of the genome. Here we present multiple analyses of several phenotypes and provide methods supporting a conceptual shift that embraces the structural dimension of genotype. We comprehensively investigate the impact of the structural dimension of genotype on (1) GWAS methods, (2) interpretation of rare LOF variants, (3) characterization of genomic architecture, and (4) implications for mapping loci involved in complex disease. Taken together, these results argue for the inclusion of a structural dimension and suggest that some portion of the "missing" heritability might be recovered through integration of the structural dimension of SNP effects on complex traits.

Genetic variation associated with euphorigenic effects of d-amphetamine is associated with diminished risk for schizophrenia and attention deficit hyperactivity disorder

Here, we extended our findings from a genome-wide association study of the euphoric response to d-amphetamine in healthy human volunteers by identifying enrichment between SNPs associated with response to d-amphetamine and SNPs associated with psychiatric disorders. We found that SNPs nominally associated (P ≤ 0.05 and P ≤ 0.01) with schizophrenia and attention deficit hyperactivity disorder were also nominally associated with d-amphetamine response. Furthermore, we found that the source of this enrichment was an excess of alleles that increased sensitivity to the euphoric effects of d-amphetamine and decreased susceptibility to schizophrenia and attention deficit hyperactivity disorder. In contrast, three negative control phenotypes (height, inflammatory bowel disease, and Parkinson disease) did not show this enrichment. Taken together, our results suggest that alleles identified using an acute challenge with a dopaminergic drug in healthy individuals can be used to identify alleles that confer risk for psychiatric disorders commonly treated with dopaminergic agonists and antagonists. More importantly, our results show the use of the enrichment approach as an alternative to stringent standards for genome-wide significance and suggest a relatively novel approach to the analysis of small cohorts in which intermediate phenotypes have been measured.

Converging Pathways in Autism Spectrum Disorders: Interplay between Synaptic Dysfunction and Immune Responses

Autism spectrum disorders (ASD) are highly heritable, yet genetically heterogeneous neurodevelopmental conditions. Recent genome-wide association and gene expression studies have provided evidence supporting the notion that the large number of genetic variants associated with ASD converge toward a core set of dysregulated biological processes. Here we review recent data demonstrating the involvement of synaptic dysfunction and abnormal immune responses in ASD, and discuss the functional interplay between the two phenomena.

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.

An eQTL mapping approach reveals that rare variants in the SEMA5A regulatory network impact autism risk

To date, genome-wide single nucleotide polymorphism (SNP) and copy number variant (CNV) association studies of autism spectrum disorders (ASDs) have led to promising signals but not to easily interpretable or translatable results. Our own genome-wide association study (GWAS) showed significant association to an intergenic SNP near Semaphorin 5A (SEMA5A) and provided evidence for reduced expression of the same gene. In a novel GWAS follow-up approach, we map an expression regulatory pathway for a GWAS candidate gene, SEMA5A, in silico by using population expression and genotype data sets. We find that the SEMA5A regulatory network significantly overlaps rare autism-specific CNVs. The SEMA5A regulatory network includes previous autism candidate genes and regions, including MACROD2, A2BP1, MCPH1, MAST4, CDH8, CADM1, FOXP1, AUTS2, MBD5, 7q21, 20p, USH2A, KIRREL3, DBF4B and RELN, among others. Our results provide: (i) a novel data-derived network implicated in autism, (ii) evidence that the same pathway seeded by an initial SNP association shows association with rare genetic variation in ASDs, (iii) a potential mechanism of action and interpretation for the previous autism candidate genes and genetic variants that fall in this network, and (iv) a novel approach that can be applied to other candidate genes for complex genetic disorders. We take a step towards better understanding of the significance of SEMA5A pathways in autism that can guide interpretation of many other genetic results in ASDs.

Individual common variants exert weak effects on the risk for autism spectrum disorders

While it is apparent that rare variation can play an important role in the genetic architecture of autism spectrum disorders (ASDs), the contribution of common variation to the risk of developing ASD is less clear. To produce a more comprehensive picture, we report Stage 2 of the Autism Genome Project genome-wide association study, adding 1301 ASD families and bringing the total to 2705 families analysed (Stages 1 and 2). In addition to evaluating the association of individual single nucleotide polymorphisms (SNPs), we also sought evidence that common variants, en masse, might affect the risk. Despite genotyping over a million SNPs covering the genome, no single SNP shows significant association with ASD or selected phenotypes at a genome-wide level. The SNP that achieves the smallest P-value from secondary analyses is rs1718101. It falls in CNTNAP2, a gene previously implicated in susceptibility for ASD. This SNP also shows modest association with age of word/phrase acquisition in ASD subjects, of interest because features of language development are also associated with other variation in CNTNAP2. In contrast, allele scores derived from the transmission of common alleles to Stage 1 cases significantly predict case status in the independent Stage 2 sample. Despite being significant, the variance explained by these allele scores was small (Vm< 1%). Based on results from individual SNPs and their en masse effect on risk, as inferred from the allele score results, it is reasonable to conclude that common variants affect the risk for ASD but their individual effects are modest.

Common genetic variants, acting additively, are a major source of risk for autism

Background

Autism spectrum disorders (ASD) are early onset neurodevelopmental syndromes typified by impairments in reciprocal social interaction and communication, accompanied by restricted and repetitive behaviors. While rare and especially de novo genetic variation are known to affect liability, whether common genetic polymorphism plays a substantial role is an open question and the relative contribution of genes and environment is contentious. It is probable that the relative contributions of rare and common variation, as well as environment, differs between ASD families having only a single affected individual (simplex) versus multiplex families who have two or more affected individuals.

Methods

By using quantitative genetics techniques and the contrast of ASD subjects to controls, we estimate what portion of liability can be explained by additive genetic effects, known as narrow-sense heritability. We evaluate relatives of ASD subjects using the same methods to evaluate the assumptions of the additive model and partition families by simplex/multiplex status to determine how heritability changes with status.

Results

By analyzing common variation throughout the genome, we show that common genetic polymorphism exerts substantial additive genetic effects on ASD liability and that simplex/multiplex family status has an impact on the identified composition of that risk. As a fraction of the total variation in liability, the estimated narrow-sense heritability exceeds 60% for ASD individuals from multiplex families and is approximately 40% for simplex families. By analyzing parents, unaffected siblings and alleles not transmitted from parents to their affected children, we conclude that the data for simplex ASD families follow the expectation for additive models closely. The data from multiplex families deviate somewhat from an additive model, possibly due to parental assortative mating.

Conclusions

Our results, when viewed in the context of results from genome-wide association studies, demonstrate that a myriad of common variants of very small effect impacts ASD liability.