Long COVID-19 and Peripheral Serotonin: A Commentary and Reconsideration

We believe there are serious problems with a recently published and highly publicized paper entitled "Serotonin reduction in post-acute sequelae of viral infection." The blood centrifugation procedure reportedly used by Wong et al would produce plasma that is substantially (over 95%) depleted of platelets. Given this, their published mean plasma serotonin values of 1.2 uM and 2.4 uM for the control/contrast groups appear to be at least 30 to 60 times too high and should be disregarded. The plasma serotonin values reported for the long COVID and viremia patients also should be disregarded, as should any comparisons to the control/contrast groups. We also note that the plasma serotonin means for the two control/contrast groups are not in good agreement. In the "Discussion" section, Wong et al state that their results tend to support the use of selective serotonin reuptake inhibitors (SSRIs) for the treatment of COVID-19, and they encourage further clinical trials of SSRIs. While they state that, "Our animal models demonstrate that serotonin levels can be restored and memory impairment reversed by precursor supplementation or SSRI treatment", it should be noted that no data are presented showing an increase or restoration in circulating serotonin with SSRI administration. In fact, one would expect a marked decline in platelet serotonin due to SSRIs' effective inhibition of the platelet serotonin transporter. Wong et al hypothesize that problems of long COVID arise from too little peripheral serotonin. However, given the frequent presence of a hyperaggregation state in long COVID, and the known augmenting effects of platelet serotonin on platelet aggregation, it is plausible to suggest that reductions in platelet serotonin might be associated with a lessening of the cardiovascular sequelae of COVID-19.

Genome-wide association study identifies 30 obsessive-compulsive disorder associated loci

Obsessive-compulsive disorder (OCD) affects ~1% of the population and exhibits a high SNP-heritability, yet previous genome-wide association studies (GWAS) have provided limited information on the genetic etiology and underlying biological mechanisms of the disorder. We conducted a GWAS meta-analysis combining 53,660 OCD cases and 2,044,417 controls from 28 European-ancestry cohorts revealing 30 independent genome-wide significant SNPs and a SNP-based heritability of 6.7%. Separate GWAS for clinical, biobank, comorbid, and self-report sub-groups found no evidence of sample ascertainment impacting our results. Functional and positional QTL gene-based approaches identified 249 significant candidate risk genes for OCD, of which 25 were identified as putatively causal, highlighting WDR6, DALRD3, CTNND1 and genes in the MHC region. Tissue and single-cell enrichment analyses highlighted hippocampal and cortical excitatory neurons, along with D1- and D2-type dopamine receptor-containing medium spiny neurons, as playing a role in OCD risk. OCD displayed significant genetic correlations with 65 out of 112 examined phenotypes. Notably, it showed positive genetic correlations with all included psychiatric phenotypes, in particular anxiety, depression, anorexia nervosa, and Tourette syndrome, and negative correlations with a subset of the included autoimmune disorders, educational attainment, and body mass index.. This study marks a significant step toward unraveling its genetic landscape and advances understanding of OCD genetics, providing a foundation for future interventions to address this debilitating disorder.

Rare coding variation provides insight into the genetic architecture and phenotypic context of autism

Some individuals with autism spectrum disorder (ASD) carry functional mutations rarely observed in the general population. We explored the genes disrupted by these variants from joint analysis of protein-truncating variants (PTVs), missense variants and copy number variants (CNVs) in a cohort of 63,237 individuals. We discovered 72 genes associated with ASD at false discovery rate (FDR) ≤ 0.001 (185 at FDR ≤ 0.05). De novo PTVs, damaging missense variants and CNVs represented 57.5%, 21.1% and 8.44% of association evidence, while CNVs conferred greatest relative risk. Meta-analysis with cohorts ascertained for developmental delay (DD) (n = 91,605) yielded 373 genes associated with ASD/DD at FDR ≤ 0.001 (664 at FDR ≤ 0.05), some of which differed in relative frequency of mutation between ASD and DD cohorts. The DD-associated genes were enriched in transcriptomes of progenitor and immature neuronal cells, whereas genes showing stronger evidence in ASD were more enriched in maturing neurons and overlapped with schizophrenia-associated genes, emphasizing that these neuropsychiatric disorders may share common pathways to risk.

Maternal Duplication 15q11-13 Syndrome with Autism Spectrum Disorder: Mood Stabilization by Carbamazepine

Objectives: Maternal 15q11-13 duplication syndrome (dup15q) is one of the most frequently observed and penetrant genetic abnormalities associated with autism spectrum disorder (ASD), and commonly presents with psychiatric symptoms and seizures. Although carbamazepine has been reported as effective in managing comorbid seizures in dup15q, it has not been reported to be used as a mood stabilizer in this population. Methods: We retrospectively reviewed the charts of five consecutive patients presenting with previously diagnosed dup15q and ASD seeking treatment for psychiatric symptoms and, in four of the patients, seizures. These were the only patients with dup15q treated with carbamazepine in the Neurodevelopmental Psychopharmacology Clinic at the University of Illinois at Chicago during the review period. Results: During treatment, carbamazepine was found to be more effective than other mood stabilizers in all five patients, and in one case a better antiepileptic. Symptoms consistent with bipolar mood disorder such as hyperactivity, impulsivity, irritability, mood lability, intrusiveness, and pressured speech were improved with carbamazepine in combination with other psychotropic medications. This improvement was greater than with other mood stabilizers, including oxcarbazepine, valproate, and lamotrigine. In one case, valproate paradoxically worsened symptoms. In three cases, anxiety was improved with carbamazepine when used in conjunction with other medications targeting anxiety. Conclusions: In treating five patients with dup15q, carbamazepine more effectively stabilized mood-related symptoms than alternative treatments. Prospective randomized controlled trials are necessary to confirm this observation.

Elevated Polygenic Burden for Autism Spectrum Disorder Is Associated With the Broad Autism Phenotype in Mothers of Individuals With Autism Spectrum Disorder

BACKGROUND

Autism spectrum disorder (ASD) is a multifactorial neurodevelopmental disorder that encompasses a complex and heterogeneous set of traits. Subclinical traits that mirror the core features of ASD, referred to as the broad autism phenotype (BAP), have been documented repeatedly in unaffected relatives and are believed to reflect underlying genetic liability to ASD. The BAP may help inform the etiology of ASD by allowing the stratification of families into more phenotypically and etiologically homogeneous subgroups. This study explores polygenic scores related to the BAP.

METHODS

Phenotypic and genotypic information were obtained from 2614 trios from the Simons Simplex Collection. Polygenic scores of ASD (ASD-PGSs) were generated across the sample to determine the shared genetic overlap between the BAP and ASD. Maternal and paternal ASD-PGSs were explored in relation to BAP traits and their child's ASD symptomatology.

RESULTS

Maternal pragmatic language was related to child's social communicative atypicalities. In fathers, rigid personality was related to increased repetitive behaviors in children. Maternal (but not paternal) ASD-PGSs were related to the pragmatic language and rigid BAP domains.

CONCLUSIONS

Associations emerged between parent and child phenotypes, with more associations emerging in mothers than in fathers. ASD-PGS associations emerged with BAP in mothers only, highlighting the potential for a female protective factor, and implicating the polygenic etiology of ASD-related phenotypes in the BAP.

Properties of beta oscillations in Dup15q syndrome

Background

Duplications of 15q11.2-q13.1 (Dup15q syndrome) are highly penetrant for autism, intellectual disability, hypotonia, and epilepsy. The 15q region harbors genes critical for brain development, particularly UBE3A and a cluster of gamma-aminobutyric acid type A receptor (GABA_AR) genes. We recently described an electrophysiological biomarker of the syndrome, characterized by excessive beta oscillations (12–30 Hz), resembling electroencephalogram (EEG) changes induced by allosteric modulation of GABA_ARs. In this follow-up study, we tested a larger cohort of children with Dup15q syndrome to comprehensively examine properties of this EEG biomarker that would inform its use in future clinical trials, specifically, its (1) relation to basic clinical features, such as age, duplication type, and epilepsy; (2) relation to behavioral characteristics, such as cognition and adaptive function; (3) stability over time; and (4) reproducibility of the signal in clinical EEG recordings.

Methods

We computed EEG power and beta peak frequency (BPF) in a cohort of children with Dup15q syndrome (N = 41, age range 9–189 months). To relate EEG parameters to clinical (study 1) and behavioral features (study 2), we examined age, duplication type, epilepsy, cognition, and daily living skills (DLS) as predictors of beta power and BPF. To evaluate stability over time (study 3), we derived the intraclass correlation coefficients (ICC) from beta power and BPF computed from children with multiple EEG recordings (N = 10, age range 18–161 months). To evaluate reproducibility in a clinical setting (study 4), we derived ICCs from beta power computed from children (N = 8, age range 19–96 months), who had undergone both research EEG and clinical EEG.

Results

The most promising relationships between EEG and clinical traits were found using BPF. BPF was predicted both by epilepsy status (R² = 0.11, p = 0.038) and the DLS component of the Vineland Adaptive Behavior Scale (R² = 0.17, p = 0.01). Beta power and peak frequency showed high stability across repeated visits (beta power ICC = 0.93, BPF ICC = 0.92). A reproducibility analysis revealed that beta power estimates are comparable between research and clinical EEG (ICC = 0.94).

Conclusions

In this era of precision health, with pharmacological and neuromodulatory therapies being developed and tested for specific genetic etiologies of neurodevelopmental disorders, quantification and examination of mechanistic biomarkers can greatly improve clinical trials. To this end, the robust beta oscillations evident in Dup15q syndrome are clinically reproducible and stable over time. With future preclinical and computational studies that will help disentangle the underlying mechanism, it is possible that this biomarker could serve as a robust measure of drug target engagement or a proximal outcome measure in future disease modifying intervention trials.

Large-Scale Exome Sequencing Study Implicates Both Developmental and Functional Changes in the Neurobiology of Autism

We present the largest exome sequencing study of autism spectrum disorder (ASD) to date (n = 35,584 total samples, 11,986 with ASD). Using an enhanced analytical framework to integrate de novo and case-control rare variation, we identify 102 risk genes at a false discovery rate of 0.1 or less. Of these genes, 49 show higher frequencies of disruptive de novo variants in individuals ascertained to have severe neurodevelopmental delay, whereas 53 show higher frequencies in individuals ascertained to have ASD; comparing ASD cases with mutations in these groups reveals phenotypic differences. Expressed early in brain development, most risk genes have roles in regulation of gene expression or neuronal communication (i.e., mutations effect neurodevelopmental and neurophysiological changes), and 13 fall within loci recurrently hit by copy number variants. In cells from the human cortex, expression of risk genes is enriched in excitatory and inhibitory neuronal lineages, consistent with multiple paths to an excitatory-inhibitory imbalance underlying ASD.

Familiality of behavioral flexibility and response inhibition deficits in autism spectrum disorder (ASD)

Background

Diminished cognitive control, including reduced behavioral flexibility and behavioral response inhibition, has been repeatedly documented in autism spectrum disorder (ASD). We evaluated behavioral flexibility and response inhibition in probands and their parents using a family trio design to determine the extent to which these cognitive control impairments represent familial traits associated with ASD.

Methods

We examined 66 individuals with ASD (probands), 135 unaffected biological parents, and 76 typically developing controls. Participants completed a probabilistic reversal learning task (PRL) and a stop-signal task (SST) to assess behavioral flexibility and response inhibition respectively. Rates of PRL and SST errors were examined across groups, within families, and in relation to clinical and subclinical traits of ASD. Based on prior findings that subclinical broader autism phenotypic (BAP) traits may co-segregate within families and reflect heritable risk factors, we also examined whether cognitive control deficits were more prominent in families in which parents showed BAP features (BAP+).

Results

Probands and parents each showed increased rates of PRL and SST errors relative to controls. Error rates across tasks were not related. SST error rates inter-correlated among probands and their parents. PRL errors were more severe in BAP+ parents and their children relative to BAP- parents and their children. For probands of BAP+ parents, PRL and SST error rates were associated with more severe social-communication abnormalities and repetitive behaviors, respectively.

Conclusion

Reduced behavioral flexibility and response inhibition are present among probands and their unaffected parents, but represent unique familial deficits associated with ASD that track with separate clinical issues. Specifically, behavioral response inhibition impairments are familial in ASD and manifest independently from parental subclinical features. In contrast, behavioral flexibility deficits are selectively present in families with BAP characteristics, suggesting they co-segregate in families with parental subclinical social, communication, and rigid personality traits. Together, these findings provide evidence that behavioral flexibility and response inhibition impairments track differentially with ASD risk mechanisms and related behavioral traits.

Clinical and neurocognitive issues associated with Bosch-Boonstra-Schaaf optic atrophy syndrome: A case study

Nuclear receptor subfamily 2 group F member 1 (NR2F1) is an orphan receptor and transcriptional regulator that is involved in neurogenesis, visual processing and development, and cortical patterning. Alterations in NR2F1 cause Bosch-Boonstra-Schaaf optic atrophy syndrome (BBSOAS), a recently described autosomal dominant disorder characterized by intellectual and developmental disabilities and optic atrophy. This study describes the clinical and neurocognitive features of an individual with a de novo nonsense variant in NR2F1 (NM_005654.5:c.82C > T, p.Gln28*), identified by whole exome sequencing. The patient was diagnosed with autism spectrum disorder (ASD) and unlike most previously reported cases, he had no developmental delay, superior verbal abilities (verbal IQ = 141), and high educational attainment despite reduced nonverbal abilities (nonverbal IQ = 63). He had optic nerve hypoplasia with minimal visual impairment as well as mild dysmorphic features. Compared to both age-matched individuals with ASD and healthy controls, the patient showed reductions in manual motor speed, accuracy of saccadic eye movements, and rates of successful behavioral response inhibition. Although the majority of previously reported cases of BBSOAS have been associated with more global intellectual dysfunction, we report on a patient with selective disruption of nonverbal abilities and superior verbal abilities.

Behavioral characterization of dup15q syndrome: Toward meaningful endpoints for clinical trials

Duplication of 15q11.2-q13.1 (dup15q syndrome) is one of the most common copy number variations associated with autism spectrum disorders (ASD) and intellectual disability (ID). As with many neurogenetic conditions, accurate behavioral assessment is challenging due to the level of impairment and heterogeneity across individuals. Large-scale phenotyping studies are necessary to inform future clinical trials in this and similar ID syndromes. This study assessed developmental and behavioral characteristics in a large cohort of children with dup15q syndrome, and examined differences based on genetic subtype and epilepsy status. Participants included 62 children (2.5-18 years). Across individuals, there was a wide range of abilities. Although adaptive behavior was strongly associated with cognitive ability, adaptive abilities were higher than cognitive scores. Measures of ASD symptoms were associated with cognitive ability, while parent report of challenging behavior was not. Both genetic subtype and epilepsy were related to degree of impairment across cognitive, language, motor, and adaptive domains. Children with isodicentric duplications and epilepsy showed the greatest impairment, while children with interstitial duplications showed the least. On average, participants with epilepsy experienced seizures over 53% of their lives, and half of children with epilepsy had infantile spasms. Parents of children with isodicentric duplications reported more concerns regarding challenging behaviors. Future trials in ID syndromes should employ a flexible set of assessments, allowing each participant to receive assessments that capture their skills. Multiple sources of information should be considered, and the impact of language and cognitive ability should be taken into consideration when interpreting results.

Maternal Serotonin Levels Are Associated With Cognitive Ability and Core Symptoms in Autism Spectrum Disorder

OBJECTIVE

The serotonin (5-hydroxytryptamine [HT]) system has long been implicated in autism spectrum disorder (ASD). Whole-blood 5-HT level (WB5-HT) is a stable, heritable biomarker that is elevated in more than 25% of children with ASD. Recent findings indicate that the maternal 5-HT system may influence embryonic neurodevelopment, but maternal WB5-HT has not been examined in relation to ASD phenotypes.

METHOD

WB5-HT levels were obtained from 181 individuals (3-27 years of age) diagnosed with ASD, 99 of their fathers, and 119 of their mothers. Standardized assessments were used to evaluate cognitive, behavioral, and language phenotypes.

RESULTS

Exploratory regression analyses found relationships between maternal WB5-HT and nonverbal IQ (NVIQ), Autism Diagnostic Interview-Revised (ADI-R) Nonverbal Communication Algorithm scores, and overall adaptive function on the Vineland Adaptive Behavior Scales-II. Latent class analysis identified a three-class structure in the assessment data, describing children with low, intermediate, and high severity across measures of behavior, cognition, and adaptive function. Mean maternal WB5-HT differed across classes, with the lowest maternal WB5-HT levels seen in the highest-severity group (Welch F_2,46.048 = 17.394, p < .001). Paternal and proband WB5-HT did not differ between classes.

CONCLUSION

Maternal WB5-HT is associated with neurodevelopmental outcomes in offspring with ASD. Prospective, longitudinal studies will be needed to better understand the relationship between the function of the maternal serotonin system during pregnancy and brain development. Further studies in animal models may be able to reveal the mechanisms underlying these findings.

Cognitive mechanisms of inhibitory control deficits in autism spectrum disorder

BACKGROUND

Inhibitory control deficits are common in autism spectrum disorder (ASD) and associated with more severe repetitive behaviors. Inhibitory control deficits may reflect slower execution of stopping processes, or a reduced ability to delay the onset of behavioral responses in contexts of uncertainty. Previous studies have documented relatively spared stopping processes in ASD, but whether inhibitory control deficits in ASD reflect failures to delay response onset has not been systematically assessed. Further, while improvements in stopping abilities and response slowing are seen through adolescence/early adulthood in health, their development in ASD is less clear.

METHODS

A stop-signal test (SST) was administered to 121 individuals with ASD and 76 age and IQ-matched healthy controls (ages 5-28). This test included 'GO trials' in which participants pressed a button when a peripheral target appeared and interleaved 'STOP trials' in which they were cued to inhibit button-presses when a stop-signal appeared at variable times following the GO cue. STOP trial accuracy, RT of the stopping process (SSRT), and reaction time (RT) slowing during GO trials were examined.

RESULTS

Relative to controls, individuals with ASD had reduced accuracy on STOP trials. SSRTs were similar across control and ASD participants, but RT slowing was reduced in patients compared to controls. Age-related increases in stopping ability and RT slowing were attenuated in ASD. Reduced stopping accuracy and RT slowing were associated with more severe repetitive behaviors in ASD.

DISCUSSION

Our findings show that inhibitory control deficits in ASD involve failures to strategically delay behavioral response onset. These results suggest that reduced preparatory behavioral control may underpin inhibitory control deficits as well as repetitive behaviors in ASD. Typical age-related improvements in inhibitory control during late childhood/early adolescence are reduced in ASD, highlighting an important developmental window during which treatments may mitigate cognitive alterations contributing to repetitive behaviors.

De novo unbalanced translocation (4p duplication/8p deletion) in a patient with autism, OCD, and overgrowth syndrome

Chromosomal abnormalities, such as unbalanced translocations and copy number variants (CNVs), are found in autism spectrum disorders (ASDs) [Sanders et al. (2011) Neuron 70: 863-885]. Many chromosomal abnormalities, including sub microscopic genomic deletions and duplications, are missed by G-banded karyotyping or Fragile X screening alone and are picked up by chromosomal microarrays [Shen et al. (2010) Pediatrics 125: e727-735]. Translocations involving chromosomes 4 and 8 are possibly the second most frequent translocation in humans and are often undetected in routine cytogenetics [Giglio et al. (2002) Circulation 102: 432-437]. Deletions of 4p16 have been associated with Wolf-Hirschhorn syndrome while 4p16 duplications have been associated with an overgrowth syndrome and mild to moderate mental retardation [Partington et al. (1997) Journal of Medical Genetics 34: 719-728]. The 8p23.3 region contains the autism candidate gene DLGAP2, which can contribute to autism when disrupted [Marshall et al. (2008) The American Journal of Human Genetics 82: 477-488] . There has been a case report of a family with autism spectrum disorder (ASD), prominent obsessional behavior, and overgrowth in patients with der (8) t (4;8) p (16;23) [Partington et al. (1997)]. This is an independent report of a male patient with autism, obsessive compulsive disorder (OCD), attention-deficit hyperactivity disorder (ADHD), and an overgrowth syndrome, whose de novo unbalanced translocation der (8) t (4;8) p (16.1→ter; 23.1→ter) was initially missed by routine cytogenetics but detected with SNP microarray, allowing higher resolution of translocation breakpoints.

Is there sexual dimorphism of hyperserotonemia in autism spectrum disorder?

Approximately 30% of individuals with autism spectrum disorder (ASD) have elevated whole blood serotonin (5-HT) levels. Genetic linkage and association studies of ASD and of whole blood 5-HT levels as a quantitative trait have revealed sexual dimorphism. Few studies have examined the presence of a sex difference on hyperserotonemia within ASD. To assess whether the rate of hyperserotonemia is different in males than in females with ASD, we measured whole blood 5-HT levels in 292 children and adolescents with ASD, the largest sample in which this biomarker has been assessed. Based upon previous work suggesting that hyperserotonemia is more common prior to puberty, we focused our analysis on the 182 pre-pubertal children with ASD. 42% of pre-pubertal participants were within the hyperserotonemia range. In this population, we found that males were significantly more likely to manifest hyperserotonemia than females (P = 0.03). As expected, no significant difference was found in the post-pubertal population. Additional work will be needed to replicate this intriguing finding and to understand whether it could potentially explain differences in patterns of ASD risk between males and females. Autism Res 2017, 10: 1417-1423. © 2017 International Society for Autism Research, Wiley Periodicals, Inc.

Leveraging blood serotonin as an endophenotype to identify de novo and rare variants involved in autism

BACKGROUND

Autism spectrum disorder (ASD) is one of the most highly heritable neuropsychiatric disorders, but underlying molecular mechanisms are still unresolved due to extreme locus heterogeneity. Leveraging meaningful endophenotypes or biomarkers may be an effective strategy to reduce heterogeneity to identify novel ASD genes. Numerous lines of evidence suggest a link between hyperserotonemia, i.e., elevated serotonin (5-hydroxytryptamine or 5-HT) in whole blood, and ASD. However, the genetic determinants of blood 5-HT level and their relationship to ASD are largely unknown.

METHODS

In this study, pursuing the hypothesis that de novo variants (DNVs) and rare risk alleles acting in a recessive mode may play an important role in predisposition of hyperserotonemia in people with ASD, we carried out whole exome sequencing (WES) in 116 ASD parent-proband trios with most (107) probands having 5-HT measurements.

RESULTS

Combined with published ASD DNVs, we identified USP15 as having recurrent de novo loss of function mutations and discovered evidence supporting two other known genes with recurrent DNVs (FOXP1 and KDM5B). Genes harboring functional DNVs significantly overlap with functional/disease gene sets known to be involved in ASD etiology, including FMRP targets and synaptic formation and transcriptional regulation genes. We grouped the probands into High-5HT and Normal-5HT groups based on normalized serotonin levels, and used network-based gene set enrichment analysis (NGSEA) to identify novel hyperserotonemia-related ASD genes based on LoF and missense DNVs. We found enrichment in the High-5HT group for a gene network module (DAWN-1) previously implicated in ASD, and this points to the TGF-β pathway and cell junction processes. Through analysis of rare recessively acting variants (RAVs), we also found that rare compound heterozygotes (CHs) in the High-5HT group were enriched for loci in an ASD-associated gene set. Finally, we carried out rare variant group-wise transmission disequilibrium tests (gTDT) and observed significant association of rare variants in genes encoding a subset of the serotonin pathway with ASD.

CONCLUSIONS

Our study identified USP15 as a novel gene implicated in ASD based on recurrent DNVs. It also demonstrates the potential value of 5-HT as an effective endophenotype for gene discovery in ASD, and the effectiveness of this strategy needs to be further explored in studies of larger sample sizes.

Whole genome sequencing resource identifies 18 new candidate genes for autism spectrum disorder

We are performing whole-genome sequencing of families with autism spectrum disorder (ASD) to build a resource (MSSNG) for subcategorizing the phenotypes and underlying genetic factors involved. Here we report sequencing of 5,205 samples from families with ASD, accompanied by clinical information, creating a database accessible on a cloud platform and through a controlled-access internet portal. We found an average of 73.8 de novo single nucleotide variants and 12.6 de novo insertions and deletions or copy number variations per ASD subject. We identified 18 new candidate ASD-risk genes and found that participants bearing mutations in susceptibility genes had significantly lower adaptive ability (P = 6 × 10^-4). In 294 of 2,620 (11.2%) of ASD cases, a molecular basis could be determined and 7.2% of these carried copy number variations and/or chromosomal abnormalities, emphasizing the importance of detecting all forms of genetic variation as diagnostic and therapeutic targets in ASD.

Pharmacogenomics of autism spectrum disorder

Autism spectrum disorder (ASD) is characterized by persistent deficits in social communication and interactions as well as restricted, repetitive behaviors and interests. Pharmacologic interventions are often needed to manage irritability, aggressive behaviors and hyperactivity. Pharmacogenomic studies have investigated genetic associations with treatment response and side effects in an attempt to better understand drug mechanisms in hopes of optimizing the balance of symptom improvement versus side effects. The majority of pharmacogenomic studies to date have focused on antipsychotics, antidepressants and stimulants that are the most commonly utilized medication classes for ASD. This review is a comprehensive examination of the existing pharmacogenomic studies in ASD highlighting the current state of knowledge regarding genetic variation influencing pharmacokinetics and pharmacodynamics, and associated clinical outcomes.

ASD and Genetic Associations with Receptors for Oxytocin and Vasopressin-AVPR1A, AVPR1B, and OXTR

Background: There are limited treatments available for autism spectrum disorder (ASD). Studies have reported significant associations between the receptor genes of oxytocin (OT) and vasopressin (AVP) and ASD diagnosis, as well as ASD-related phenotypes. Researchers have also found the manipulation of these systems affects social and repetitive behaviors, core characteristics of ASD. Consequently, research involving the oxytocin/vasopressin pathways as intervention targets has increased. Therefore, further examination into the relationship between these neuropeptides and ASD was undertaken. In this study, we examined associations between variants in the receptor genes of vasopressin (AVPR1A, AVPR1B), oxytocin (OXTR), and ASD diagnosis along with related subphenotypes.

Methods: Probands were assessed using Autism Diagnostic Interview-Revised, Autism Diagnostic Observation Schedule, and clinical DSM-IV-TR criteria. Single nucleotide polymorphisms (SNPs) in AVPR1B and OXTR, and microsatellites in AVPR1A were genotyped in ~200 families with a proband with ASD. Family-based association testing (FBAT) was utilized to determine associations between variants and ASD. Haplotypes composed of OXTR SNPs (i.e., rs53576-rs2254298-rs2268493) were also analyzed due to previously published associations.

Results: Using the additive inheritance model in FBAT we found associations between AVPR1B SNPs (rs28632197, p = 0.005, rs35369693, p = 0.025) and diagnosis. As in other studies, OXTR rs2268493 (p = 0.050) was associated with diagnosis. rs2268493 was also associated with ASD subphenotypes of social withdrawal (p = 0.013) and Insistence on Sameness (p = 0.039). Further analyses demonstrated that the haplotype, rs2254298–rs2268493 was found to be significantly associated with diagnosis (A-T; p = 0.026). FBAT was also used to analyze AVPR1A microsatellites (RS1 and RS3). Both length variants were found to be associated with restrictive, repetitive behaviors, but not overall diagnosis. Correction for multiple comparisons was performed for SNPs tested in each gene region, only AVPR1B SNPs remained significantly associated with ASD diagnosis.

Conclusions: Autism is a heterogeneous disorder with many genes and pathways that contribute to its development. SNPs and microsatellites in the receptor genes of OT and AVP are associated with ASD diagnosis and measures of social behavior as well as restricted repetitive behaviors. We reported a novel association with ASD and AVPR1B SNPs. Understanding of genotype-phenotype relationships may be helpful in the development of pharmacological interventions for the OT/AVP system.

STAMS: STRING-assisted module search for genome wide association studies and application to autism

Motivation: Analyzing genome wide association data in the context of biological pathways helps us understand how genetic variation influences phenotype and increases power to find associations. However, the utility of pathway-based analysis tools is hampered by undercuration and reliance on a distribution of signal across all of the genes in a pathway. Methods that combine genome wide association results with genetic networks to infer the key phenotype-modulating subnetworks combat these issues, but have primarily been limited to network definitions with yes/no labels for gene-gene interactions. A recent method (EW_dmGWAS) incorporates a biological network with weighted edge probability by requiring a secondary phenotype-specific expression dataset. In this article, we combine an algorithm for weighted-edge module searching and a probabilistic interaction network in order to develop a method, STAMS, for recovering modules of genes with strong associations to the phenotype and probable biologic coherence. Our method builds on EW_dmGWAS but does not require a secondary expression dataset and performs better in six test cases.

Results: We show that our algorithm improves over EW_dmGWAS and standard gene-based analysis by measuring precision and recall of each method on separately identified associations. In the Wellcome Trust Rheumatoid Arthritis study, STAMS-identified modules were more enriched for separately identified associations than EW_dmGWAS (STAMS P-value 3.0 × 10⁻⁴; EW_dmGWAS- P-value = 0.8). We demonstrate that the area under the Precision-Recall curve is 5.9 times higher with STAMS than EW_dmGWAS run on the Wellcome Trust Type 1 Diabetes data.

Availability and Implementation: STAMS is implemented as an R package and is freely available at https://simtk.org/projects/stams.

Contact:rbaltman@stanford.edu

Supplementary information:Supplementary data are available at Bioinformatics online.

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.

Commentary on "Platelet Studies in Autism Spectrum Disorder Patients and First-Degree Relatives"

We comment on the recent report entitled "Platelet Studies in Autism Spectrum Disorder Patients and First-Degree Relatives" [Molecular Autism 2015;6:57]. We find it commendable that the authors have investigated platelet factors potentially involved in the well-replicated observation of platelet hyperserotonemia in autism. However, we believe the results need a fuller discussion in the context of prior studies, think that certain aspects of the interpretation need to be reassessed, and attempt to provide a framework for further research in this area.

Confirmation of the Factor Structure and Measurement Invariance of the Children's Scale of Hostility and Aggression: Reactive/Proactive in Clinic-Referred Children With and Without Autism Spectrum Disorder

OBJECTIVE

The measurement of aggression in its different forms (e.g., physical and verbal) and functions (e.g., impulsive and instrumental) is given little attention in subjects with developmental disabilities (DD). In this study, we confirm the factor structure of the Children's Scale for Hostility and Aggression: Reactive/Proactive (C-SHARP) and demonstrate measurement invariance (consistent performance across clinical groups) between clinic-referred groups with and without autism spectrum disorder (ASD). We also provide evidence of the construct validity of the C-SHARP.

METHODS

Caregivers provided C-SHARP, Child Behavior Checklist (CBCL), and Proactive/Reactive Rating Scale (PRRS) ratings for 644 children, adolescents, and young adults 2-21 years of age. Five types of measurement invariance were evaluated within a confirmatory factor analytic framework. Associations among the C-SHARP, CBCL, and PRRS were explored.

RESULTS

The factor structure of the C-SHARP had a good fit to the data from both groups, and strict measurement invariance between ASD and non-ASD groups was demonstrated (i.e., equivalent structure, factor loadings, item intercepts and residuals, and latent variance/covariance between groups). The C-SHARP Problem Scale was more strongly associated with CBCL Externalizing than with CBCL Internalizing, supporting its construct validity. Subjects classified with the PRRS as both Reactive and Proactive had significantly higher C-SHARP Proactive Scores than those classified as Reactive only, who were rated significantly higher than those classified by the PRRS as Neither Reactive nor Proactive. A similar pattern was observed for the C-SHARP Reactive Score.

CONCLUSIONS

This study provided evidence of the validity of the C-SHARP through confirmation of its factor structure and its relationship with more established scales. The demonstration of measurement invariance demonstrates that differences in C-SHARP factor scores were the result of differences in the construct rather than to error or unmeasured/nuisance variables. These data suggest that the C-SHARP is useful for quantifying subtypes of aggressive behavior in children, adolescents, and young adults with DD.

Separating Family-Level and Direct Exposure Effects of Smoking During Pregnancy on Offspring Externalizing Symptoms: Bridging the Behavior Genetic and Behavior Teratologic Divide

Maternal smoking during pregnancy (MSDP) has been robustly associated with externalizing problems and their developmental precursors in offspring in studies using behavioral teratologic designs (Wakschlag et al., Am J Public Health 92(6):966-974, 2002; Espy et al., Dev Psychol 47(1):153-169, 2011). In contrast, the use of behavior genetic approaches has shown that the effects commonly attributed to MSDP can be explained by family-level variables (D'Onofrio et al., Dev Psychopathol 20(01):139-164, 2008). Reconciling these conflicting findings requires integration of these study designs. We utilize longitudinal data on a preschool proband and his/her sibling from the Midwest Infant Development Study-Preschool (MIDS-P) to test for teratologic and family level effects of MSDP. We find considerable variation in prenatal smoking patterns both within and across pregnancies within families, indicating that binary smoking measures are not sufficiently capturing exposure. Structural equation models indicate that both conduct disorder and oppositional defiant disorder symptoms showed unique effects of MSDP over and above family level effects. Blending high quality exposure measurement with a within-family design suggests that it is premature to foreclose the possibility of a teratologic effect of MSDP on externalizing problems. Implications and recommendations for future studies are discussed.

The phenotypic manifestations of rare genic CNVs in autism spectrum disorder

Significant evidence exists for the association between copy number variants (CNVs) and Autism Spectrum Disorder (ASD); however, most of this work has focused solely on the diagnosis of ASD. There is limited understanding of the impact of CNVs on the 'sub-phenotypes' of ASD. The objective of this paper is to evaluate associations between CNVs in differentially brain expressed (DBE) genes or genes previously implicated in ASD/intellectual disability (ASD/ID) and specific sub-phenotypes of ASD. The sample consisted of 1590 cases of European ancestry from the Autism Genome Project (AGP) with a diagnosis of an ASD and at least one rare CNV impacting any gene and a core set of phenotypic measures, including symptom severity, language impairments, seizures, gait disturbances, intelligence quotient (IQ) and adaptive function, as well as paternal and maternal age. Classification analyses using a non-parametric recursive partitioning method (random forests) were employed to define sets of phenotypic characteristics that best classify the CNV-defined groups. There was substantial variation in the classification accuracy of the two sets of genes. The best variables for classification were verbal IQ for the ASD/ID genes, paternal age at birth for the DBE genes and adaptive function for de novo CNVs. CNVs in the ASD/ID list were primarily associated with communication and language domains, whereas CNVs in DBE genes were related to broader manifestations of adaptive function. To our knowledge, this is the first study to examine the associations between sub-phenotypes and CNVs genome-wide in ASD. This work highlights the importance of examining the diverse sub-phenotypic manifestations of CNVs in ASD, including the specific features, comorbid conditions and clinical correlates of ASD that comprise underlying characteristics of the disorder.

Significant neuronal soma volume deficit in the limbic system in subjects with 15q11.2-q13 duplications

Introduction

Autism is diagnosed in numerous genetic and genomic developmental disorders associated with an overlap in high-risk genes and loci that underlie intellectual disability (ID) and epilepsy. The aim of this stereological study of neuronal soma volume in 25 brain structures and their subdivisions in eight individuals 9 to 26 years of age who were diagnosed with chromosome 15q11.2-13.1 duplication syndrome [dup(15)], autism, ID and epilepsy; eight age-matched subjects diagnosed with autism of unknown etiology (idiopathic autism) and seven control individuals was to establish whether defects of neuronal soma growth are a common denominator of developmental pathology in idiopathic and syndromic autism and how genetic modifications alter the trajectory of neuronal soma growth in dup(15) autism.

Results

Application of the Nucleator software to estimate neuronal size revealed significant neuronal soma volume deficits in 11 of 25 structures and their subregions (44 %) in subjects diagnosed with dup(15) autism, including consistent neuronal soma volume deficits in the limbic system (sectors CA2, 3 and 4 in Ammon’s horn, the second and third layers of the entorhinal cortex and in the amygdala), as well as in the thalamus, nucleus accumbens, external globus pallidus, and Ch3 nucleus in the magnocellular basal complex, and in the inferior olive in the brainstem. The second feature distinguishing dup(15) autism was persistent neuronal soma deficits in adolescents and young adults, whereas in idiopathic autism, neuronal volume deficit is most prominent in 4- to 8-year-old children but affects only a few brain regions in older subjects.

Conclusions

This study demonstrates that alterations in the trajectory of neuronal growth throughout the lifespan are a core pathological features of idiopathic and syndromic autism. However, dup(15) causes persistent neuronal volume deficits in adolescence and adulthood, with prominent neuronal growth deficits in all major compartments of the limbic system. The more severe neuronal nuclear and cytoplasic volume deficits in syndromic autism found in this study and the more severe focal developmental defects in the limbic system in dup(15) previously reported in this cohort may contribute to the high prevalence of early onset intractable epilepsy and sudden unexpected death in epilepsy.

An ontology for Autism Spectrum Disorder (ASD) to infer ASD phenotypes from Autism Diagnostic Interview-Revised data

OBJECTIVE

Our goal is to create an ontology that will allow data integration and reasoning with subject data to classify subjects, and based on this classification, to infer new knowledge on Autism Spectrum Disorder (ASD) and related neurodevelopmental disorders (NDD). We take a first step toward this goal by extending an existing autism ontology to allow automatic inference of ASD phenotypes and Diagnostic & Statistical Manual of Mental Disorders (DSM) criteria based on subjects' Autism Diagnostic Interview-Revised (ADI-R) assessment data.

MATERIALS AND METHODS

Knowledge regarding diagnostic instruments, ASD phenotypes and risk factors was added to augment an existing autism ontology via Ontology Web Language class definitions and semantic web rules. We developed a custom Protégé plugin for enumerating combinatorial OWL axioms to support the many-to-many relations of ADI-R items to diagnostic categories in the DSM. We utilized a reasoner to infer whether 2642 subjects, whose data was obtained from the Simons Foundation Autism Research Initiative, meet DSM-IV-TR (DSM-IV) and DSM-5 diagnostic criteria based on their ADI-R data.

RESULTS

We extended the ontology by adding 443 classes and 632 rules that represent phenotypes, along with their synonyms, environmental risk factors, and frequency of comorbidities. Applying the rules on the data set showed that the method produced accurate results: the true positive and true negative rates for inferring autistic disorder diagnosis according to DSM-IV criteria were 1 and 0.065, respectively; the true positive rate for inferring ASD based on DSM-5 criteria was 0.94.

DISCUSSION

The ontology allows automatic inference of subjects' disease phenotypes and diagnosis with high accuracy.

CONCLUSION

The ontology may benefit future studies by serving as a knowledge base for ASD. In addition, by adding knowledge of related NDDs, commonalities and differences in manifestations and risk factors could be automatically inferred, contributing to the understanding of ASD pathophysiology.

Pharmacogenetic Study of Serotonin Transporter and 5HT2A Genotypes in Autism

OBJECTIVE

The purpose of this study was to determine whether polymorphisms in the serotonin transporter (SLC6A4) and serotonin-2A receptor (HTR2A) genes are associated with response to escitalopram in patients with autism spectrum disorder (ASD).

METHODS

Forty-four participants with ASD were enrolled in a 6 week, forced titration, open label examination of the selective serotonin reuptake inhibitor (SSRI) escitalopram. Doses increased at weekly intervals starting at 2.5mg daily with a maximum possible dose of 20 mg daily achieved by the end of the study. If adverse events were experienced, participants subsequently received the previously tolerated dose for the duration of study. SLC6A4 (5-HTTLPR) and HTR2A (rs7997012) genotype groups were assessed in relation to treatment outcomes and drug doses.

RESULTS

Insistence on sameness and irritability symptoms significantly improved over the course of the 6 week treatment period (p<0.0001) in this open-label trial. There were no significant differences observed in the rate of symptom improvement over time across genotype groups. Similarly, dosing trajectory was not significantly associated with genotype groups.

CONCLUSIONS

Previous studies have identified SLC6A4 and HTR2A associations with SSRI response in patients with depression and 5-HTTLPR (SLC6A4) associations with escitalopram response in ASD. We did not observe evidence for similar relationships in this ASD study.

Cognitive set shifting deficits and their relationship to repetitive behaviors in autism spectrum disorder

The neurocognitive impairments associated with restricted and repetitive behaviors (RRBs) in autism spectrum disorder (ASD) are not yet clear. Prior studies indicate that individuals with ASD show reduced cognitive flexibility, which could reflect difficulty shifting from a previously learned response pattern or a failure to maintain a new response set. We examined different error types on a test of set-shifting completed by 60 individuals with ASD and 55 age- and nonverbal IQ-matched controls. Individuals with ASD were able to initially shift sets, but they exhibited difficulty maintaining new response sets. Difficulty with set maintenance was related to increased severity of RRBs. General difficulty maintaining new response sets and a heightened tendency to revert to old preferences may contribute to RRBs.

Integrin β3 Haploinsufficiency Modulates Serotonin Transport and Antidepressant-Sensitive Behavior in Mice

Converging lines of evidence have identified genetic interactions between the serotonin transporter (SERT) gene and ITGB3, which encodes the β3 subunit that forms the αIIbβ3 and αvβ3 integrin receptor complexes. Here we examine the consequences of haploinsufficiency in the mouse integrin β3 subunit gene (Itgb3) on SERT function and selective 5-hydroxytryptamine (5-HT) reuptake inhibitor (SSRI) effectiveness in vivo. Biochemical fractionation studies and immunofluorescent staining of murine brain slices reveal that αvβ3 receptors and SERTs are enriched in presynaptic membranes from several brain regions and that αvβ3 colocalizes with a subpopulation of SERT-containing synapses in raphe nuclei. Notably, we establish that loss of a single allele of Itgb3 in murine neurons is sufficient to decrease 5-HT uptake by SERT in midbrain synaptosomes. Pharmacological assays to elucidate the αvβ3-mediated mechanism of reduced SERT function indicate that decreased integrin β3 subunit expression scales down the population size of active SERT molecules and, as a consequence, lowers the effective dose of SSRIs. These data are consistent with the existence of a subpopulation of SERTs that are tightly modulated by integrin αvβ3 and significantly contribute to global SERT function at 5-HT synapses in the midbrain. Importantly, our screen of a normal human population for single nucleotide polymorphisms in human ITGB3 identified a variant associated with reductions in integrin β3 expression levels that parallel our mouse findings. Thus, polymorphisms in human ITGB3 may contribute to the differential responsiveness of select patients to SSRIs.

A genome-wide association study of autism using the Simons Simplex Collection: Does reducing phenotypic heterogeneity in autism increase genetic homogeneity?

BACKGROUND

Phenotypic heterogeneity in autism has long been conjectured to be a major hindrance to the discovery of genetic risk factors, leading to numerous attempts to stratify children based on phenotype to increase power of discovery studies. This approach, however, is based on the hypothesis that phenotypic heterogeneity closely maps to genetic variation, which has not been tested. Our study examines the impact of subphenotyping of a well-characterized autism spectrum disorder (ASD) sample on genetic homogeneity and the ability to discover common genetic variants conferring liability to ASD.

METHODS

Genome-wide genotypic data of 2576 families from the Simons Simplex Collection were analyzed in the overall sample and phenotypic subgroups defined on the basis of diagnosis, IQ, and symptom profiles. We conducted a family-based association study, as well as estimating heritability and evaluating allele scores for each phenotypic subgroup.

RESULTS

Association analyses revealed no genome-wide significant association signal. Subphenotyping did not increase power substantially. Moreover, allele scores built from the most associated single nucleotide polymorphisms, based on the odds ratio in the full sample, predicted case status in subsets of the sample equally well and heritability estimates were very similar for all subgroups.

CONCLUSIONS

In genome-wide association analysis of the Simons Simplex Collection sample, reducing phenotypic heterogeneity had at most a modest impact on genetic homogeneity. Our results are based on a relatively small sample, one with greater homogeneity than the entire population; if they apply more broadly, they imply that analysis of subphenotypes is not a productive path forward for discovering genetic risk variants in ASD.

CYP2A6 Longitudinal Effects in Young Smokers

INTRODUCTION

The present study sought to identify time-dependent within-participant effects of CYP2A6 genotypes on smoking frequency and nicotine dependence in young smokers.

METHODS

Predicted nicotine metabolic rate based on CYP2A6 diplotypes (CYP2A6 diplotype predicted rate [CDPR]) was partitioned into Normal, Intermediate, and Slow categories using a metabolism metric. Growth-curve models characterized baseline and longitudinal CDPR effects with data from eight longitudinal assessments during a 6-year period (from approximately age 16-22) in young smokers of European descent (N = 296, 57% female) who had smoked less than 100 cigarettes lifetime at baseline and more than that amount by Year 6. Phenotypes were number of days smoked during the previous 30 days and a youth version of the Nicotine Dependence Syndrome Scale (NDSS). A zero-inflated Poisson growth-curve model was used to account for the preponderance of zero days smoked.

RESULTS

At baseline, Intermediate CDPR was a risk factor relative to both Normal and Slow CDPR for smoking frequency and the NDSS. Slow CDPR was associated with the highest probability of smoking discontinuation at baseline. However, due to CDPR time trend differences, by young adulthood these baseline effects had been reordered such that the greatest risks for smoking frequency and the NDSS were associated with Normal CDPR.

CONCLUSIONS

Reduced metabolism CYP2A6 genotypes are associated with both risk and protective effects in novice smokers. However, differences in the time-by-CDPR effects result in a reordering of genotype effects such that normal metabolism becomes the risk variant by young adulthood, as has been reliably reported in older smokers.

A haplotype-based framework for group-wise transmission/disequilibrium tests for rare variant association analysis

MOTIVATION

A major focus of current sequencing studies for human genetics is to identify rare variants associated with complex diseases. Aside from reduced power of detecting associated rare variants, controlling for population stratification is particularly challenging for rare variants. Transmission/disequilibrium tests (TDT) based on family designs are robust to population stratification and admixture, and therefore provide an effective approach to rare variant association studies to eliminate spurious associations. To increase power of rare variant association analysis, gene-based collapsing methods become standard approaches for analyzing rare variants. Existing methods that extend this strategy to rare variants in families usually combine TDT statistics at individual variants and therefore lack the flexibility of incorporating other genetic models.

RESULTS

In this study, we describe a haplotype-based framework for group-wise TDT (gTDT) that is flexible to encompass a variety of genetic models such as additive, dominant and compound heterozygous (CH) (i.e. recessive) models as well as other complex interactions. Unlike existing methods, gTDT constructs haplotypes by transmission when possible and inherently takes into account the linkage disequilibrium among variants. Through extensive simulations we showed that type I error was correctly controlled for rare variants under all models investigated, and this remained true in the presence of population stratification. Under a variety of genetic models, gTDT showed increased power compared with the single marker TDT. Application of gTDT to an autism exome sequencing data of 118 trios identified potentially interesting candidate genes with CH rare variants.

AVAILABILITY AND IMPLEMENTATION

We implemented gTDT in C++ and the source code and the detailed usage are available on the authors' website (https://medschool.vanderbilt.edu/cgg).

CONTACT

bingshan.li@vanderbilt.edu or wei.chen@chp.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Cross-disorder genome-wide analyses suggest a complex genetic relationship between Tourette's syndrome and OCD

OBJECTIVE

Obsessive-compulsive disorder (OCD) and Tourette's syndrome are highly heritable neurodevelopmental disorders that are thought to share genetic risk factors. However, the identification of definitive susceptibility genes for these etiologically complex disorders remains elusive. The authors report a combined genome-wide association study (GWAS) of Tourette's syndrome and OCD.

METHOD

The authors conducted a GWAS in 2,723 cases (1,310 with OCD, 834 with Tourette's syndrome, 579 with OCD plus Tourette's syndrome/chronic tics), 5,667 ancestry-matched controls, and 290 OCD parent-child trios. GWAS summary statistics were examined for enrichment of functional variants associated with gene expression levels in brain regions. Polygenic score analyses were conducted to investigate the genetic architecture within and across the two disorders.

RESULTS

Although no individual single-nucleotide polymorphisms (SNPs) achieved genome-wide significance, the GWAS signals were enriched for SNPs strongly associated with variations in brain gene expression levels (expression quantitative loci, or eQTLs), suggesting the presence of true functional variants that contribute to risk of these disorders. Polygenic score analyses identified a significant polygenic component for OCD (p=2×10(-4)), predicting 3.2% of the phenotypic variance in an independent data set. In contrast, Tourette's syndrome had a smaller, nonsignificant polygenic component, predicting only 0.6% of the phenotypic variance (p=0.06). No significant polygenic signal was detected across the two disorders, although the sample is likely underpowered to detect a modest shared signal. Furthermore, the OCD polygenic signal was significantly attenuated when cases with both OCD and co-occurring Tourette's syndrome/chronic tics were included in the analysis (p=0.01).

CONCLUSIONS

Previous work has shown that Tourette's syndrome and OCD have some degree of shared genetic variation. However, the data from this study suggest that there are also distinct components to the genetic architectures of these two disorders. Furthermore, OCD with co-occurring Tourette's syndrome/chronic tics may have different underlying genetic susceptibility compared with OCD alone.

Self-injury in autism spectrum disorder: an effect of serotonin transporter gene promoter variants

Self-injurious behavior in autism spectrum disorder (ASD) has been associated with lower whole blood serotonin levels and the role of serotonin transporter gene promoter region (5HTTLPR) polymorphisms is of interest because of their effects on transporter functioning. This study examined the association between self-injurious behavior in ASD and allelic frequencies of 5HTTLPR. Sixty-four children and adolescents with ASD who were not taking serotonergic medication at the time of the assessment were included in the analysis. Self-injury was assessed using the Autism Diagnostic Interview-Revised (ADI-R) and whole blood serotonin levels were measured using high-pressure liquid chromatography (HPLC) with fluorometic detection. DNA was extracted from saliva and PCR amplified with fluorescent primers. Self-injury significantly increased with the number of La alleles of the 5HTTLPR and decreased with the number of Lg alleles. Self-injury in ASD may be associated with a specific genotype of the serotonin transporter gene promoter region. Future studies should continue to explore subgroups to clarify the underlying clinical and genetic heterogeneity in ASD.

UGT1A and UGT2B genetic variation alters nicotine and nitrosamine glucuronidation in european and african american smokers

BACKGROUND

Identifying sources of variation in the nicotine and nitrosamine metabolic inactivation pathways is important to understanding the relationship between smoking and cancer risk. Numerous UGT1A and UGT2B enzymes are implicated in nicotine and nitrosamine metabolism in vitro; however, little is known about their roles in vivo.

METHODS

Within UGT1A1, UGT1A4, UGT1A9, UGT2B7, UGT2B10, and UGT2B17, 47 variants were genotyped, including UGT2B10*2 and UGT2B17*2. The association between variation in these UGTs and glucuronidation activity within European and African American current smokers (n = 128), quantified as urinary ratios of the glucuronide over unconjugated compound for nicotine, cotinine, trans-3'-hydroxycotinine, and 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanol (NNAL), was investigated in regression models assuming a dominant effect of variant alleles.

RESULTS

Correcting for multiple testing, three UGT2B10 variants were associated with cotinine glucuronidation, rs2331559 and rs11726322 in European Americans and rs835309 in African Americans (P ≤ 0.0002). Additional variants predominantly in UGT2B10 were nominally associated with nicotine (P = 0.008-0.04) and cotinine (P = <0.001-0.02) glucuronidation in both ethnicities in addition to UGT2B10*2 in European Americans (P = 0.01, P < 0.001). UGT2B17*2 (P = 0.03) in European Americans and UGT2B7 variants (P = 0.02-0.04) in African Americans were nominally associated with 3HC glucuronidation. UGT1A (P = 0.007-0.01), UGT2B10 (P = 0.02), and UGT2B7 (P = 0.02-0.03) variants in African Americans were nominally associated with NNAL glucuronidation.

CONCLUSIONS

Findings from this initial in vivo study support a role for multiple UGTs in the glucuronidation of tobacco-related compounds in vivo, in particular UGT2B10 and cotinine glucuronidation.

IMPACT

Findings also provide insight into ethnic differences in glucuronidation activity, which could be contributing to ethnic disparities in the risk for smoking-related cancers. Cancer Epidemiol Biomarkers Prev; 24(1); 94-104. ©2014 AACR.

Consensus Genotyper for Exome Sequencing (CGES): improving the quality of exome variant genotypes

MOTIVATION

The development of cost-effective next-generation sequencing methods has spurred the development of high-throughput bioinformatics tools for detection of sequence variation. With many disparate variant-calling algorithms available, investigators must ask, 'Which method is best for my data?' Machine learning research has shown that so-called ensemble methods that combine the output of multiple models can dramatically improve classifier performance. Here we describe a novel variant-calling approach based on an ensemble of variant-calling algorithms, which we term the Consensus Genotyper for Exome Sequencing (CGES). CGES uses a two-stage voting scheme among four algorithm implementations. While our ensemble method can accept variants generated by any variant-calling algorithm, we used GATK2.8, SAMtools, FreeBayes and Atlas-SNP2 in building CGES because of their performance, widespread adoption and diverse but complementary algorithms.

RESULTS

We apply CGES to 132 samples sequenced at the Hudson Alpha Institute for Biotechnology (HAIB, Huntsville, AL) using the Nimblegen Exome Capture and Illumina sequencing technology. Our sample set consisted of 40 complete trios, two families of four, one parent-child duo and two unrelated individuals. CGES yielded the fewest total variant calls (N(CGES) = 139° 897), the highest Ts/Tv ratio (3.02), the lowest Mendelian error rate across all genotypes (0.028%), the highest rediscovery rate from the Exome Variant Server (EVS; 89.3%) and 1000 Genomes (1KG; 84.1%) and the highest positive predictive value (PPV; 96.1%) for a random sample of previously validated de novo variants. We describe these and other quality control (QC) metrics from consensus data and explain how the CGES pipeline can be used to generate call sets of varying quality stringency, including consensus calls present across all four algorithms, calls that are consistent across any three out of four algorithms, calls that are consistent across any two out of four algorithms or a more liberal set of all calls made by any algorithm.

AVAILABILITY AND IMPLEMENTATION

To enable accessible, efficient and reproducible analysis, we implement CGES both as a stand-alone command line tool available for download in GitHub and as a set of Galaxy tools and workflows configured to execute on parallel computers.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Saccadic eye movement abnormalities in autism spectrum disorder indicate dysfunctions in cerebellum and brainstem

Background

Individuals with autism spectrum disorder (ASD) show atypical scan paths during social interaction and when viewing faces, and recent evidence suggests that they also show abnormal saccadic eye movement dynamics and accuracy when viewing less complex and non-social stimuli. Eye movements are a uniquely promising target for studies of ASD as their spatial and temporal characteristics can be measured precisely and the brain circuits supporting them are well-defined. Control of saccade metrics is supported by discrete circuits within the cerebellum and brainstem - two brain regions implicated in magnetic resonance (MR) morphometry and histopathological studies of ASD. The functional integrity of these distinct brain systems can be examined by evaluating different parameters of visually-guided saccades.

Methods

A total of 65 participants with ASD and 43 healthy controls, matched on age (between 6 and 44-years-old), gender and nonverbal IQ made saccades to peripheral targets. To examine the influence of attentional processes, blocked gap and overlap trials were presented. We examined saccade latency, accuracy and dynamics, as well as the trial-to-trial variability of participants’ performance.

Results

Saccades of individuals with ASD were characterized by reduced accuracy, elevated variability in accuracy across trials, and reduced peak velocity and prolonged duration. In addition, their saccades took longer to accelerate to peak velocity, with no alteration in the duration of saccade deceleration. Gap/overlap effects on saccade latencies were similar across groups, suggesting that visual orienting and attention systems are relatively spared in ASD. Age-related changes did not differ across groups.

Conclusions

Deficits precisely and consistently directing eye movements suggest impairment in the error-reducing function of the cerebellum in ASD. Atypical increases in the duration of movement acceleration combined with lower peak saccade velocities implicate pontine nuclei, specifically suggesting reduced excitatory activity in burst cells that drive saccades relative to inhibitory activity in omnipause cells that maintain stable fixation. Thus, our findings suggest that both cerebellar and brainstem abnormalities contribute to altered sensorimotor control in ASD.

Electronic supplementary material

The online version of this article (doi:10.1186/2040-2392-5-47) contains supplementary material, which is available to authorized users.

A framework for the interpretation of de novo mutation in human disease

Spontaneously arising (‘de novo’) mutations play an important role in medical genetics. For diseases with extensive locus heterogeneity – such as autism spectrum disorders (ASDs) – the signal from de novo mutations (DNMs) is distributed across many genes, making it difficult to distinguish disease-relevant mutations from background variation. We provide a statistical framework for the analysis of DNM excesses per gene and gene set by calibrating a model of de novo mutation. We applied this framework to DNMs collected from 1,078 ASD trios and – while affirming a significant role for loss-of-function (LoF) mutations – found no excess of de novo LoF mutations in cases with IQ above 100, suggesting that the role of DNMs in ASD may reside in fundamental neurodevelopmental processes. We also used our model to identify ~1,000 genes that are significantly lacking functional coding variation in non-ASD samples and are enriched for de novo LoF mutations identified in ASD cases.

Copy number variation in obsessive-compulsive disorder and tourette syndrome: a cross-disorder study

OBJECTIVE

Obsessive-compulsive disorder (OCD) and Tourette syndrome (TS) are heritable neurodevelopmental disorders with a partially shared genetic etiology. This study represents the first genome-wide investigation of large (>500 kb), rare (<1%) copy number variants (CNVs) in OCD and the largest genome-wide CNV analysis in TS to date.

METHOD

The primary analyses used a cross-disorder design for 2,699 case patients (1,613 ascertained for OCD, 1,086 ascertained for TS) and 1,789 controls. Parental data facilitated a de novo analysis in 348 OCD trios.

RESULTS

Although no global CNV burden was detected in the cross-disorder analysis or in secondary, disease-specific analyses, there was a 3.3-fold increased burden of large deletions previously associated with other neurodevelopmental disorders (p = .09). Half of these neurodevelopmental deletions were located in a single locus, 16p13.11 (5 case patient deletions: 0 control deletions, p = .08 in the current study, p = .025 compared to published controls). Three 16p13.11 deletions were confirmed de novo, providing further support for the etiological significance of this region. The overall OCD de novo rate was 1.4%, which is intermediate between published rates in controls (0.7%) and in individuals with autism or schizophrenia (2-4%).

CONCLUSION

Several converging lines of evidence implicate 16p13.11 deletions in OCD, with weaker evidence for a role in TS. The trend toward increased overall neurodevelopmental CNV burden in TS and OCD suggests that deletions previously associated with other neurodevelopmental disorders may also contribute to these phenotypes.

Dopamine transporter genotype and stimulant dose-response in youth with attention-deficit/hyperactivity disorder

OBJECTIVES

This study seeks to determine if variation in the dopamine transporter gene (SLC6A3/DAT1) moderates the dose-response effects of long-acting dexmethylphenidate (D-MPH) and mixed amphetamine salts (MAS) in children with attention-deficit/hyperactivity disorder (ADHD).

METHODS

Fifty-six children and adolescents (mean age=11.7±2.2) participated in a double-blind, two period crossover, dose-response study with a randomized placebo week in each 4 week drug period. Each period consisted of sequential week-long exposures to three dose levels (10, 20, 25-30 mg, depending upon weight) of D-MPH or MAS.

RESULTS

Doses of 10-20 mg of either D-MPH or MAS had little to no effect on hyperactivity-impulsivity and total ADHD symptom scores in subjects with the 9/9 genotype; this was in contrast to the dose-response curves of subjects with either the 10/10 or 10/9 genotype.

CONCLUSIONS

ADHD youth with the 9/9 genotype may require higher stimulant doses to achieve adequate symptom control.

Modest impact on risk for autism spectrum disorder of rare copy number variants at 15q11.2, specifically breakpoints 1 to 2

The proximal region of chromosome 15 is one of the genomic hotspots for copy number variants (CNVs). Among the rearrangements observed in this region, CNVs from the interval between the common breakpoints 1 and 2 (BP1 and BP2) have been reported cosegregating with autism spectrum disorder (ASD). Although evidence supporting an association between BP1-BP2 CNVs and autism accumulates, the magnitude of the effect of BP1-BP2 CNVs remains elusive, posing a great challenge to recurrence-risk counseling. To gain further insight into their pathogenicity for ASD, we estimated the penetrance of the BP1-BP2 CNVs for ASD as well as their effects on ASD-related phenotypes in a well-characterized ASD sample (n = 2525 families). Transmission disequilibrium test revealed significant preferential transmission only for the duplicated chromosome in probands (20T:9NT). The penetrance of the BP1-BP2 CNVs for ASD was low, conferring additional risks of 0.3% (deletion) and 0.8% (duplication). Stepwise regression analyses suggest a greater effect of the CNVs on ASD-related phenotype in males and when maternally inherited. Taken together, the results are consistent with BP1-BP2 CNVs as risk factors for autism. However, their effect is modest, more akin to that seen for common variants. To be consistent with the current American College of Medical Genetics guidelines for interpretation of postnatal CNV, the BP1-BP2 deletion and duplication CNVs would probably best be classified as variants of uncertain significance (VOUS): they appear to have an impact on risk, but one so modest that these CNVs do not merit pathogenic status.

The Autism Simplex Collection: an international, expertly phenotyped autism sample for genetic and phenotypic analyses

BACKGROUND

There is an urgent need for expanding and enhancing autism spectrum disorder (ASD) samples, in order to better understand causes of ASD.

METHODS

In a unique public-private partnership, 13 sites with extensive experience in both the assessment and diagnosis of ASD embarked on an ambitious, 2-year program to collect samples for genetic and phenotypic research and begin analyses on these samples. The program was called The Autism Simplex Collection (TASC). TASC sample collection began in 2008 and was completed in 2010, and included nine sites from North America and four sites from Western Europe, as well as a centralized Data Coordinating Center.

RESULTS

Over 1,700 trios are part of this collection, with DNA from transformed cells now available through the National Institute of Mental Health (NIMH). Autism Diagnostic Interview-Revised (ADI-R) and Autism Diagnostic Observation Schedule-Generic (ADOS-G) measures are available for all probands, as are standardized IQ measures, Vineland Adaptive Behavioral Scales (VABS), the Social Responsiveness Scale (SRS), Peabody Picture Vocabulary Test (PPVT), and physical measures (height, weight, and head circumference). At almost every site, additional phenotypic measures were collected, including the Broad Autism Phenotype Questionnaire (BAPQ) and Repetitive Behavior Scale-Revised (RBS-R), as well as the non-word repetition scale, Communication Checklist (Children's or Adult), and Aberrant Behavior Checklist (ABC). Moreover, for nearly 1,000 trios, the Autism Genome Project Consortium (AGP) has carried out Illumina 1 M SNP genotyping and called copy number variation (CNV) in the samples, with data being made available through the National Institutes of Health (NIH). Whole exome sequencing (WES) has been carried out in over 500 probands, together with ancestry matched controls, and this data is also available through the NIH. Additional WES is being carried out by the Autism Sequencing Consortium (ASC), where the focus is on sequencing complete trios. ASC sequencing for the first 1,000 samples (all from whole-blood DNA) is complete and data will be released in 2014. Data is being made available through NIH databases (database of Genotypes and Phenotypes (dbGaP) and National Database for Autism Research (NDAR)) with DNA released in Dist 11.0. Primary funding for the collection, genotyping, sequencing and distribution of TASC samples was provided by Autism Speaks and the NIH, including the National Institute of Mental Health (NIMH) and the National Human Genetics Research Institute (NHGRI).

CONCLUSIONS

TASC represents an important sample set that leverages expert sites. Similar approaches, leveraging expert sites and ongoing studies, represent an important path towards further enhancing available ASD samples.

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⁻⁵) 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⁻¹⁵, ∼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.

Measuring anxiety as a treatment endpoint in youth with autism spectrum disorder

Despite the high rate of anxiety in individuals with autism spectrum disorder (ASD), measuring anxiety in ASD is fraught with uncertainty. This is due, in part, to incomplete consensus on the manifestations of anxiety in this population. Autism Speaks assembled a panel of experts to conduct a systematic review of available measures for anxiety in youth with ASD. To complete the review, the panel held monthly conference calls and two face-to-face meetings over a fourteen-month period. Thirty eight published studies were reviewed and ten assessment measures were examined: four were deemed appropriate for use in clinical trials, although with conditions; three were judged to be potentially appropriate, while three were considered not useful for clinical trials assessing anxiety. Despite recent advances, additional relevant, reliable and valid outcome measures are needed to evaluate treatments for anxiety in ASD.

A deletion involving CD38 and BST1 results in a fusion transcript in a patient with autism and asthma

CD38 encodes a ligand in the oxytocin signaling pathway. Some single nucleotide polymorphisms in this gene have been associated with low serum oxytocin levels in autism spectrum disorder (ASD) patients. Oxytocin disruption has been hypothesized to account for features of ASD, including impaired communication and social behavior, based on animal studies. Recent human studies have shown administration of oxytocin improving emotion recognition, promoting social behavior, and improving auditory processing of social stimuli in ASD patients. In addition to its role in oxytocin signaling, CD38 is involved in the regulation of calcium concentration in airway smooth muscle with impairment of CD38 being implicated in airway diseases like asthma. While a number of studies have implicated rare chromosomal deletions and duplications in helping determine genetic risk for autism, there are to our knowledge no reports describing rearrangements involving CD38 or deletions in patients with ASD. Here, we present two sisters diagnosed with autism and with features of regression-previously acquired speech lost in the second year of life. The younger sister, who also had asthma, inherited a maternal deletion of 4p15.32 that results in a BST1-CD38 fusion transcript. Their mother's deletion was mosaic and she was not affected. Although further work is required to assess functional consequences of the fusion transcript, we hypothesize that the proband's deletion may have served as a risk factor for autism that, when combined with other susceptibility variants, resulted in a more severe presentation than her sister.

Environmental and state-level regulatory factors affect the incidence of autism and intellectual disability

Many factors affect the risks for neurodevelopmental maladies such as autism spectrum disorders (ASD) and intellectual disability (ID). To compare environmental, phenotypic, socioeconomic and state-policy factors in a unified geospatial framework, we analyzed the spatial incidence patterns of ASD and ID using an insurance claims dataset covering nearly one third of the US population. Following epidemiologic evidence, we used the rate of congenital malformations of the reproductive system as a surrogate for environmental exposure of parents to unmeasured developmental risk factors, including toxins. Adjusted for gender, ethnic, socioeconomic, and geopolitical factors, the ASD incidence rates were strongly linked to population-normalized rates of congenital malformations of the reproductive system in males (an increase in ASD incidence by 283% for every percent increase in incidence of malformations, 95% CI: [91%, 576%], p<6×10(-5)). Such congenital malformations were barely significant for ID (94% increase, 95% CI: [1%, 250%], p = 0.0384). Other congenital malformations in males (excluding those affecting the reproductive system) appeared to significantly affect both phenotypes: 31.8% ASD rate increase (CI: [12%, 52%], p<6×10(-5)), and 43% ID rate increase (CI: [23%, 67%], p<6×10(-5)). Furthermore, the state-mandated rigor of diagnosis of ASD by a pediatrician or clinician for consideration in the special education system was predictive of a considerable decrease in ASD and ID incidence rates (98.6%, CI: [28%, 99.99%], p = 0.02475 and 99% CI: [68%, 99.99%], p = 0.00637 respectively). Thus, the observed spatial variability of both ID and ASD rates is associated with environmental and state-level regulatory factors; the magnitude of influence of compound environmental predictors was approximately three times greater than that of state-level incentives. The estimated county-level random effects exhibited marked spatial clustering, strongly indicating existence of as yet unidentified localized factors driving apparent disease incidence. Finally, we found that the rates of ASD and ID at the county level were weakly but significantly correlated (Pearson product-moment correlation 0.0589, p = 0.00101), while for females the correlation was much stronger (0.197, p<2.26×10(-16)).

A survey of seizures and current treatments in 15q duplication syndrome

OBJECTIVE

Seizures are common in individuals with duplications of chromosome 15q11.2-q13 (Dup15q). The goal of this study was to examine the phenotypes and treatments of seizures in Dup15q in a large population.

METHODS

A detailed electronic survey was conducted through the Dup15q Alliance containing comprehensive questions regarding seizures and their treatments in Dup15q.

RESULTS

There were 95 responses from Dup15q families. For the 83 with idic(15), 63% were reported to have seizures, of which 81% had multiple seizure types and 42% had infantile spasms. Other common seizure types were tonic-clonic, atonic, myoclonic, and focal. Only 3 of 12 individuals with int dup(15) had seizures. Broad spectrum antiepileptic drugs (AEDs) were the most effective medications, but carbamazepine and oxcarbazepine were also effective, although typical benzodiazepines were relatively ineffective. There was a 24% response rate (>90% seizure reduction) to the first AED tried. For those with infantile spasms, adrenocorticotropic hormone (ACTH) was more effective than vigabatrin.

SIGNIFICANCE

This is the largest study assessing seizures in Duplication 15q syndrome, but because this was a questionnaire-based study with a low return rate, it is susceptible to bias. Seizures are common in idic(15) and typically difficult to control, often presenting with infantile spasms and progressing to a Lennox-Gastaut-type syndrome. Seizures in those with int dup(15) are less common, with a frequency similar to the general autism population. In addition to broad spectrum AED, medications such as carbamazepine and oxcarbazepine are also relatively effective in controlling seizures in this population, suggesting a possible multifocal etiology, which may also explain the high rate of infantile spasms. Our small sample suggests a relative lack of efficacy of vigabatrin and other γ-aminobutyric acid (GABA)ergic medications, such as typical benzodiazepines, which may be attributable to abnormal GABAergic transmission resulting from the duplication of a cluster of GABAβ3 receptor genes in the 15q11.2-13 region.

Aggression in children with autism spectrum disorders and a clinic-referred comparison group

A gap exists in the literature regarding aggression in autism spectrum disorders and how this behavior compares to other groups. In this multisite study, the Children's Scale for Hostility and Aggression: Reactive/Proactive and the Aggression subscale of the Child Behavior Checklist were rated for 414 children with autism spectrum disorder (autistic disorder, 69%; pervasive developmental disorder not otherwise specified, 24%; Asperger's disorder, 7%) and 243 clinic-referred children without autism spectrum disorder, aged 1-21 years (mean age about 7 years). Participants were not selected for aggressive behavior. Relative to the comparison group, children with autism spectrum disorder were reported to have less aggression and were more likely to be rated as reactive rather than proactive. Among all subjects, sex was not associated with aggression; higher IQ/adaptive behavior and older age were associated with more sophisticated types of aggression, while lower scores on IQ, adaptive behavior, and communication measures were associated with more physical aggression. The interaction between demographic variables and diagnosis was significant only for age: younger but not older children with autism spectrum disorder showed less aggression than clinic-referred controls.