Inherited and De Novo Genetic Risk for Autism Impacts Shared Networks.

We performed a comprehensive assessment of rare inherited variation in autism spectrum disorder (ASD) by analyzing whole-genome sequences of 2,308 individuals from families with multiple affected children. We implicate 69 genes in ASD risk, including 24 passing genome-wide Bonferroni correction and 16 new ASD risk genes, most supported by rare inherited variants, a substantial extension of previous findings. Biological pathways enriched for genes harboring inherited variants represent cytoskeletal organization and ion transport, which are distinct from pathways implicated in previous studies. Nevertheless, the de novo and inherited genes contribute to a common protein-protein interaction network. We also identified structural variants (SVs) affecting non-coding regions, implicating recurrent deletions in the promoters of DLG2 and NR3C2. Loss of nr3c2 function in zebrafish disrupts sleep and social function, overlapping with human ASD-related phenotypes. These data support the utility of studying multiplex families in ASD and are available through the Hartwell Autism Research and Technology portal.

Integrative functional genomic analyses implicate specific molecular pathways and circuits in autism.

Genetic studies have identified dozens of autism spectrum disorder (ASD) susceptibility genes, raising two critical questions: (1) do these genetic loci converge on specific biological processes, and (2) where does the phenotypic specificity of ASD arise, given its genetic overlap with intellectual disability (ID)? To address this, we mapped ASD and ID risk genes onto coexpression networks representing developmental trajectories and transcriptional profiles representing fetal and adult cortical laminae. ASD genes tightly coalesce in modules that implicate distinct biological functions during human cortical development, including early transcriptional regulation and synaptic development. Bioinformatic analyses suggest that translational regulation by FMRP and transcriptional coregulation by common transcription factors connect these processes. At a circuit level, ASD genes are enriched in superficial cortical layers and glutamatergic projection neurons. Furthermore, we show that the patterns of ASD and ID risk genes are distinct, providing a biological framework for further investigating the pathophysiology of ASD.

The contributions of rare inherited and polygenic risk to ASD in multiplex families.

Autism spectrum disorder (ASD) has a complex genetic architecture involving contributions from both de novo and inherited variation. Few studies have been designed to address the role of rare inherited variation or its interaction with common polygenic risk in ASD. Here, we performed whole-genome sequencing of the largest cohort of multiplex families to date, consisting of 4,551 individuals in 1,004 families having two or more autistic children. Using this study design, we identify seven previously unrecognized ASD risk genes supported by a majority of rare inherited variants, finding support for a total of 74 genes in our cohort and a total of 152 genes after combined analysis with other studies. Autistic children from multiplex families demonstrate an increased burden of rare inherited protein-truncating variants in known ASD risk genes. We also find that ASD polygenic score (PGS) is overtransmitted from nonautistic parents to autistic children who also harbor rare inherited variants, consistent with combinatorial effects in the offspring, which may explain the reduced penetrance of these rare variants in parents. We also observe that in addition to social dysfunction, language delay is associated with ASD PGS overtransmission. These results are consistent with an additive complex genetic risk architecture of ASD involving rare and common variation and further suggest that language delay is a core biological feature of ASD.

Author Correction: Genome-wide changes in lncRNA, splicing, and regional gene expression patterns in autism.

Change history: In this Letter, the labels for splicing events A3SS and A5SS were swapped in column D of Supplementary Table 3a and b. This has been corrected online.

Impact of autism genetic risk on brain connectivity: a mechanism for the female protective effect.

The biological mechanisms underlying the greater prevalence of autism spectrum disorder in males than females remain poorly understood. One hypothesis posits that this female protective effect arises from genetic load for autism spectrum disorder differentially impacting male and female brains. To test this hypothesis, we investigated the impact of cumulative genetic risk for autism spectrum disorder on functional brain connectivity in a balanced sample of boys and girls with autism spectrum disorder and typically developing boys and girls (127 youth, ages 8-17). Brain connectivity analyses focused on the salience network, a core intrinsic functional connectivity network which has previously been implicated in autism spectrum disorder. The effects of polygenic risk on salience network functional connectivity were significantly modulated by participant sex, with genetic load for autism spectrum disorder influencing functional connectivity in boys with and without autism spectrum disorder but not girls. These findings support the hypothesis that autism spectrum disorder risk genes interact with sex differential processes, thereby contributing to the male bias in autism prevalence and proposing an underlying neurobiological mechanism for the female protective effect.

Genome-wide changes in lncRNA, splicing, and regional gene expression patterns in autism.

Autism spectrum disorder (ASD) involves substantial genetic contributions. These contributions are profoundly heterogeneous but may converge on common pathways that are not yet well understood. Here, through post-mortem genome-wide transcriptome analysis of the largest cohort of samples analysed so far, to our knowledge, we interrogate the noncoding transcriptome, alternative splicing, and upstream molecular regulators to broaden our understanding of molecular convergence in ASD. Our analysis reveals ASD-associated dysregulation of primate-specific long noncoding RNAs (lncRNAs), downregulation of the alternative splicing of activity-dependent neuron-specific exons, and attenuation of normal differences in gene expression between the frontal and temporal lobes. Our data suggest that SOX5, a transcription factor involved in neuron fate specification, contributes to this reduction in regional differences. We further demonstrate that a genetically defined subtype of ASD, chromosome 15q11.2-13.1 duplication syndrome (dup15q), shares the core transcriptomic signature observed in idiopathic ASD. Co-expression network analysis reveals that individuals with ASD show age-related changes in the trajectory of microglial and synaptic function over the first two decades, and suggests that genetic risk for ASD may influence changes in regional cortical gene expression. Our findings illustrate how diverse genetic perturbations can lead to phenotypic convergence at multiple biological levels in a complex neuropsychiatric disorder.

Rare Inherited and De Novo CNVs Reveal Complex Contributions to ASD Risk in Multiplex Families.

Rare mutations, including copy-number variants (CNVs), contribute significantly to autism spectrum disorder (ASD) risk. Although their importance has been established in families with only one affected child (simplex families), the contribution of both de novo and inherited CNVs to ASD in families with multiple affected individuals (multiplex families) is less well understood. We analyzed 1,532 families from the Autism Genetic Resource Exchange (AGRE) to assess the impact of de novo and rare CNVs on ASD risk in multiplex families. We observed a higher burden of large, rare CNVs, including inherited events, in individuals with ASD than in their unaffected siblings (odds ratio [OR] = 1.7), but the rate of de novo events was significantly lower than in simplex families. In previously characterized ASD risk loci, we identified 49 CNVs, comprising 24 inherited events, 19 de novo events, and 6 events of unknown inheritance, a significant enrichment in affected versus control individuals (OR = 3.3). In 21 of the 30 families (71%) in whom at least one affected sibling harbored an established ASD major risk CNV, including five families harboring inherited CNVs, the CNV was not shared by all affected siblings, indicating that other risk factors are contributing. We also identified a rare risk locus for ASD and language delay at chromosomal region 2q24 (implicating NR4A2) and another lower-penetrance locus involving inherited deletions and duplications of WWOX. The genetic architecture in multiplex families differs from that in simplex families and is complex, warranting more complete genetic characterization of larger multiplex ASD cohorts.

Genome-wide survey of copy number variants finds MAPT duplications in progressive supranuclear palsy.

BACKGROUND: Progressive supranuclear palsy is a neurodegenerative tauopathy manifesting clinically as a progressive akinetic-rigid syndrome. In this study, we sought to identify genetic variants influencing PSP susceptibility through a genome-wide association analysis of a cohort of well-characterized patients who had participated in the Neuroprotection and Natural History in Parkinson Plus Syndromes and Blood Brain Barrier in Parkinson Plus Syndromes studies. METHODS: We genotyped single-nucleotide polymorphisms in 283 PSP cases from the United Kingdom, Germany, and France and compared these with genotypes from 4472 controls. Copy number variants were identified from genotyping data. RESULTS: We observed associations on chromosome 17 within or close to the MAPT gene and explored the genetic architecture at this locus. We confirmed the previously reported association of rs1768208 in the MOBP gene (P = 3.29 × 10-13 ) and rs1411478 in STX6 (P = 3.45 × 10-10 ). The population-attributable risk from the MAPT, MOBP, and STX6 single-nucleotide polymorphisms was found to be 0.37, 0.26, and 0.08, respectively. In addition, we found 2 instances of copy number variants spanning the MAPT gene in patients with PSP. These copy number variants include tau but few other genes within the chromosome 17 haplotype region, providing additional support for the direct pathogenicity of MAPT in PSP. CONCLUSIONS: Clinicians should also be aware of MAPT duplication as a possible genetic cause of PSP, especially in patients presenting with young age at onset. © 2019 International Parkinson and Movement Disorder Society.

Evidence for α-synuclein prions causing multiple system atrophy in humans with parkinsonism.

Prions are proteins that adopt alternative conformations that become self-propagating; the PrP(Sc) prion causes the rare human disorder Creutzfeldt-Jakob disease (CJD). We report here that multiple system atrophy (MSA) is caused by a different human prion composed of the α-synuclein protein. MSA is a slowly evolving disorder characterized by progressive loss of autonomic nervous system function and often signs of parkinsonism; the neuropathological hallmark of MSA is glial cytoplasmic inclusions consisting of filaments of α-synuclein. To determine whether human α-synuclein forms prions, we examined 14 human brain homogenates for transmission to cultured human embryonic kidney (HEK) cells expressing full-length, mutant human α-synuclein fused to yellow fluorescent protein (α-syn140*A53T-YFP) and TgM83(+/-) mice expressing α-synuclein (A53T). The TgM83(+/-) mice that were hemizygous for the mutant transgene did not develop spontaneous illness; in contrast, the TgM83(+/+) mice that were homozygous developed neurological dysfunction. Brain extracts from 14 MSA cases all transmitted neurodegeneration to TgM83(+/-) mice after incubation periods of ∼120 d, which was accompanied by deposition of α-synuclein within neuronal cell bodies and axons. All of the MSA extracts also induced aggregation of α-syn*A53T-YFP in cultured cells, whereas none of six Parkinson's disease (PD) extracts or a control sample did so. Our findings argue that MSA is caused by a unique strain of α-synuclein prions, which is different from the putative prions causing PD and from those causing spontaneous neurodegeneration in TgM83(+/+) mice. Remarkably, α-synuclein is the first new human prion to be identified, to our knowledge, since the discovery a half century ago that CJD was transmissible.

Convergence of genes and cellular pathways dysregulated in autism spectrum disorders.

Rare copy-number variation (CNV) is an important source of risk for autism spectrum disorders (ASDs). We analyzed 2,446 ASD-affected families and confirmed an excess of genic deletions and duplications in affected versus control groups (1.41-fold, p = 1.0 × 10(-5)) and an increase in affected subjects carrying exonic pathogenic CNVs overlapping known loci associated with dominant or X-linked ASD and intellectual disability (odds ratio = 12.62, p = 2.7 × 10(-15), ∼3% of ASD subjects). Pathogenic CNVs, often showing variable expressivity, included rare de novo and inherited events at 36 loci, implicating ASD-associated genes (CHD2, HDAC4, and GDI1) previously linked to other neurodevelopmental disorders, as well as other genes such as SETD5, MIR137, and HDAC9. Consistent with hypothesized gender-specific modulators, females with ASD were more likely to have highly penetrant CNVs (p = 0.017) and were also overrepresented among subjects with fragile X syndrome protein targets (p = 0.02). Genes affected by de novo CNVs and/or loss-of-function single-nucleotide variants converged on networks related to neuronal signaling and development, synapse function, and chromatin regulation.

Transcriptomic analysis of autistic brain reveals convergent molecular pathology.

Autism spectrum disorder (ASD) is a common, highly heritable neurodevelopmental condition characterized by marked genetic heterogeneity. Thus, a fundamental question is whether autism represents an aetiologically heterogeneous disorder in which the myriad genetic or environmental risk factors perturb common underlying molecular pathways in the brain. Here, we demonstrate consistent differences in transcriptome organization between autistic and normal brain by gene co-expression network analysis. Remarkably, regional patterns of gene expression that typically distinguish frontal and temporal cortex are significantly attenuated in the ASD brain, suggesting abnormalities in cortical patterning. We further identify discrete modules of co-expressed genes associated with autism: a neuronal module enriched for known autism susceptibility genes, including the neuronal specific splicing factor A2BP1 (also known as FOX1), and a module enriched for immune genes and glial markers. Using high-throughput RNA sequencing we demonstrate dysregulated splicing of A2BP1-dependent alternative exons in the ASD brain. Moreover, using a published autism genome-wide association study (GWAS) data set, we show that the neuronal module is enriched for genetically associated variants, providing independent support for the causal involvement of these genes in autism. In contrast, the immune-glial module showed no enrichment for autism GWAS signals, indicating a non-genetic aetiology for this process. Collectively, our results provide strong evidence for convergent molecular abnormalities in ASD, and implicate transcriptional and splicing dysregulation as underlying mechanisms of neuronal dysfunction in this disorder.

Genome-wide transcriptome profiling reveals the functional impact of rare de novo and recurrent CNVs in autism spectrum disorders.

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.

Evaluating potential for whole-genome studies in Kosrae, an isolated population in Micronesia.

Whole-genome association studies are predicted to be especially powerful in isolated populations owing to increased linkage disequilibrium (LD) and decreased allelic diversity, but this possibility has not been empirically tested. We compared genome-wide data on 113,240 SNPs typed on 30 trios from the Pacific island of Kosrae to the same markers typed in the 270 samples from the International HapMap Project. The extent of LD is longer and haplotype diversity is lower in Kosrae than in the HapMap populations. More than 98% of Kosraen haplotypes are present in HapMap populations, indicating that HapMap will be useful for genetic studies on Kosrae. The long-range LD around common alleles and limited diversity result in improved efficiency in genetic studies in this population and augments the power to detect association of 'hidden SNPs'.

DASH: a method for identical-by-descent haplotype mapping uncovers association with recent variation.

Rare variants affecting phenotype pose a unique challenge for human genetics. Although genome-wide association studies have successfully detected many common causal variants, they are underpowered in identifying disease variants that are too rare or population-specific to be imputed from a general reference panel and thus are poorly represented on commercial SNP arrays. We set out to overcome these challenges and detect association between disease and rare alleles using SNP arrays by relying on long stretches of genomic sharing that are identical by descent. We have developed an algorithm, DASH, which builds upon pairwise identical-by-descent shared segments to infer clusters of individuals likely to be sharing a single haplotype. DASH constructs a graph with nodes representing individuals and links on the basis of such segments spanning a locus and uses an iterative minimum cut algorithm to identify densely connected components. We have applied DASH to simulated data and diverse GWAS data sets by constructing haplotype clusters and testing them for association. In simulations we show this approach to be significantly more powerful than single-marker testing in an isolated population that is from Kosrae, Federated States of Micronesia and has abundant IBD, and we provide orthogonal information for rare, recent variants in the outbred Wellcome Trust Case-Control Consortium (WTCCC) data. In both cohorts, we identified a number of haplotype associations, five such loci in the WTCCC data and ten in the isolated, that were conditionally significant beyond any individual nearby markers. We have replicated one of these loci in an independent European cohort and identified putative structural changes in low-pass whole-genome sequence of the cluster carriers.

Increased power of mixed models facilitates association mapping of 10 loci for metabolic traits in an isolated population.

The potential benefits of using population isolates in genetic mapping, such as reduced genetic, phenotypic and environmental heterogeneity, are offset by the challenges posed by the large amounts of direct and cryptic relatedness in these populations confounding basic assumptions of independence. We have evaluated four representative specialized methods for association testing in the presence of relatedness; (i) within-family (ii) within- and between-family and (iii) mixed-models methods, using simulated traits for 2906 subjects with known genome-wide genotype data from an extremely isolated population, the Island of Kosrae, Federated States of Micronesia. We report that mixed models optimally extract association information from such samples, demonstrating 88% power to rank the true variant as among the top 10 genome-wide with 56% achieving genome-wide significance, a >80% improvement over the other methods, and demonstrate that population isolates have similar power to non-isolate populations for observing variants of known effects. We then used the mixed-model method to reanalyze data for 17 published phenotypes relating to metabolic traits and electrocardiographic measures, along with another 8 previously unreported. We replicate nine genome-wide significant associations with known loci of plasma cholesterol, high-density lipoprotein, low-density lipoprotein, triglycerides, thyroid stimulating hormone, homocysteine, C-reactive protein and uric acid, with only one detected in the previous analysis of the same traits. Further, we leveraged shared identity-by-descent genetic segments in the region of the uric acid locus to fine-map the signal, refining the known locus by a factor of 4. Finally, we report a novel associations for height (rs17629022, P< 2.1 × 10⁻8).

Prospects for Leveling the Playing Field for Black Children With Autism.

Among the many race-based health disparities that have persistently plagued the US population,1 the disproportionate burden of adverse neurodevelopmental outcomes to Black children affected by autism spectrum disorder (ASD) is particularly devastating given its major lifelong consequences. Recently, in 3 successive reports from the Autism and Developmental Disabilities Monitoring (ADDM) program of the US Centers for Disease Control and Prevention (CDC) (birth cohort years 2014, 2016, and 2018), we and our collaborators reported that although the prevalence of community-diagnosed ASD had equalized for Black and non-Hispanic White (NHW) children in the United States, there has persisted a pronounced racial disparity in the proportion of ASD-affected children with comorbid intellectual disability (ID), on the order of 50% for Black children with ASD vs 20% for White children with ASD.2 Here, we provide data to support the following: much earlier diagnosis is possible; early diagnosis alone is not likely to close the ID comorbidity disparity; and judicious efforts over care as usual are necessary to ensure that Black children have access to timely implementation of developmental therapy, for which we observed promising associations with improved cognitive and adaptive outcomes in our sample.

Systematic haplotype analysis resolves a complex plasma plant sterol locus on the Micronesian Island of Kosrae.

Pinpointing culprit causal variants along signal peaks of genome-wide association studies (GWAS) is challenging. To overcome confounding effects of multiple independent variants at such a locus and narrow the interval for causal allele capture, we developed an approach that maps local shared haplotypes harboring a putative causal variant. We demonstrate our method in an extreme isolate founder population, the pacific Island of Kosrae. We analyzed plasma plant sterol (PPS) levels, a surrogate measure of cholesterol absorption from the intestine, where previous studies have implicated 2p21 mutations in the ATP binding cassette subfamily G members 5 or 8 (ABCG5 or ABCG8) genes. We have previously reported that 11.1% of the islanders are carriers of a frameshift ABCG8 mutation increasing PPS levels in carriers by 50%. GWAS adjusted for this mutation revealed genomewide significant signals along 11 Mb around it. To fine-map this signal, we detected pairwise identity-by-descent haplotypes using our tool GERMLINE and implemented a clustering algorithm to identify haplotypes shared across multiple samples with their unique shared boundaries. A single 526-kb haplotype mapped strongly to PPS levels, dramatically refining the mapped interval. This haplotype spans the ABCG5/ABCG8 genes, is carried by 1.8% of the islanders, and results in a striking 100% increase of PPS in carriers. Resequencing of ABCG5 in these carriers found a D450H missense mutation along the associated haplotype. These findings exemplify the power of haplotype analysis for mapping mutations in isolated populations and specifically for dissecting effects of multiple variants of the same locus.

Genome-wide association studies in an isolated founder population from the Pacific Island of Kosrae.

It has been argued that the limited genetic diversity and reduced allelic heterogeneity observed in isolated founder populations facilitates discovery of loci contributing to both Mendelian and complex disease. A strong founder effect, severe isolation, and substantial inbreeding have dramatically reduced genetic diversity in natives from the island of Kosrae, Federated States of Micronesia, who exhibit a high prevalence of obesity and other metabolic disorders. We hypothesized that genetic drift and possibly natural selection on Kosrae might have increased the frequency of previously rare genetic variants with relatively large effects, making these alleles readily detectable in genome-wide association analysis. However, mapping in large, inbred cohorts introduces analytic challenges, as extensive relatedness between subjects violates the assumptions of independence upon which traditional association test statistics are based. We performed genome-wide association analysis for 15 quantitative traits in 2,906 members of the Kosrae population, using novel approaches to manage the extreme relatedness in the sample. As positive controls, we observe association to known loci for plasma cholesterol, triglycerides, and C-reactive protein and to a compelling candidate loci for thyroid stimulating hormone and fasting plasma glucose. We show that our study is well powered to detect common alleles explaining >/=5% phenotypic variance. However, no such large effects were observed with genome-wide significance, arguing that even in such a severely inbred population, common alleles typically have modest effects. Finally, we show that a majority of common variants discovered in Caucasians have indistinguishable effect sizes on Kosrae, despite the major differences in population genetics and environment.

Timing of the Diagnosis of Autism in African American Children.

OBJECTIVES: African American (AA) children affected by autism spectrum disorder (ASD) experience delays in diagnosis and obstacles to service access, as well as a disproportionate burden of intellectual disability (ID) as documented in surveillance data recently published by the US Centers for Disease Control and Prevention. Our objective in this study was to analyze data from the largest-available repository of diagnostic and phenotypic information on AA children with ASD, and to explore the wide variation in outcome within the cohort as a function of sociodemographic risk and specific obstacles to service access for the purpose of informing a national approach to resolution of these disparities. METHODS: Parents of 584 AA children with autism consecutively enrolled in the Autism Genetic Resource Exchange across 4 US data collection sites completed event history calendar interviews of the diagnostic odysseys for their children with ASD. These data were examined in relation to developmental outcomes of the children with autism and their unaffected siblings. RESULTS: The average age of ASD diagnosis was 64.9 months (±49.6), on average 42.3 months (±45.1) after parents' first concerns about their children's development. The relationship between timing of diagnosis and ASD severity was complex, and ID comorbidity was not predicted in a straightforward manner by familial factors associated with cognitive variation in the general population. CONCLUSIONS: These findings document significant opportunity to expedite diagnosis, the need to further understand causes of ID comorbidity, and the necessity to identify effective approaches to the resolution of disparities in severity-of-outcome for AA children with autism.

Imaging-genetics of sex differences in ASD: distinct effects of OXTR variants on brain connectivity.

Autism spectrum disorder (ASD) is more prevalent in males than in females, but the neurobiological mechanisms that give rise to this sex-bias are poorly understood. The female protective hypothesis suggests that the manifestation of ASD in females requires higher cumulative genetic and environmental risk relative to males. Here, we test this hypothesis by assessing the additive impact of several ASD-associated OXTR variants on reward network resting-state functional connectivity in males and females with and without ASD, and explore how genotype, sex, and diagnosis relate to heterogeneity in neuroendophenotypes. Females with ASD who carried a greater number of ASD-associated risk alleles in the OXTR gene showed greater functional connectivity between the nucleus accumbens (NAcc; hub of the reward network) and subcortical brain areas important for motor learning. Relative to males with ASD, females with ASD and higher OXTR risk-allele-dosage showed increased connectivity between the NAcc, subcortical regions, and prefrontal brain areas involved in mentalizing. This increased connectivity between NAcc and prefrontal cortex mirrored the relationship between genetic risk and brain connectivity observed in neurotypical males showing that, under increased OXTR genetic risk load, females with ASD and neurotypical males displayed increased connectivity between reward-related brain regions and prefrontal cortex. These results indicate that females with ASD differentially modulate the effects of increased genetic risk on brain connectivity relative to males with ASD, providing new insights into the neurobiological mechanisms through which the female protective effect may manifest.