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.

Examining the latent structure and correlates of sensory reactivity in autism: a multi-site integrative data analysis by the autism sensory research consortium

BACKGROUND

Differences in responding to sensory stimuli, including sensory hyperreactivity (HYPER), hyporeactivity (HYPO), and sensory seeking (SEEK) have been observed in autistic individuals across sensory modalities, but few studies have examined the structure of these "supra-modal" traits in the autistic population.

METHODS

Leveraging a combined sample of 3868 autistic youth drawn from 12 distinct data sources (ages 3-18 years and representing the full range of cognitive ability), the current study used modern psychometric and meta-analytic techniques to interrogate the latent structure and correlates of caregiver-reported HYPER, HYPO, and SEEK within and across sensory modalities. Bifactor statistical indices were used to both evaluate the strength of a "general response pattern" factor for each supra-modal construct and determine the added value of "modality-specific response pattern" scores (e.g., Visual HYPER). Bayesian random-effects integrative data analysis models were used to examine the clinical and demographic correlates of all interpretable HYPER, HYPO, and SEEK (sub)constructs.

RESULTS

All modality-specific HYPER subconstructs could be reliably and validly measured, whereas certain modality-specific HYPO and SEEK subconstructs were psychometrically inadequate when measured using existing items. Bifactor analyses supported the validity of a supra-modal HYPER construct (ωH = .800) but not a supra-modal HYPO construct (ωH = .653), and supra-modal SEEK models suggested a more limited version of the construct that excluded some sensory modalities (ωH = .800; 4/7 modalities). Modality-specific subscales demonstrated significant added value for all response patterns. Meta-analytic correlations varied by construct, although sensory features tended to correlate most with other domains of core autism features and co-occurring psychiatric symptoms (with general HYPER and speech HYPO demonstrating the largest numbers of practically significant correlations).

LIMITATIONS

Conclusions may not be generalizable beyond the specific pool of items used in the current study, which was limited to caregiver report of observable behaviors and excluded multisensory items that reflect many "real-world" sensory experiences.

CONCLUSION

Of the three sensory response patterns, only HYPER demonstrated sufficient evidence for valid interpretation at the supra-modal level, whereas supra-modal HYPO/SEEK constructs demonstrated substantial psychometric limitations. For clinicians and researchers seeking to characterize sensory reactivity in autism, modality-specific response pattern scores may represent viable alternatives that overcome many of these limitations.

Examining the Latent Structure and Correlates of Sensory Reactivity in Autism: A Multi-site Integrative Data Analysis by the Autism Sensory Research Consortium

Background Differences in responding to sensory stimuli, including sensory hyperreactivity (HYPER), hyporeactivity (HYPO), and sensory seeking (SEEK) have been observed in autistic individuals across sensory modalities, but few studies have examined the structure of these "supra-modal" traits in the autistic population. Methods Leveraging a combined sample of 3,868 autistic youth drawn from 12 distinct data sources (ages 3-18 years and representing the full range of cognitive ability), the current study used modern psychometric and meta-analytic techniques to interrogate the latent structure and correlates of caregiver-reported HYPER, HYPO, and SEEK within and across sensory modalities. Bifactor statistical indices were used to both evaluate the strength of a "general response pattern" factor for each supra-modal construct and determine the added value of "modality-specific response pattern" scores (e.g., Visual HYPER). Bayesian random-effects integrative data analysis models were used to examine the clinical and demographic correlates of all interpretable HYPER, HYPO and SEEK (sub)constructs. Results All modality-specific HYPER subconstructs could be reliably and validly measured, whereas certain modality-specific HYPO and SEEK subconstructs were psychometrically inadequate when measured using existing items. Bifactor analyses unambiguously supported the validity of a supra-modal HYPER construct (ω H  = .800), whereas a coherent supra-modal HYPO construct was not supported (ω H  = .611), and supra-modal SEEK models suggested a more limited version of the construct that excluded some sensory modalities (ω H  = .799; 4/7 modalities). Within each sensory construct, modality-specific subscales demonstrated substantial added value beyond the supra-modal score. Meta-analytic correlations varied by construct, although sensory features tended to correlate most strongly with other domains of core autism features and co-occurring psychiatric symptoms. Certain subconstructs within the HYPO and SEEK domains were also associated with lower adaptive behavior scores. Limitations: Conclusions may not be generalizable beyond the specific pool of items used in the current study, which was limited to parent-report of observable behaviors and excluded multisensory items that reflect many "real-world" sensory experiences. Conclusion Psychometric issues may limit the degree to which some measures of supra-modal HYPO/SEEK can be interpreted. Depending on the research question at hand, modality-specific response pattern scores may represent a valid alternative method of characterizing sensory reactivity in autism.

A phenome-wide association study of polygenic scores for attention deficit hyperactivity disorder across two genetic ancestries in electronic health record data

Testing the association between genetic scores for Attention Deficit Hyperactivity Disorder (ADHD) and health conditions, can help us better understand its complex etiology. Electronic health records linked to genetic data provide an opportunity to test whether genetic scores for ADHD correlate with ADHD and additional health outcomes in a health care context across different age groups. We generated polygenic scores (ADHD-PGS) trained on summary statistics from the latest genome-wide association study of ADHD (N = 55,374) and applied them to genome-wide data from 12,383 unrelated individuals of African-American ancestry and 66,378 unrelated individuals of European ancestry from the Vanderbilt Biobank. Overall, only Tobacco use disorder (TUD) was associated with ADHD-PGS in the African-American ancestry group (Odds ratio [95% confidence intervals] = 1.23[1.16-1.31], p = 9.3 × 10-09 ). Eighty-six phenotypes were associated with ADHD-PGS in the European ancestry individuals, including ADHD (OR[95%CIs] = 1.22[1.16-1.29], p = 3.6 × 10-10 ), and TUD (OR[95%CIs] = 1.22[1.19-1.25], p = 2.8 × 10-46 ). We then stratified outcomes by age (ages 0-11, 12-18, 19-25, 26-40, 41-60, and 61-100). Our results suggest that ADHD polygenic scores are associated with ADHD diagnoses early in life and with an increasing number of health conditions throughout the lifespan (even in the absence of ADHD diagnosis). This study reinforces the utility of applying trait-specific PGSs to biobank data, and performing exploratory sensitivity analyses, to probe relationships among clinical conditions.

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.

Calculating genetic risk for dysfunction in pleiotropic biological processes using whole exome sequencing data

BACKGROUND

Numerous genes are implicated in autism spectrum disorder (ASD). ASD encompasses a wide-range and severity of symptoms and co-occurring conditions; however, the details of how genetic variation contributes to phenotypic differences are unclear. This creates a challenge for translating genetic evidence into clinically useful knowledge. Sleep disturbances are particularly prevalent co-occurring conditions in ASD, and genetics may inform treatment. Identifying convergent mechanisms with evidence for dysfunction that connect ASD and sleep biology could help identify better treatments for sleep disturbances in these individuals.

METHODS

To identify mechanisms that influence risk for ASD and co-occurring sleep disturbances, we analyzed whole exome sequence data from individuals in the Simons Simplex Collection (n = 2380). We predicted protein damaging variants (PDVs) in genes currently implicated in either ASD or sleep duration in typically developing children. We predicted a network of ASD-related proteins with direct evidence for interaction with sleep duration-related proteins encoded by genes with PDVs. Overrepresentation analyses of Gene Ontology-defined biological processes were conducted on the resulting gene set. We calculated the likelihood of dysfunction in the top overrepresented biological process. We then tested if scores reflecting genetic dysfunction in the process were associated with parent-reported sleep duration.

RESULTS

There were 29 genes with PDVs in the ASD dataset where variation was reported in the literature to be associated with both ASD and sleep duration. A network of 108 proteins encoded by ASD and sleep duration candidate genes with PDVs was identified. The mechanism overrepresented in PDV-containing genes that encode proteins in the interaction network with the most evidence for dysfunction was cerebral cortex development (GO:0,021,987). Scores reflecting dysfunction in this process were associated with sleep durations; the largest effects were observed in adolescents (p = 4.65 × 10^-3).

CONCLUSIONS

Our bioinformatic-driven approach detected a biological process enriched for genes encoding a protein-protein interaction network linking ASD gene products with sleep duration gene products where accumulation of potentially damaging variants in individuals with ASD was associated with sleep duration as reported by the parents. Specifically, genetic dysfunction impacting development of the cerebral cortex may affect sleep by disrupting sleep homeostasis which is evidenced to be regulated by this brain region. Future functional assessments and objective measurements of sleep in adolescents with ASD could provide the basis for more informed treatment of sleep problems in these individuals.

Characterizing sleep disorders in an autism-specific collection of electronic health records

OBJECTIVE/BACKGROUND

Sleep problems are common in people on the autism spectrum. This study reviews one detailed approach to querying the electronic health record (EHR) in a large tertiary care center.

PATIENTS/METHODS

We developed methods for identifying people on the autism spectrum and defined their sleep problems using the key words, "sleep" or "melatonin", or International Classification of Diseases (ICD) codes. We examined treatment responses of these individuals to melatonin supplementation.

RESULTS

Sleep problems were documented in 86% of patients with ages ranging from 6 to 30 years old. Our specific keyword search yielded more patients with sleep diagnoses than ICD codes alone. About two-thirds of patients who received melatonin supplementation reported benefit from its use.

CONCLUSIONS

Our study provides a framework for using deidentified medical records to characterize sleep, a common co-occurring condition, in people on the autism spectrum. Using specific keywords could be helpful in future work that queries the EHR.

Psychometric validation and refinement of the Interoception Sensory Questionnaire (ISQ) in adolescents and adults on the autism spectrum

BACKGROUND

Individuals on the autism spectrum are reported to display alterations in interoception, the sense of the internal state of the body. The Interoception Sensory Questionnaire (ISQ) is a 20-item self-report measure of interoception specifically intended to measure this construct in autistic people. The psychometrics of the ISQ, however, have not previously been evaluated in a large sample of autistic individuals.

METHODS

Using confirmatory factor analysis, we evaluated the latent structure of the ISQ in a large online sample of adults on the autism spectrum and found that the unidimensional model fit the data poorly. Using misspecification analysis to identify areas of local misfit and item response theory to investigate the appropriateness of the seven-point response scale, we removed redundant items and collapsed the response options to put forth a novel eight-item, five-response choice ISQ.

RESULTS

The revised, five-response choice ISQ (ISQ-8) showed much improved fit while maintaining high internal reliability. Differential item functioning (DIF) analyses indicated that the items of the ISQ-8 were answered in comparable ways by autistic adolescents and adults and across multiple other sociodemographic groups.

LIMITATIONS

Our results were limited by the fact that we did not collect data for typically developing controls, preventing the analysis of DIF by diagnostic status. Additionally, while this study proposes a new 5-response scale for the ISQ-8, our data were not collected using this method; thus, the psychometric properties for the revised version of this instrument require further investigation.

CONCLUSION

The ISQ-8 shows promise as a reliable and valid measure of interoception in adolescents and adults on the autism spectrum, but additional work is needed to examine its psychometrics in this population. A free online score calculator has been created to facilitate the use of ISQ-8 latent trait scores for further studies of autistic adolescents and adults (available at https://asdmeasures.shinyapps.io/ISQ_score/ ).

An Automated Functional Annotation Pipeline That Rapidly Prioritizes Clinically Relevant Genes for Autism Spectrum Disorder

Human genetic studies have implicated more than a hundred genes in Autism Spectrum Disorder (ASD). Understanding how variation in implicated genes influence expression of co-occurring conditions and drug response can inform more effective, personalized approaches for treatment of individuals with ASD. Rapidly translating this information into the clinic requires efficient algorithms to sort through the myriad of genes implicated by rare gene-damaging single nucleotide and copy number variants, and common variation detected in genome-wide association studies (GWAS). To pinpoint genes that are more likely to have clinically relevant variants, we developed a functional annotation pipeline. We defined clinical relevance in this project as any ASD associated gene with evidence indicating a patient may have a complex, co-occurring condition that requires direct intervention (e.g., sleep and gastrointestinal disturbances, attention deficit hyperactivity, anxiety, seizures, depression), or is relevant to drug development and/or approaches to maximizing efficacy and minimizing adverse events (i.e., pharmacogenomics). Starting with a list of all candidate genes implicated in all manifestations of ASD (i.e., idiopathic and syndromic), this pipeline uses databases that represent multiple lines of evidence to identify genes: (1) expressed in the human brain, (2) involved in ASD-relevant biological processes and resulting in analogous phenotypes in mice, (3) whose products are targeted by approved pharmaceutical compounds or possessing pharmacogenetic variation and (4) whose products directly interact with those of genes with variants recommended to be tested for by the American College of Medical Genetics (ACMG). Compared with 1000 gene sets, each with a random selection of human protein coding genes, more genes in the ASD set were annotated for each category evaluated (p ≤ 1.99 × 10⁻²). Of the 956 ASD-implicated genes in the full set, 18 were flagged based on evidence in all categories. Fewer genes from randomly drawn sets were annotated in all categories (x = 8.02, sd = 2.56, p = 7.75 × 10⁻⁴). Notably, none of the prioritized genes are represented among the 59 genes compiled by the ACMG, and 78% had a pathogenic or likely pathogenic variant in ClinVar. Results from this work should rapidly prioritize potentially actionable results from genetic studies and, in turn, inform future work toward clinical decision support for personalized care based on genetic testing.

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.

The Drosophila Gene Sulfateless Modulates Autism-Like Behaviors

Major challenges to identifying genes that contribute to autism spectrum disorder (ASD) risk include the availability of large ASD cohorts, the contribution of many genes overall, and small effect sizes attributable to common gene variants. An alternative approach is to use a model organism to detect alleles that impact ASD-relevant behaviors and ask whether homologous human genes infer ASD risk. Here we utilized the Drosophila genetic reference panel (DGRP) as a tool to probe for perturbation in naturally occurring behaviors in Drosophila melanogaster that are analogous to three behavior domains: impaired social communication, social reciprocity and repetitive behaviors or restricted interests. Using 40 of the available DGRP lines, we identified single nucleotide polymorphisms (SNPs) in or near genes controlling these behavior domains, including ASD gene orthologs (neurexin 4 and neuroligin 2), an intellectual disability (ID) gene homolog (kirre), and a gene encoding a heparan sulfate (HS) modifying enzyme called sulfateless (sfl). SNPs in sfl were associated with all three ASD-like behaviors. Using RNAi knock-down of neuronal sfl expression, we observed significant changes in expressive and receptive communication during mating, decreased grooming behavior, and increased social spacing. These results suggest a role for HS proteoglycan synthesis and/or modification in normal social communication, repetitive behavior, and social interaction in flies. Finally, using the DGRP to directly identify genetic effects relevant to a neuropsychiatric disorder further demonstrates the utility of the Drosophila system in the discovery of genes relevant to human disease.

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.

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.

Shorter sleep duration is associated with social impairment and comorbidities in ASD

Sleep disturbance, particularly insomnia, is common in children with autism spectrum disorders (ASD). Furthermore, disturbed sleep affects core symptoms and other related comorbidities. Understanding the causes and consequences of sleep disturbances in children with ASD is an important step toward mitigating these symptoms. To better understand the connection between sleep duration and ASD severity, we analyzed ASD-related symptoms using the Autism Diagnostic Interview-Revised (ADI-R), Autism Diagnostic Observation Schedule (ADOS), IQ scores, and parent reports of the average amount of time slept per night that were available in the medical histories of 2,714 children with ASD in the Simons Simplex Collection (SSC). The mean (SD) sleep duration was 555 minutes. Sleep duration and severity of core ASD symptoms were negatively correlated, and sleep duration and IQ scores were positively correlated. Regression results indicated that more severe social impairment, primarily a failure to develop peer relationships, is the core symptom most strongly associated with short sleep duration. Furthermore, increased severity for numerous maladaptive behaviors assessed on the Child Behavior Checklist, as well as reports of attention deficit disorder, depressive disorder, and obsessive compulsive disorder were associated with short sleep duration. Severity scores for social/communication impairment and restricted and repetitive behaviors (RRB) were increased, and IQ scores were decreased, for children reported to sleep ≤420 minutes per night (lower 5th percentile) compared to children sleeping ≥660 minutes (upper 95th percentile). Our results indicate that reduced amounts of sleep are related to more severe symptoms in children with ASD. Autism Res 2017. © 2017 International Society for Autism Research, Wiley Periodicals, Inc. Autism Res 2017, 10: 1221-1238. © 2017 International Society for Autism Research, Wiley Periodicals, Inc.

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.

Zn(2+) reverses functional deficits in a de novo dopamine transporter variant associated with autism spectrum disorder

Our laboratory recently characterized a novel autism spectrum disorder (ASD)-associated de novo missense mutation in the human dopamine transporter (hDAT) gene SLC6A3 (hDAT T356M). This hDAT variant exhibits dysfunctional forward and reverse transport properties that may contribute to DA dysfunction in ASD. Here, we report that Zn²⁺ reverses, at least in part, the functional deficits of ASD-associated hDAT variant T356M. These data suggest that the molecular mechanism targeted by Zn²⁺ to restore partial function in hDAT T356M may be a novel therapeutic target to rescue functional deficits in hDAT variants associated with ASD.

Analysis of CHRNA7 rare variants in autism spectrum disorder susceptibility

Chromosome 15q13.3 recurrent microdeletions are causally associated with a wide range of phenotypes, including autism spectrum disorder (ASD), seizures, intellectual disability, and other psychiatric conditions. Whether the reciprocal microduplication is pathogenic is less certain. CHRNA7, encoding for the alpha7 subunit of the neuronal nicotinic acetylcholine receptor, is considered the likely culprit gene in mediating neurological phenotypes in 15q13.3 deletion cases. To assess if CHRNA7 rare variants confer risk to ASD, we performed copy number variant analysis and Sanger sequencing of the CHRNA7 coding sequence in a sample of 135 ASD cases. Sequence variation in this gene remains largely unexplored, given the existence of a fusion gene, CHRFAM7A, which includes a nearly identical partial duplication of CHRNA7. Hence, attempts to sequence coding exons must distinguish between CHRNA7 and CHRFAM7A, making next-generation sequencing approaches unreliable for this purpose. A CHRNA7 microduplication was detected in a patient with autism and moderate cognitive impairment; while no rare damaging variants were identified in the coding region, we detected rare variants in the promoter region, previously described to functionally reduce transcription. This study represents the first sequence variant analysis of CHRNA7 in a sample of idiopathic autism.

Rare autism-associated variants implicate syntaxin 1 (STX1 R26Q) phosphorylation and the dopamine transporter (hDAT R51W) in dopamine neurotransmission and behaviors

Background

Syntaxin 1 (STX1) is a presynaptic plasma membrane protein that coordinates synaptic vesicle fusion. STX1 also regulates the function of neurotransmitter transporters, including the dopamine (DA) transporter (DAT). The DAT is a membrane protein that controls DA homeostasis through the high-affinity re-uptake of synaptically released DA.

Methods

We adopt newly developed animal models and state-of-the-art biophysical techniques to determine the contribution of the identified gene variants to impairments in DA neurotransmission observed in autism spectrum disorder (ASD).

Outcomes

Here, we characterize two independent autism-associated variants in the genes that encode STX1 and the DAT. We demonstrate that each variant dramatically alters DAT function. We identify molecular mechanisms that converge to inhibit reverse transport of DA and DA-associated behaviors. These mechanisms involve decreased phosphorylation of STX1 at Ser14 mediated by casein kinase 2 as well as a reduction in STX1/DAT interaction. These findings point to STX1/DAT interactions and STX1 phosphorylation as key regulators of DA homeostasis.

Interpretation

We determine the molecular identity and the impact of these variants with the intent of defining DA dysfunction and associated behaviors as possible complications of ASD.

•

We report two independent autism-associated variants in syntaxin and the dopamine transporter.

•

The variants in syntaxin and dopamine transporter each impair reverse transport of dopamine.

•

Dysregulation of dopamine neurotransmission may represent a complication of autism spectrum disorder.

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.

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.

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.

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.

Pro32Pro33 mutations in the integrin β3 PSI domain result in αIIbβ3 priming and enhanced adhesion: reversal of the hypercoagulability phenotype by the Src inhibitor SKI-606

The plasma-membrane integrin αIIbβ3 (CD41/CD61, GPIIbIIIa) is a major functional receptor in platelets during clotting. A common isoform of integrin β3, Leu33Pro is associated with enhanced platelet function and increased risk for coronary thrombosis and stroke, although these findings remain controversial. To better understand the molecular mechanisms by which this sequence variation modifies platelet function, we produced transgenic knockin mice expressing a Pro32Pro33 integrin β3. Consistent with reports utilizing human platelets, we found significantly reduced bleeding and clotting times, as well as increased in vivo thrombosis, in Pro32Pro33 homozygous mice. These alterations paralleled increases in platelet attachment and spreading onto fibrinogen resulting from enhanced integrin αIIbβ3 function. Activation with protease-activated receptor 4- activating peptide, the main thrombin signaling receptor in mice, showed no significant difference in activation of Pro32Pro33 mice as compared with controls, suggesting that inside-out signaling remains intact. However, under unstimulated conditions, the Pro32Pro33 mutation led to elevated Src phosphorylation, facilitated by increased talin interactions with the β3 cytoplasmic domain, indicating that the αIIbβ3 intracellular domains are primed for activation while the ligand-binding domain remains unchanged. Acute dosing of animals with a Src inhibitor was sufficient to rescue the clotting phenotype in knockin mice to wild-type levels. Together, our data establish that the Pro32Pro33 structural alteration modifies the function of integrin αIIbβ3, priming the integrin for outside-in signaling, ultimately leading to hypercoagulability. Furthermore, our data may support a novel approach to antiplatelet therapy by Src inhibition where hemostasis is maintained while reducing risk for cardiovascular disease.

Brain and peripheral pharmacokinetics of levodopa in the cynomolgus monkey following administration of opicapone, a third generation nitrocatechol COMT inhibitor

OBJECTIVE

The present study aimed at evaluating the effect of opicapone, a third generation nitrocatechol catechol-O-methyltransferase (COMT) inhibitor, on the systemic and central bioavailability of 3,4-dihydroxy-l-phenylalanine (levodopa) and related metabolites in the cynomolgus monkey.

METHODS

Four monkeys, implanted with guiding cannulas for microdialysis probes, in the substantia nigra, dorsal striatum and prefrontal cortex, were randomized in two groups that received, in a crossover design, vehicle or 100 mg/kg opicapone for 14 days. Twenty-three hours after last administration of vehicle or opicapone, animals were challenged with levodopa/benserazide (12/3 mg/kg). Extracellular dialysate and blood samples were collected over 360 min (at 30 min intervals) for the assays of catecholamine and COMT activity.

RESULTS

Opicapone increased levodopa systemic exposure by 2-fold not changing Cmax values and reduced both 3-O-methyldopa (3-OMD) exposure and Cmax values by 5-fold. These changes were accompanied by ∼76-84% reduction in erythrocyte COMT activity. In dorsal striatum and substantia nigra, opicapone increased levodopa exposure by 1.7- and 1.4-fold, respectively, reducing 3-OMD exposure by 5- and 7-fold respectively. DOPAC exposure was increased by 4-fold in the substantia nigra. In the prefrontal cortex, opicapone increased levodopa exposure and reduced 3-OMD levels by 2.3- and 2.4-fold, respectively.

CONCLUSIONS

Opicapone behaved as long-acting COMT inhibitor that markedly increased systemic and central levodopa bioavailability. Opicapone is a strong candidate to fill the unmet need for COMT inhibitors that lead to more sustained levodopa levels in Parkinson's disease patients.

Adjusting head circumference for covariates in autism: clinical correlates of a highly heritable continuous trait

BACKGROUND

Brain development follows a different trajectory in children with autism spectrum disorders (ASD) than in typically developing children. A proxy for neurodevelopment could be head circumference (HC), but studies assessing HC and its clinical correlates in ASD have been inconsistent. This study investigates HC and clinical correlates in the Simons Simplex Collection cohort.

METHODS

We used a mixed linear model to estimate effects of covariates and the deviation from the expected HC given parental HC (genetic deviation). After excluding individuals with incomplete data, 7225 individuals in 1891 families remained for analysis. We examined the relationship between HC/genetic deviation of HC and clinical parameters.

RESULTS

Gender, age, height, weight, genetic ancestry, and ASD status were significant predictors of HC (estimate of the ASD effect = .2 cm). HC was approximately normally distributed in probands and unaffected relatives, with only a few outliers. Genetic deviation of HC was also normally distributed, consistent with a random sampling of parental genes. Whereas larger HC than expected was associated with ASD symptom severity and regression, IQ decreased with the absolute value of the genetic deviation of HC.

CONCLUSIONS

Measured against expected values derived from covariates of ASD subjects, statistical outliers for HC were uncommon. HC is a strongly heritable trait, and population norms for HC would be far more accurate if covariates including genetic ancestry, height, and age were taken into account. The association of diminishing IQ with absolute deviation from predicted HC values suggests HC could reflect subtle underlying brain development and warrants further investigation.

Genetic background modulates phenotypes of serotonin transporter Ala56 knock-in mice

Background

Previously, we identified multiple, rare serotonin (5-HT) transporter (SERT) variants in children with autism spectrum disorder (ASD). Although in our study the SERT Ala56 variant was over-transmitted to ASD probands, it was also seen in some unaffected individuals, suggesting that associated ASD risk is influenced by the epistatic effects of other genetic variation. Subsequently, we established that mice expressing the SERT Ala56 variant on a 129S6/S4 genetic background display multiple biochemical, physiological and behavioral changes, including hyperserotonemia, altered 5-HT receptor sensitivity, and altered social, communication, and repetitive behavior. Here we explore the effects of genetic background on SERT Ala56 knock-in phenotypes.

Methods

To explore the effects of genetic background, we backcrossed SERT Ala56 mice on the 129 background into a C57BL/6 (B6) background to achieve congenic B6 SERT Ala56 mice, and assessed autism-relevant behavior, including sociability, ultrasonic vocalizations, and repetitive behavior in the home cage, as well as serotonergic phenotypes, including whole blood serotonin levels and serotonin receptor sensitivity.

Results

One consistent phenotype between the two strains was performance in the tube test for dominance, where mutant mice displayed a greater tendency to withdraw from a social encounter in a narrow tube as compared to wildtype littermate controls. On the B6 background, mutant pup ultrasonic vocalizations were significantly increased, in contrast to decreased vocalizations seen previously on the 129 background. Several phenotypes seen on the 129 background were reduced or absent when the mutation was placed on the B6 background, including hyperserotonemia, 5-HT receptor hypersensivity, and repetitive behavior.

Conclusions

Our findings provide a cogent example of how epistatic interactions can modulate the impact of functional genetic variation and suggest that some aspects of social behavior may be especially sensitive to changes in SERT function. Finally, these results provide a platform for the identification of genes that may modulate the risk of ASD in humans.

Integrated model of de novo and inherited genetic variants yields greater power to identify risk genes

De novo mutations affect risk for many diseases and disorders, especially those with early-onset. An example is autism spectrum disorders (ASD). Four recent whole-exome sequencing (WES) studies of ASD families revealed a handful of novel risk genes, based on independent de novo loss-of-function (LoF) mutations falling in the same gene, and found that de novo LoF mutations occurred at a twofold higher rate than expected by chance. However successful these studies were, they used only a small fraction of the data, excluding other types of de novo mutations and inherited rare variants. Moreover, such analyses cannot readily incorporate data from case-control studies. An important research challenge in gene discovery, therefore, is to develop statistical methods that accommodate a broader class of rare variation. We develop methods that can incorporate WES data regarding de novo mutations, inherited variants present, and variants identified within cases and controls. TADA, for Transmission And De novo Association, integrates these data by a gene-based likelihood model involving parameters for allele frequencies and gene-specific penetrances. Inference is based on a Hierarchical Bayes strategy that borrows information across all genes to infer parameters that would be difficult to estimate for individual genes. In addition to theoretical development we validated TADA using realistic simulations mimicking rare, large-effect mutations affecting risk for ASD and show it has dramatically better power than other common methods of analysis. Thus TADA's integration of various kinds of WES data can be a highly effective means of identifying novel risk genes. Indeed, application of TADA to WES data from subjects with ASD and their families, as well as from a study of ASD subjects and controls, revealed several novel and promising ASD candidate genes with strong statistical support.

The genetic underpinnings of autism spectrum disorder (ASD) have proven difficult to determine, despite a wealth of evidence for genetic causes and ongoing effort to identify genes. Recently investigators sequenced the coding regions of the genomes from ASD children along with their unaffected parents (ASD trios) and identified numerous new candidate genes by pinpointing spontaneously occurring (de novo) mutations in the affected offspring. A gene with a severe (de novo) mutation observed in more than one individual is immediately implicated in ASD; however, the majority of severe mutations are observed only once per gene. These genes create a short list of candidates, and our results suggest about 50% are true risk genes. To strengthen our inferences, we develop a novel statistical method (TADA) that utilizes inherited variation transmitted to affected offspring in conjunction with (de novo) mutations to identify risk genes. Through simulations we show that TADA dramatically increases power. We apply this approach to nearly 1000 ASD trios and 2000 subjects from a case-control study and identify several promising genes. Through simulations and application we show that TADA's integration of sequencing data can be a highly effective means of identifying risk genes.

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

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

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

Background

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

Methods

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

Results

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

Conclusions

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

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

Background

Autism spectrum disorder is a severe early onset neurodevelopmental disorder with high heritability but significant heterogeneity. Traditional genome-wide approaches to test for an association of common variants with autism susceptibility risk have met with limited success. However, novel methods to identify moderate risk alleles in attainable sample sizes are now gaining momentum.

Methods

In this study, we utilized publically available genome-wide association study data from the Autism Genome Project and annotated the results (P <0.001) for expression quantitative trait loci present in the parietal lobe (GSE35977), cerebellum (GSE35974) and lymphoblastoid cell lines (GSE7761). We then performed a test of enrichment by comparing these results to simulated data conditioned on minor allele frequency to generate an empirical P-value indicating statistically significant enrichment of expression quantitative trait loci in top results from the autism genome-wide association study.

Results

Our findings show a global enrichment of brain expression quantitative trait loci, but not lymphoblastoid cell line expression quantitative trait loci, among top single nucleotide polymorphisms from an autism genome-wide association study. Additionally, the data implicates individual genes SLC25A12, PANX1 and PANX2 as well as pathways previously implicated in autism.

Conclusions

These findings provide supportive rationale for the use of annotation-based approaches to genome-wide association studies.

A novel approach of homozygous haplotype sharing identifies candidate genes in autism spectrum disorder

Autism spectrum disorder (ASD) is a highly heritable disorder of complex and heterogeneous aetiology. It is primarily characterized by altered cognitive ability including impaired language and communication skills and fundamental deficits in social reciprocity. Despite some notable successes in neuropsychiatric genetics, overall, the high heritability of ASD (~90%) remains poorly explained by common genetic risk variants. However, recent studies suggest that rare genomic variation, in particular copy number variation, may account for a significant proportion of the genetic basis of ASD. We present a large scale analysis to identify candidate genes which may contain low-frequency recessive variation contributing to ASD while taking into account the potential contribution of population differences to the genetic heterogeneity of ASD. Our strategy, homozygous haplotype (HH) mapping, aims to detect homozygous segments of identical haplotype structure that are shared at a higher frequency amongst ASD patients compared to parental controls. The analysis was performed on 1,402 Autism Genome Project trios genotyped for 1 million single nucleotide polymorphisms (SNPs). We identified 25 known and 1,218 novel ASD candidate genes in the discovery analysis including CADM2, ABHD14A, CHRFAM7A, GRIK2, GRM3, EPHA3, FGF10, KCND2, PDZK1, IMMP2L and FOXP2. Furthermore, 10 of the previously reported ASD genes and 300 of the novel candidates identified in the discovery analysis were replicated in an independent sample of 1,182 trios. Our results demonstrate that regions of HH are significantly enriched for previously reported ASD candidate genes and the observed association is independent of gene size (odds ratio 2.10). Our findings highlight the applicability of HH mapping in complex disorders such as ASD and offer an alternative approach to the analysis of genome-wide association data.

A multisite study of the clinical diagnosis of different autism spectrum disorders

CONTEXT

Best-estimate clinical diagnoses of specific autism spectrum disorders (autistic disorder, pervasive developmental disorder-not otherwise specified, and Asperger syndrome) have been used as the diagnostic gold standard, even when information from standardized instruments is available.

OBJECTIVE

To determine whether the relationships between behavioral phenotypes and clinical diagnoses of different autism spectrum disorders vary across 12 university-based sites.

DESIGN

Multisite observational study collecting clinical phenotype data (diagnostic, developmental, and demographic) for genetic research. Classification trees were used to identify characteristics that predicted diagnosis across and within sites.

SETTING

Participants were recruited through 12 university-based autism service providers into a genetic study of autism.

PARTICIPANTS

A total of 2102 probands (1814 male probands) between 4 and 18 years of age (mean [SD] age, 8.93 [3.5] years) who met autism spectrum criteria on the Autism Diagnostic Interview-Revised and the Autism Diagnostic Observation Schedule and who had a clinical diagnosis of an autism spectrum disorder.

MAIN OUTCOME MEASURE

Best-estimate clinical diagnoses predicted by standardized scores from diagnostic, cognitive, and behavioral measures.

RESULTS

Although distributions of scores on standardized measures were similar across sites, significant site differences emerged in best-estimate clinical diagnoses of specific autism spectrum disorders. Relationships between clinical diagnoses and standardized scores, particularly verbal IQ, language level, and core diagnostic features, varied across sites in weighting of information and cutoffs.

CONCLUSIONS

Clinical distinctions among categorical diagnostic subtypes of autism spectrum disorders were not reliable even across sites with well-documented fidelity using standardized diagnostic instruments. Results support the move from existing subgroupings of autism spectrum disorders to dimensional descriptions of core features of social affect and fixated, repetitive behaviors, together with characteristics such as language level and cognitive function.

Accuracy of phenotyping children with autism based on parent report: what specifically do we gain phenotyping "rapidly"?

Autism spectrum disorder (ASD) is considered among the most heritable of all neurodevelopmental and psychiatric disorders, but identification of etiologically significant genetic markers and risk variants has been hampered by a lack of sufficiently large samples. Rapid phenotyping procedures, where self-report measures are used instead of extensive clinical assessment, have been proposed as methods for amassing large genetic databases due to their hypothesized time-efficiency and affordability. We assessed the diagnostic accuracy of potential rapid phenotyping procedures using the Social Communication Questionnaire and the Social Responsiveness Scale in a sample of 333 children who also received extensive phenotypic assessments. While the rapid phenotyping measures were able to accurately identify a large number of children with ASD, they also frequently failed to differentiate children with ASD from children with other complex neurobehavioral profiles. These data support the continued need of expert clinical validation in combination with rapid phenotyping procedures in order to accurately amass large-scale genetic collections of children with ASD.

Use of array CGH to detect exonic copy number variants throughout the genome in autism families detects a novel deletion in TMLHE

Autism is a neurodevelopmental disorder with increasing evidence of heterogeneous genetic etiology including de novo and inherited copy number variants (CNVs). We performed array comparative genomic hybridization using a custom Agilent 1 M oligonucleotide array intended to cover 197 332 unique exons in RefSeq genes; 98% were covered by at least one probe and 95% were covered by three or more probes with the focus on detecting relatively small CNVs that would implicate a single protein-coding gene. The study group included 99 trios from the Simons Simplex Collection. The analysis identified and validated 55 potentially pathogenic CNVs, categorized as de novo autosomal heterozygous, inherited homozygous autosomal, complex autosomal and hemizygous deletions on the X chromosome of probands. Twenty percent (11 of 55) of these CNV calls were rare when compared with the Database of Genomic Variants. Thirty-six percent (20 of 55) of the CNVs were also detected in the same samples in an independent analysis using the 1 M Illumina single-nucleotide polymorphism array. Findings of note included a common and sometimes homozygous 61 bp exonic deletion in SLC38A10, three CNVs found in lymphoblast-derived DNA but not present in whole-blood derived DNA and, most importantly, in a male proband, an exonic deletion of the TMLHE (trimethyllysine hydroxylase epsilon) that encodes the first enzyme in the biosynthesis of carnitine. Data for CNVs present in lymphoblasts but absent in fresh blood DNA suggest that these represent clonal outgrowth of individual B cells with pre-existing somatic mutations rather than artifacts arising in cell culture. GEO accession number GSE23765 (http://www.ncbi.nlm.nih.gov/geo/, date last accessed on 30 August 2011). Genboree accession: http://genboree.org/java-bin/gbrowser.jsp?refSeqId=1868&entryPointId=chr17&from=53496072&to=53694382&isPublic=yes, date last accessed on 30 August 2011.

Colocalization and regulated physical association of presynaptic serotonin transporters with A₃ adenosine receptors

Activation of A₃ adenosine receptors (A₃ARs) rapidly enhances the activity of antidepressant-sensitive serotonin (5-HT) transporters (SERTs) in vitro, ex vivo, and in vivo. A₃AR agonist stimulation of SERT activity is lost in A₃AR knockout mice. A₃AR-stimulated SERT activity is mediated by protein kinase G1 (PKGI)- and p38 mitogen-activated protein kinase (MAPK)-linked pathways that support, respectively, enhanced SERT surface expression and catalytic activation. The mechanisms by which A₃ARs target SERTs among other potential effectors is unknown. Here we present evidence that A₃ARs are coexpressed with SERT in midbrain serotonergic neurons and form a physical complex in A₃AR/hSERT cotransfected cells. Treatment of A₃AR/SERT-cotransfected Chinese hamster ovary cells with the A₃AR agonist N⁶-(3-iodobenzyl)-N-methyl-5'-carbamoyladenosine (1 μM, 10 min), conditions previously reported to increase SERT surface expression and 5-HT uptake activity, enhanced the abundance of A₃AR/SERT complexes in a PKGI-dependent manner. Cotransfection of SERT with L90V-A₃AR, a hyperfunctional coding variant identified in subjects with autism spectrum disorder, resulted in a prolonged recovery of receptor/transporter complexes after A₃AR activation. Because PKGI and nitric-oxide synthetase are required for A₃AR stimulation of SERT activity, and proteins PKGI and NOS both form complexes with SERT, our findings suggest a mechanism by which signaling pathways coordinating A₃AR signaling to SERT can be spatially restricted and regulated, as well as compromised by neuropsychiatric disorders.

Association of oxytocin receptor (OXTR) gene variants with multiple phenotype domains of autism spectrum disorder

Autism spectrum disorder (ASD) is characterized by core deficits in social behavior, communication, and behavioral flexibility. Several lines of evidence indicate that oxytocin, signaling through its receptor (OXTR), is important in a wide range of social behaviors. In attempts to determine whether genetic variations in the oxytocin signaling system contribute to ASD susceptibility, seven recent reports indicated association of common genetic polymorphisms in the OXTR gene with ASD. Each involved relatively small sample sizes (57 to 436 families) and, where it was examined, failed to identify association of OXTR polymorphisms with measures of social behavior in individuals with ASD. We report genetic association analysis of 25 markers spanning the OXTR locus in 1,238 pedigrees including 2,333 individuals with ASD. Association of three markers previously implicated in ASD susceptibility, rs2268493 (P = 0.043), rs1042778 (P = 0.037), and rs7632287 (P = 0.016), was observed. Further, these genetic markers were associated with multiple core ASD phenotypes, including social domain dysfunction, measured by standardized instruments used to diagnose and describe ASD. The data suggest association of OXTR genetic polymorphisms with ASD, although the results should be interpreted with caution because none of the significant associations would survive appropriate correction for multiple comparisons. However, the current findings of association in a large independent cohort are consistent with previous results, and the biological plausibility of participation of the oxytocin signaling system in modulating social disruptions characteristic of ASD, suggest that functional polymorphisms of OXTR may contribute to ASD risk in a subset of families.

Genetic analysis of biological pathway data through genomic randomization

Genome Wide Association Studies (GWAS) are a standard approach for large-scale common variation characterization and for identification of single loci predisposing to disease. However, due to issues of moderate sample sizes and particularly multiple testing correction, many variants of smaller effect size are not detected within a single allele analysis framework. Thus, small main effects and potential epistatic effects are not consistently observed in GWAS using standard analytical approaches that consider only single SNP alleles. Here, we propose unique methodology that aggregates variants of interest (for example, genes in a biological pathway) using GWAS results. Multiple testing and type I error concerns are minimized using empirical genomic randomization to estimate significance. Randomization corrects for common pathway-based analysis biases, such as SNP coverage and density, linkage disequilibrium, gene size and pathway size. Pathway Analysis by Randomization Incorporating Structure (PARIS) applies this randomization and in doing so directly accounts for linkage disequilibrium effects. PARIS is independent of association analysis method and is thus applicable to GWAS datasets of all study designs. Using the KEGG database as an example, we apply PARIS to the publicly available Autism Genetic Resource Exchange GWAS dataset, revealing pathways with a significant enrichment of positive association results.

Parent-of-origin effects of the serotonin transporter gene associated with autism

A promoter-linked insertion/deletion polymorphism of the serotonin transporter gene (SLC6A4) has been implicated in autism spectrum disorders (ASDs) in numerous family based association studies. However, the results of these investigations have been inconsistent in that both the long and short alleles have been shown to be over-transmitted to affected offspring. In order to further elucidate the relationship between the 5-HTTLPR variant and autism risk, we undertook a thorough study of parent-of-origin effects, maternal genotype effects, and offspring genotype effects in a sample of affected offspring from the Autism Genetic Resource Exchange (AGRE). Both the overall autism phenotype and measures of autism behaviors from the Autism Diagnostic Interview-Revised [Lord et al. (1994); J Autism Dev Disord 24(5): 659–685] were considered. We found evidence of over-transmission (risk allele short, P = 0.012), maternal effects (risk allele long, P = 0.035), and parent-of-origin effects (risk allele short from mother, P = 0.018) of the 5-HTTLPR variant in the AGRE sample. Population- and gender-specific effects were also explored as associations may be heterogeneous across populations and sexes. Parent-of-origin effects of the variant were associated with maternally inherited copies of the short allele that resulted in more impaired overall level of language (P = 0.04). Our study was conducted to further investigate the 5-HTTLPR risk variants by identifying allelic associations that may be population-specific, phenotype-specific, or conferred by maternal or parent-of-origin effects. In light of conflicting observations from previous studies, these are just a few of the possible explanations that deserve attention.

Maternal transmission of a rare GABRB3 signal peptide variant is associated with autism

Maternal 15q11-q13 duplication is the most common copy number variant in autism, accounting for ∼1-3% of cases. The 15q11-q13 region is subject to epigenetic regulation, and genomic copy number losses and gains cause genomic disorders in a parent-of-origin-specific manner. One 15q11-q13 locus encodes the GABA(A) receptor β3 subunit gene (GABRB3), which has been implicated by several studies in both autism and absence epilepsy, and the co-morbidity of epilepsy in autism is well established. We report that maternal transmission of a GABRB3 signal peptide variant (P11S), previously implicated in childhood absence epilepsy, is associated with autism. An analysis of wild-type and mutant β3 subunit-containing α1β3γ2 or α3β3γ2 GABA(A) receptors shows reduced whole-cell current and decreased β3 subunit protein on the cell surface due to impaired intracellular β3 subunit processing. We thus provide the first evidence of an association between a specific GABA(A) receptor defect and autism, direct evidence that this defect causes synaptic dysfunction that is autism relevant and the first maternal risk effect in the 15q11-q13 autism duplication region that is linked to a coding variant.

Rare familial 16q21 microdeletions under a linkage peak implicate cadherin 8 (CDH8) in susceptibility to autism and learning disability

Background

Autism spectrum disorder (ASD) is characterised by impairments in social communication and by a pattern of repetitive behaviours, with learning disability (LD) typically seen in up to 70% of cases. A recent study using the PPL statistical framework identified a novel region of genetic linkage on chromosome 16q21 that is limited to ASD families with LD.

Methods

In this study, two families with autism and/or LD are described which harbour rare >1.6 Mb microdeletions located within this linkage region. The deletion breakpoints are mapped at base-pair resolution and segregation analysis is performed using a combination of 1M single nucleotide polymorphism (SNP) technology, array comparative genomic hybridisation (CGH), long-range PCR, and Sanger sequencing. The frequency of similar genomic variants in control subjects is determined through analysis of published SNP array data. Expression of CDH8, the only gene disrupted by these microdeletions, is assessed using reverse transcriptase PCR and in situ hybridisation analysis of 9 week human embryos.

Results

The deletion of chr16: 60 025 584–61 667 839 was transmitted to three of three boys with autism and LD and none of four unaffected siblings, from their unaffected mother. In a second family, an overlapping deletion of chr16: 58 724 527–60 547 472 was transmitted to an individual with severe LD from his father with moderate LD. No copy number variations (CNVs) disrupting CDH8 were observed in 5023 controls. Expression analysis indicates that the two CDH8 isoforms are present in the developing human cortex.

Conclusion

Rare familial 16q21 microdeletions and expression analysis implicate CDH8 in susceptibility to autism and LD.

A genome-wide scan for common alleles affecting risk for autism

Although autism spectrum disorders (ASDs) have a substantial genetic basis, most of the known genetic risk has been traced to rare variants, principally copy number variants (CNVs). To identify common risk variation, the Autism Genome Project (AGP) Consortium genotyped 1558 rigorously defined ASD families for 1 million single-nucleotide polymorphisms (SNPs) and analyzed these SNP genotypes for association with ASD. In one of four primary association analyses, the association signal for marker rs4141463, located within MACROD2, crossed the genome-wide association significance threshold of P < 5 × 10⁻⁸. When a smaller replication sample was analyzed, the risk allele at rs4141463 was again over-transmitted; yet, consistent with the winner's curse, its effect size in the replication sample was much smaller; and, for the combined samples, the association signal barely fell below the P < 5 × 10⁻⁸ threshold. Exploratory analyses of phenotypic subtypes yielded no significant associations after correction for multiple testing. They did, however, yield strong signals within several genes, KIAA0564, PLD5, POU6F2, ST8SIA2 and TAF1C.

Functional impact of global rare copy number variation in autism spectrum disorders

The autism spectrum disorders (ASDs) are a group of conditions characterized by impairments in reciprocal social interaction and communication, and the presence of restricted and repetitive behaviours. Individuals with an ASD vary greatly in cognitive development, which can range from above average to intellectual disability. Although ASDs are known to be highly heritable ( approximately 90%), the underlying genetic determinants are still largely unknown. Here we analysed the genome-wide characteristics of rare (<1% frequency) copy number variation in ASD using dense genotyping arrays. When comparing 996 ASD individuals of European ancestry to 1,287 matched controls, cases were found to carry a higher global burden of rare, genic copy number variants (CNVs) (1.19 fold, P = 0.012), especially so for loci previously implicated in either ASD and/or intellectual disability (1.69 fold, P = 3.4 x 10(-4)). Among the CNVs there were numerous de novo and inherited events, sometimes in combination in a given family, implicating many novel ASD genes such as SHANK2, SYNGAP1, DLGAP2 and the X-linked DDX53-PTCHD1 locus. We also discovered an enrichment of CNVs disrupting functional gene sets involved in cellular proliferation, projection and motility, and GTPase/Ras signalling. Our results reveal many new genetic and functional targets in ASD that may lead to final connected pathways.

Association of MET with social and communication phenotypes in individuals with autism spectrum disorder

Autism is a complex neurodevelopmental disorder diagnosed by impairments in social interaction, communication, and behavioral flexibility. Autism is highly heritable, but it is not known whether a genetic risk factor contributes to all three core domains of the disorder or autism results from the confluence of multiple genetic risk factors for each domain. We and others reported previously association of variants in the gene encoding the MET receptor tyrosine kinase in five independent samples. We further described enriched association of the MET promoter variant rs1858830 C allele in families with co-occurring autism and gastrointestinal conditions. To test the contribution of this functional MET promoter variant to the domains of autism, we analyzed its association with quantitative scores derived from three instruments used to diagnose and describe autism phenotypes: the Autism Diagnostic Interview-Revised (ADI-R), the Autism Diagnostic Observation Schedule (ADOS), and both the parent and the teacher report forms of the Social Responsiveness Scale (SRS). In 748 individuals from 367 families, the transmission of the MET C allele from parent to child was consistently associated with both social and communication phenotypes of autism. Stratification by gastrointestinal conditions revealed a similar pattern of association with both social and communication phenotypes in 242 individuals with autism from 118 families with co-occurring gastrointestinal conditions, but a lack of association with any domain in 181 individuals from 96 families with ASD and no co-occurring gastrointestinal condition. These data indicate that the MET C allele influences at least two of the three domains of the autism triad.

Microduplications of 16p11.2 are associated with schizophrenia

Recurrent microdeletions and microduplications of a 600 kb genomic region of chromosome 16p11.2 have been implicated in childhood-onset developmental disorders1-3. Here we report the strong association of 16p11.2 microduplications with schizophrenia in two large cohorts. In the primary sample, the microduplication was detected in 12/1906 (0.63%) cases and 1/3971 (0.03%) controls (P=1.2×10^-5, OR=25.8). In the replication sample, the microduplication was detected in 9/2645 (0.34%) cases and 1/2420 (0.04%) controls (P=0.022, OR=8.3). For the series combined, microduplication of 16p11.2 was associated with 14.5-fold increased risk of schizophrenia (95% C.I. [3.3, 62]). A meta-analysis of multiple psychiatric disorders showed a significant association of the microduplication with schizophrenia, bipolar disorder and autism. The reciprocal microdeletion was associated only with autism and developmental disorders. Analysis of patient clinical data showed that head circumference was significantly larger in patients with the microdeletion compared with patients with the microduplication (P = 0.0007). Our results suggest that the microduplication of 16p11.2 confers substantial risk for schizophrenia and other psychiatric disorders, whereas the reciprocal microdeletion is associated with contrasting clinical features.