CWAS-Plus: estimating category-wide association of rare noncoding variation from whole-genome sequencing data with cell-type-specific functional data.

Variants in cis-regulatory elements link the noncoding genome to human pathology; however, detailed analytic tools for understanding the association between cell-level brain pathology and noncoding variants are lacking. CWAS-Plus, adapted from a Python package for category-wide association testing (CWAS), enhances noncoding variant analysis by integrating both whole-genome sequencing (WGS) and user-provided functional data. With simplified parameter settings and an efficient multiple testing correction method, CWAS-Plus conducts the CWAS workflow 50 times faster than CWAS, making it more accessible and user-friendly for researchers. Here, we used a single-nuclei assay for transposase-accessible chromatin with sequencing to facilitate CWAS-guided noncoding variant analysis at cell-type-specific enhancers and promoters. Examining autism spectrum disorder WGS data (n = 7280), CWAS-Plus identified noncoding de novo variant associations in transcription factor binding sites within conserved loci. Independently, in Alzheimer's disease WGS data (n = 1087), CWAS-Plus detected rare noncoding variant associations in microglia-specific regulatory elements. These findings highlight CWAS-Plus's utility in genomic disorders and scalability for processing large-scale WGS data and in multiple-testing corrections. CWAS-Plus and its user manual are available at https://github.com/joonan-lab/cwas/ and https://cwas-plus.readthedocs.io/en/latest/, respectively.

Sex-Differential Gene Expression in Developing Human Cortex and Its Intersection With Autism Risk Pathways.

Background: Sex-differential biology may contribute to the consistently male-biased prevalence of autism spectrum disorder (ASD). Gene expression differences between males and females in the brain can indicate possible molecular and cellular mechanisms involved, although transcriptomic sex differences during human prenatal cortical development have been incompletely characterized, primarily due to small sample sizes. Methods: We performed a meta-analysis of sex-differential expression and co-expression network analysis in 2 independent bulk RNA sequencing datasets generated from cortex of 273 prenatal donors without known neuropsychiatric disorders. To assess the intersection between neurotypical sex differences and neuropsychiatric disorder biology, we tested for enrichment of ASD-associated risk genes and expression changes, neuropsychiatric disorder risk genes, and cell type markers within identified sex-differentially expressed genes (sex-DEGs) and sex-differential co-expression modules. Results: We identified 101 significant sex-DEGs, including Y-chromosome genes, genes impacted by X-chromosome inactivation, and autosomal genes. Known ASD risk genes, implicated by either common or rare variants, did not preferentially overlap with sex-DEGs. We identified 1 male-specific co-expression module enriched for immune signaling that is unique to 1 input dataset. Conclusions: Sex-differential gene expression is limited in prenatal human cortex tissue, although meta-analysis of large datasets allows for the identification of sex-DEGs, including autosomal genes that encode proteins involved in neural development. Lack of sex-DEG overlap with ASD risk genes in the prenatal cortex suggests that sex-differential modulation of ASD symptoms may occur in other brain regions, at other developmental stages, or in specific cell types, or may involve mechanisms that act downstream from mutation-carrying genes.

Males are more commonly diagnosed with autism spectrum disorder than females, and sex differences in brain development may contribute to this difference. Here, we define differences in gene expression patterns between males and females in human prenatal brain tissue from 273 donors to identify 101 genes that are expressed at different levels in males and females and gene sets that show sex-specific expression correlations. Genes with autism-associated DNA variants and genes with altered expression in autism do not preferentially overlap with sex-differential genes, suggesting that sex-differential biology may influence autism risk mechanisms in other brain regions, at other developmental stages, or in specific cell types.

CWAS-Plus: Estimating category-wide association of rare noncoding variation from whole-genome sequencing data with cell-type-specific functional data.

Variants in cis-regulatory elements link the noncoding genome to human brain pathology; however, detailed analytic tools for understanding the association between cell-level brain pathology and noncoding variants are lacking. CWAS-Plus, adapted from a Python package for category-wide association testing (CWAS) employs both whole-genome sequencing and user-provided functional data to enhance noncoding variant analysis, with a faster and more efficient execution of the CWAS workflow. Here, we used single-nuclei assay for transposase-accessible chromatin with sequencing to facilitate CWAS-guided noncoding variant analysis at cell-type specific enhancers and promoters. Examining autism spectrum disorder whole-genome sequencing data (n = 7,280), CWAS-Plus identified noncoding de novo variant associations in transcription factor binding sites within conserved loci. Independently, in Alzheimer's disease whole-genome sequencing data (n = 1,087), CWAS-Plus detected rare noncoding variant associations in microglia-specific regulatory elements. These findings highlight CWAS-Plus's utility in genomic disorders and scalability for processing large-scale whole-genome sequencing data and in multiple-testing corrections. CWAS-Plus and its user manual are available at https://github.com/joonan-lab/cwas/ and https://cwas-plus.readthedocs.io/en/latest/, respectively.

Sex differences in autism spectrum disorders.

PURPOSE OF REVIEW: A strong male bias in autism spectrum disorder (ASD) prevalence has been observed with striking consistency, but no mechanism has yet to definitively account for this sex difference. This review explores the current status of epidemiological, genetic, and neuroendocrinological work addressing ASD prevalence and liability in males and females, so as to frame the major issues necessary to pursue a more complete understanding of the biological basis for sex-differential risk. RECENT FINDINGS: Recent studies continue to report a male bias in ASD prevalence, but also suggest that sex differences in phenotypic presentation, including fewer restricted and repetitive behaviors and externalizing behavioral problems in females, may contribute to this bias. Genetic studies demonstrate that females are protected from the effects of heritable and de-novo ASD risk variants, and compelling work suggests that sex chromosomal genes and/or sex hormones, especially testosterone, may modulate the effects of genetic variation on the presentation of an autistic phenotype. SUMMARY: ASDs affect females less frequently than males, and several sex-differential genetic and hormonal factors may contribute. Future work to determine the mechanisms by which these factors confer risk and protection to males and females is essential.

Species-specific FMRP regulation of RACK1 is critical for prenatal cortical development.

Fragile X messenger ribonucleoprotein 1 protein (FMRP) deficiency leads to fragile X syndrome (FXS), an autism spectrum disorder. The role of FMRP in prenatal human brain development remains unclear. Here, we show that FMRP is important for human and macaque prenatal brain development. Both FMRP-deficient neurons in human fetal cortical slices and FXS patient stem cell-derived neurons exhibit mitochondrial dysfunctions and hyperexcitability. Using multiomics analyses, we have identified both FMRP-bound mRNAs and FMRP-interacting proteins in human neurons and unveiled a previously unknown role of FMRP in regulating essential genes during human prenatal development. We demonstrate that FMRP interaction with CNOT1 maintains the levels of receptor for activated C kinase 1 (RACK1), a species-specific FMRP target. Genetic reduction of RACK1 leads to both mitochondrial dysfunctions and hyperexcitability, resembling FXS neurons. Finally, enhancing mitochondrial functions rescues deficits of FMRP-deficient cortical neurons during prenatal development, demonstrating targeting mitochondrial dysfunction as a potential treatment.

Neural Transcriptomic Analysis of Sex Differences in Autism Spectrum Disorder: Current Insights and Future Directions.

Autism spectrum disorder (ASD) is consistently diagnosed 3 to 5 times more frequently in males than females, a dramatically sex-biased prevalence that suggests the involvement of sex-differential biological factors in modulating risk. The genomic scale of transcriptomic analyses of human brain tissue can provide an unbiased approach for identifying genes and associated functional processes at the intersection of sex-differential and ASD-impacted neurobiology. Several studies characterizing gene expression changes in the ASD brain have been published in recent years with increasing sample size and cellular resolution. These studies report several convergent patterns across data sets and genetically heterogeneous samples in the ASD brain, including elevated expression of gene sets associated with glial and immune function, and reduced expression of gene sets associated with neuronal and synaptic functions. Assessment of neurotypical cortex tissue has reported parallel patterns by sex, with male-elevated expression of overlapping sets of glial/immune-related genes and female-biased expression of neuron-associated genes, suggesting potential roles for these cell types in sex-differential ASD risk mechanisms. However, validating and further exploring these mechanisms is challenged by the available data, as existing studies of ASD brain include a limited number of female ASD donors and focus predominantly on cortex regions not known to show pronounced sex-differential morphology or function. With this review, we summarize convergent findings from several landmark studies of the transcriptome in ASD brain and their relationship to sex-differential gene expression, and we discuss limitations and remaining questions regarding transcriptomic analysis of sex differences in ASD.

SUPERGNOVA: local genetic correlation analysis reveals heterogeneous etiologic sharing of complex traits.

Local genetic correlation quantifies the genetic similarity of complex traits in specific genomic regions. However, accurate estimation of local genetic correlation remains challenging, due to linkage disequilibrium in local genomic regions and sample overlap across studies. We introduce SUPERGNOVA, a statistical framework to estimate local genetic correlations using summary statistics from genome-wide association studies. We demonstrate that SUPERGNOVA outperforms existing methods through simulations and analyses of 30 complex traits. In particular, we show that the positive yet paradoxical genetic correlation between autism spectrum disorder and cognitive performance could be explained by two etiologically distinct genetic signatures with bidirectional local genetic correlations.

Developmental dynamics of voltage-gated sodium channel isoform expression in the human and mouse brain.

BACKGROUND: Genetic variants in the voltage-gated sodium channels SCN1A, SCN2A, SCN3A, and SCN8A are leading causes of epilepsy, developmental delay, and autism spectrum disorder. The mRNA splicing patterns of all four genes vary across development in the rodent brain, including mutually exclusive copies of the fifth protein-coding exon detected in the neonate (5N) and adult (5A). A second pair of mutually exclusive exons is reported in SCN8A only (18N and 18A). We aimed to quantify the expression of individual exons in the developing human brain. METHODS: RNA-seq data from 783 human brain samples across development were analyzed to estimate exon-level expression. Developmental changes in exon utilization were validated by assessing intron splicing. Exon expression was also estimated in RNA-seq data from 58 developing mouse neocortical samples. RESULTS: In the mature human neocortex, exon 5A is consistently expressed at least 4-fold higher than exon 5N in all four genes. For SCN2A, SCN3A, and SCN8A, a brain-wide synchronized 5N to 5A transition occurs between 24 post-conceptual weeks (2nd trimester) and 6 years of age. In mice, the equivalent 5N to 5A transition begins at or before embryonic day 15.5. In SCN8A, over 90% of transcripts in the mature human cortex include exon 18A. Early in fetal development, most transcripts include 18N or skip both 18N and 18A, with a transition to 18A inclusion occurring from 13 post-conceptual weeks to 6 months of age. No other protein-coding exons showed comparably dynamic developmental trajectories. CONCLUSIONS: Exon usage in SCN1A, SCN2A, SCN3A, and SCN8A changes dramatically during human brain development. These splice isoforms, which alter the biophysical properties of the encoded channels, may account for some of the observed phenotypic differences across development and between specific variants. Manipulation of the proportion of splicing isoforms at appropriate stages of development may act as a therapeutic strategy for specific mutations or even epilepsy in general.

Whole-Genome and RNA Sequencing Reveal Variation and Transcriptomic Coordination in the Developing Human Prefrontal Cortex.

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.

Publisher Correction: Whole genome sequencing in psychiatric disorders: the WGSPD consortium.

In the version of this article initially published, the consortium authorship and corresponding authors were not presented correctly. In the PDF and print versions, the Whole Genome Sequencing for Psychiatric Disorders (WGSPD) consortium was missing from the author list at the beginning of the paper, where it should have appeared as the seventh author; it was present in the author list at the end of the paper, but the footnote directing readers to the Supplementary Note for a list of members was missing. In the HTML version, the consortium was listed as the last author instead of as the seventh, and the line directing readers to the Supplementary Note for a list of members appeared at the end of the paper under Author Information but not in association with the consortium name itself. Also, this line stated that both member names and affiliations could be found in the Supplementary Note; in fact, only names are given. In all versions of the paper, the corresponding author symbols were attached to A. Jeremy Willsey, Steven E. Hyman, Anjene M. Addington and Thomas Lehner; they should have been attached, respectively, to Steven E. Hyman, Anjene M. Addington, Thomas Lehner and Nelson B. Freimer. As a result of this shift, the respective contact links in the HTML version did not lead to the indicated individuals. The errors have been corrected in the HTML and PDF versions of the article.

Genome-wide de novo risk score implicates promoter variation in autism spectrum disorder.

Whole-genome sequencing (WGS) has facilitated the first genome-wide evaluations of the contribution of de novo noncoding mutations to complex disorders. Using WGS, we identified 255,106 de novo mutations among sample genomes from members of 1902 quartet families in which one child, but not a sibling or their parents, was affected by autism spectrum disorder (ASD). In contrast to coding mutations, no noncoding functional annotation category, analyzed in isolation, was significantly associated with ASD. Casting noncoding variation in the context of a de novo risk score across multiple annotation categories, however, did demonstrate association with mutations localized to promoter regions. We found that the strongest driver of this promoter signal emanates from evolutionarily conserved transcription factor binding sites distal to the transcription start site. These data suggest that de novo mutations in promoter regions, characterized by evolutionary and functional signatures, contribute to ASD.

Integrative functional genomic analysis of human brain development and neuropsychiatric risks.

To broaden our understanding of human neurodevelopment, we profiled transcriptomic and epigenomic landscapes across brain regions and/or cell types for the entire span of prenatal and postnatal development. Integrative analysis revealed temporal, regional, sex, and cell type-specific dynamics. We observed a global transcriptomic cup-shaped pattern, characterized by a late fetal transition associated with sharply decreased regional differences and changes in cellular composition and maturation, followed by a reversal in childhood-adolescence, and accompanied by epigenomic reorganizations. Analysis of gene coexpression modules revealed relationships with epigenomic regulation and neurodevelopmental processes. Genes with genetic associations to brain-based traits and neuropsychiatric disorders (including MEF2C, SATB2, SOX5, TCF4, and TSHZ3) converged in a small number of modules and distinct cell types, revealing insights into neurodevelopment and the genomic basis of neuropsychiatric risks.

An analytical framework for whole-genome sequence association studies and its implications for autism spectrum disorder.

Genomic association studies of common or rare protein-coding variation have established robust statistical approaches to account for multiple testing. Here we present a comparable framework to evaluate rare and de novo noncoding single-nucleotide variants, insertion/deletions, and all classes of structural variation from whole-genome sequencing (WGS). Integrating genomic annotations at the level of nucleotides, genes, and regulatory regions, we define 51,801 annotation categories. Analyses of 519 autism spectrum disorder families did not identify association with any categories after correction for 4,123 effective tests. Without appropriate correction, biologically plausible associations are observed in both cases and controls. Despite excluding previously identified gene-disrupting mutations, coding regions still exhibited the strongest associations. Thus, in autism, the contribution of de novo noncoding variation is probably modest in comparison to that of de novo coding variants. Robust results from future WGS studies will require large cohorts and comprehensive analytical strategies that consider the substantial multiple-testing burden.

Identification of Developmental and Behavioral Markers Associated With Genetic Abnormalities in Autism Spectrum Disorder.

OBJECTIVE: Aside from features associated with risk of neurogenetic syndromes in general (e.g., cognitive impairment), limited progress has been made in identifying phenotype-genotype relationships in autism spectrum disorder (ASD). The objective of this study was to extend work in the Simons Simplex Collection by comparing the phenotypic profiles of ASD probands with or without identified de novo loss of function mutations or copy number variants in high-confidence ASD-associated genes or loci. METHOD: Analyses preemptively accounted for documented differences in sex and IQ in affected individuals with de novo mutations by matching probands with and without these genetic events on sex, IQ, and age before comparing them on multiple behavioral domains. RESULTS: Children with de novo mutations (N=112) had a greater likelihood of motor delay during early development (later age at walking), but they were less impaired on certain measures of ASD core symptoms (parent-rated social communication abnormalities and clinician-rated diagnostic certainty about ASD) in later childhood. These children also showed relative strengths in verbal and language abilities, including a smaller discrepancy between nonverbal and verbal IQ and a greater likelihood of having achieved fluent language (i.e., regular use of complex sentences). CONCLUSIONS: Children with ASD with de novo mutations may exhibit a "muted" symptom profile with respect to social communication and language deficits relative to those with ASD with no identified genetic abnormalities. Such findings suggest that examining early milestone differences and standardized testing results may be helpful in etiologic efforts, and potentially in clinical differentiation of various subtypes of ASD, but only if developmental and demographic variables are properly accounted for first.

The role of sex-differential biology in risk for autism spectrum disorder.

Autism spectrum disorder (ASD) is a developmental condition that affects approximately four times as many males as females, a strong sex bias that has not yet been fully explained. Understanding the causes of this biased prevalence may highlight novel avenues for treatment development that could benefit patients with diverse genetic backgrounds, and the expertise of sex differences researchers will be invaluable in this endeavor. In this review, I aim to assess current evidence pertaining to the sex difference in ASD prevalence and to identify outstanding questions and remaining gaps in our understanding of how males come to be more frequently affected and/or diagnosed with ASD. Though males consistently outnumber females in ASD prevalence studies, prevalence estimates generated using different approaches report male/female ratios of variable magnitude that suggest that ascertainment or diagnostic biases may contribute to the male skew in ASD. Here, I present the different methods applied and implications of their findings. Additionally, even as prevalence estimations challenge the degree of male bias in ASD, support is growing for the long-proposed female protective effect model of ASD risk, and I review the relevant results from recurrence rate, quantitative trait, and genetic analyses. Lastly, I describe work investigating several sex-differential biological factors and pathways that may be responsible for females' protection and/or males' increased risk predicted by the female protective effect model, including sex steroid hormone exposure and regulation and sex-differential activity of certain neural cell types. However, much future work from both the ASD and sex differences research communities will be required to flesh out our understanding of how these factors act to influence the developing brain and modulate ASD risk.

Gene expression in human brain implicates sexually dimorphic pathways in autism spectrum disorders.

Autism spectrum disorder (ASD) is more prevalent in males, and the mechanisms behind this sex-differential risk are not fully understood. Two competing, but not mutually exclusive, hypotheses are that ASD risk genes are sex-differentially regulated, or alternatively, that they interact with characteristic sexually dimorphic pathways. Here we characterized sexually dimorphic gene expression in multiple data sets from neurotypical adult and prenatal human neocortical tissue, and evaluated ASD risk genes for evidence of sex-biased expression. We find no evidence for systematic sex-differential expression of ASD risk genes. Instead, we observe that genes expressed at higher levels in males are significantly enriched for genes upregulated in post-mortem autistic brain, including astrocyte and microglia markers. This suggests that it is not sex-differential regulation of ASD risk genes, but rather naturally occurring sexually dimorphic processes, potentially including neuron-glial interactions, that modulate the impact of risk variants and contribute to the sex-skewed prevalence of ASD.

Recurrence rates provide evidence for sex-differential, familial genetic liability for autism spectrum disorders in multiplex families and twins.

BACKGROUND: Autism spectrum disorders (ASDs) are more prevalent in males, suggesting a multiple threshold liability model in which females are, on average, protected by sex-differential mechanisms. Under this model, autistic females are predicted to carry a more penetrant risk variant load than males and to share this greater genetic liability with their siblings. However, reported ASD recurrence rates have not demonstrated significantly increased risk to siblings of affected girls. Here, we characterize recurrence patterns in multiplex families from the Autism Genetics Resource Exchange (AGRE) to determine if risk in these families follows a female protective model. METHODS: We assess recurrence rates and quantitative traits in full siblings from 1,120 multiplex nuclear families and concordance rates in 305 twin pairs from AGRE. We consider the first two affected children per family, and one randomly selected autistic twin per pair, as probands. We then compare recurrence rates and phenotypes between males and females and between twin pairs or families with at least one female proband (female-containing (FC)) versus those with only male probands (male-only (MO)). RESULTS: Among children born after two probands, we observe significantly higher recurrence in males (47.5%) than in females (21.1%; relative risk, RR = 2.25; adjusted P = 6.22e-08) and in siblings of female (44.3%) versus siblings of male probands (30.4%; RR = 1.46; adj. P = 0.036). This sex-differential recurrence is also robust in dizygotic twin pairs (males = 61.5%, females = 19.1%; RR = 3.23; adj. P = 7.66e-09). Additionally, we find a significant negative relationship between interbirth interval and ASD recurrence that is driven by children in MO families. CONCLUSIONS: By classifying families as MO or FC using two probands instead of one, we observe significant recurrence rate differences between families harboring sex-differential familial liability. However, a significant sex difference in risk to children within FC families suggests that female protective mechanisms are still operative in families carrying high genetic risk loads. Furthermore, the male-specific relationship between shorter interbirth intervals and increased ASD risk is consistent with a potentially greater contribution from environmental factors in males versus higher genetic risk in affected females and their families. Understanding the mechanisms driving these sex-differential risk profiles will be useful for treatment development and prevention.

Social responsiveness, an autism endophenotype: genomewide significant linkage to two regions on chromosome 8.

OBJECTIVE: Autism spectrum disorder is characterized by deficits in social function and the presence of repetitive and restrictive behaviors. Following a previous test of principle, the authors adopted a quantitative approach to discovering genes contributing to the broader autism phenotype by using social responsiveness as an endophenotype for autism spectrum disorder. METHOD: Linkage analyses using scores from the Social Responsiveness Scale were performed in 590 families from the Autism Genetic Resource Exchange, a largely multiplex autism spectrum disorder cohort. Regional and genomewide association analyses were performed to search for common variants contributing to social responsiveness. RESULTS: Social Responsiveness Scale scores were unimodally distributed in male offspring from multiplex autism families, in contrast with a bimodal distribution observed in female offspring. In correlated analyses differing by Social Responsiveness Scale respondent, genomewide significant linkage for social responsiveness was identified at chr8p21.3 (multipoint LOD=4.11; teacher/parent scores) and chr8q24.22 (multipoint LOD=4.54; parent-only scores), respectively. Genomewide or linkage-directed association analyses did not detect common variants contributing to social responsiveness. CONCLUSIONS: The sex-differential distributions of Social Responsiveness Scale scores in multiplex autism families likely reflect mechanisms contributing to the sex ratio for autism observed in the general population and form a quantitative signature of reduced penetrance of inherited liability to autism spectrum disorder among females. The identification of two strong loci for social responsiveness validates the endophenotype approach for the identification of genetic variants contributing to complex traits such as autism spectrum disorder. While causal mutations have yet to be identified, these findings are consistent with segregation of rare genetic variants influencing social responsiveness and underscore the increasingly recognized role of rare inherited variants in the genetic architecture of autism spectrum disorder.

Insights into Autism Spectrum Disorder Genomic Architecture and Biology from 71 Risk Loci.

Analysis of de novo CNVs (dnCNVs) from the full Simons Simplex Collection (SSC) (N = 2,591 families) replicates prior findings of strong association with autism spectrum disorders (ASDs) and confirms six risk loci (1q21.1, 3q29, 7q11.23, 16p11.2, 15q11.2-13, and 22q11.2). The addition of published CNV data from the Autism Genome Project (AGP) and exome sequencing data from the SSC and the Autism Sequencing Consortium (ASC) shows that genes within small de novo deletions, but not within large dnCNVs, significantly overlap the high-effect risk genes identified by sequencing. Alternatively, large dnCNVs are found likely to contain multiple modest-effect risk genes. Overall, we find strong evidence that de novo mutations are associated with ASD apart from the risk for intellectual disability. Extending the transmission and de novo association test (TADA) to include small de novo deletions reveals 71 ASD risk loci, including 6 CNV regions (noted above) and 65 risk genes (FDR ≤ 0.1).