Primary complex motor stereotypies are associated with de novo damaging DNA coding mutations that identify KDM5B as a risk gene.

Motor stereotypies are common in children with autism spectrum disorder (ASD), intellectual disability, or sensory deprivation, as well as in typically developing children ("primary" stereotypies, pCMS). The precise pathophysiological mechanism for motor stereotypies is unknown, although genetic etiologies have been suggested. In this study, we perform whole-exome DNA sequencing in 129 parent-child trios with pCMS and 853 control trios (118 cases and 750 controls after quality control). We report an increased rate of de novo predicted-damaging DNA coding variants in pCMS versus controls, identifying KDM5B as a high-confidence risk gene and estimating 184 genes conferring risk. Genes harboring de novo damaging variants in pCMS probands show significant overlap with those in Tourette syndrome, ASD, and those in ASD probands with high versus low stereotypy scores. An exploratory analysis of these pCMS gene expression patterns finds clustering within the cortex and striatum during early mid-fetal development. Exploratory gene ontology and network analyses highlight functional convergence in calcium ion transport, demethylation, cell signaling, cell cycle and development. Continued sequencing of pCMS trios will identify additional risk genes and provide greater insights into biological mechanisms of stereotypies across diagnostic boundaries.

Linear discriminant analysis of phenotypic data for classifying autism spectrum disorder by diagnosis and sex.

Autism Spectrum Disorder (ASD) is a developmental condition characterized by social and communication differences. Recent research suggests ASD affects 1-in-44 children in the United States. ASD is diagnosed more commonly in males, though it is unclear whether this diagnostic disparity is a result of a biological predisposition or limitations in diagnostic tools, or both. One hypothesis centers on the 'female protective effect,' which is the theory that females are biologically more resistant to the autism phenotype than males. In this examination, phenotypic data were acquired and combined from four leading research institutions and subjected to multivariate linear discriminant analysis. A linear discriminant model was trained on the training set and then deployed on the test set to predict group membership. Multivariate analyses of variance were performed to confirm the significance of the overall analysis, and individual analyses of variance were performed to confirm the significance of each of the resulting linear discriminant axes. Two discriminant dimensions were identified between the groups: a dimension separating groups by the diagnosis of ASD (LD1: 87% of variance explained); and a dimension reflective of a diagnosis-by-sex interaction (LD2: 11% of variance explained). The strongest discriminant coefficients for the first discriminant axis divided the sample in domains with known differences between ASD and comparison groups, such as social difficulties and restricted repetitive behavior. The discriminant coefficients for the second discriminant axis reveal a more nuanced disparity between boys with ASD and girls with ASD, including executive functioning and high-order behavioral domains as the dominant discriminators. These results indicate that phenotypic differences between males and females with and without ASD are identifiable using parent report measures, which could be utilized to provide additional specificity to the diagnosis of ASD in female patients, potentially leading to more targeted clinical strategies and therapeutic interventions. The study helps to isolate a phenotypic basis for future empirical work on the female protective effect using neuroimaging, EEG, and genomic methodologies.

A neurogenetic analysis of female autism.

Females versus males are less frequently diagnosed with autism spectrum disorder (ASD), and while understanding sex differences is critical to delineating the systems biology of the condition, female ASD is understudied. We integrated functional MRI and genetic data in a sex-balanced sample of ASD and typically developing youth (8-17 years old) to characterize female-specific pathways of ASD risk. Our primary objectives were to: (i) characterize female ASD (n = 45) brain response to human motion, relative to matched typically developing female youth (n = 45); and (ii) evaluate whether genetic data could provide further insight into the potential relevance of these brain functional differences. For our first objective we found that ASD females showed markedly reduced response versus typically developing females, particularly in sensorimotor, striatal, and frontal regions. This difference between ASD and typically developing females does not resemble differences between ASD (n = 47) and typically developing males (n = 47), even though neural response did not significantly differ between female and male ASD. For our second objective, we found that ASD females (n = 61), versus males (n = 66), showed larger median size of rare copy number variants containing gene(s) expressed in early life (10 postconceptual weeks to 2 years) in regions implicated by the typically developing female > female functional MRI contrast. Post hoc analyses suggested this difference was primarily driven by copy number variants containing gene(s) expressed in striatum. This striatal finding was reproducible among n = 2075 probands (291 female) from an independent cohort. Together, our findings suggest that striatal impacts may contribute to pathways of risk in female ASD and advocate caution in drawing conclusions regarding female ASD based on male-predominant cohorts.

De Novo Damaging DNA Coding Mutations Are Associated With Obsessive-Compulsive Disorder and Overlap With Tourette's Disorder and Autism.

BACKGROUND: Obsessive-compulsive disorder (OCD) is a debilitating neuropsychiatric disorder with a genetic risk component, yet identification of high-confidence risk genes has been challenging. In recent years, risk gene discovery in other complex psychiatric disorders has been achieved by studying rare de novo (DN) coding variants. METHODS: We performed whole-exome sequencing in 222 OCD parent-child trios (184 trios after quality control), comparing DN variant frequencies with 777 previously sequenced unaffected trios. We estimated the contribution of DN mutations to OCD risk and the number of genes involved. Finally, we looked for gene enrichment in other datasets and canonical pathways. RESULTS: DN likely gene disrupting and predicted damaging missense variants are enriched in OCD probands (rate ratio, 1.52; p = .0005) and contribute to risk. We identified 2 high-confidence risk genes, each containing 2 DN damaging variants in unrelated probands: CHD8 and SCUBE1. We estimate that 34% of DN damaging variants in OCD contribute to risk and that DN damaging variants in approximately 335 genes contribute to risk in 22% of OCD cases. Furthermore, genes harboring DN damaging variants in OCD are enriched for those reported in neurodevelopmental disorders, particularly Tourette's disorder and autism spectrum disorder. An exploratory network analysis reveals significant functional connectivity and enrichment in canonical pathways, biological processes, and disease networks. CONCLUSIONS: Our findings show a pathway toward systematic gene discovery in OCD via identification of DN damaging variants. Sequencing larger cohorts of OCD parent-child trios will reveal more OCD risk genes and will provide needed insights into underlying disease biology.

PAC1R Genotype to Phenotype Correlations in Autism Spectrum Disorder.

Amygdala dysfunction has been implicated in numerous neurodevelopmental disorders, including autism spectrum disorder (ASD). Previous studies in mice and humans, respectively, have linked Pac1r/PAC1R function to social behavior and PTSD-susceptibility. Based on this connection to social and emotional processing and the central role played by the amygdala in ASD, we examined a putative role for PAC1R in social deficits in ASD and determined the pattern of gene expression in the developing mouse and human amygdala. We reveal that Pac1r/PAC1R is expressed in both mouse and human amygdala from mid-neurogenesis through early postnatal stages, critical time points when altered brain trajectories are hypothesized to unfold in ASD. We further find that parents of autistic children carrying a previously identified PTSD-risk genotype (CC) report greater reciprocal social deficits compared to those carrying the non-risk GC genotype. Additionally, by exploring resting-state functional connectivity differences in a subsample of the larger behavioral sample, we find higher functional connectivity between the amygdala and right middle temporal gyrus in individuals with the CC risk genotype. Thus, using multimodal approaches, our data reveal that the amygdala-expressed PAC1R gene may be linked to severity of ASD social phenotype and possible alterations in brain connectivity, therefore potentially acting as a modifier of amygdala-related phenotypes. Autism Res 2019, 12: 200-211 © 2018 International Society for Autism Research, Wiley Periodicals, Inc. LAY SUMMARY: In this multimodal study across mouse and human, we examined expression patterns of Pac1r/PAC1R, a gene implicated in social behavior, and further explored whether a previously identified human PTSD-linked mutation in PAC1R can predict brain connectivity and social deficits in ASD. We find that PAC1R is highly expressed in the both the mouse and human amygdala. Furthermore, our human data suggest that PAC1R genotype is linked to severity of social deficits and functional amygdala connectivity in ASD.

Neurogenetic analysis of childhood disintegrative disorder.

BACKGROUND: Childhood disintegrative disorder (CDD) is a rare form of autism spectrum disorder (ASD) of unknown etiology. It is characterized by late-onset regression leading to significant intellectual disability (ID) and severe autism. Although there are phenotypic differences between CDD and other forms of ASD, it is unclear if there are neurobiological differences. METHODS: We pursued a multidisciplinary study of CDD (n = 17) and three comparison groups: low-functioning ASD (n = 12), high-functioning ASD (n = 50), and typically developing (n = 26) individuals. We performed whole-exome sequencing (WES), copy number variant (CNV), and gene expression analyses of CDD and, on subsets of each cohort, non-sedated functional magnetic resonance imaging (fMRI) while viewing socioemotional (faces) and non-socioemotional (houses) stimuli and eye tracking while viewing emotional faces. RESULTS: We observed potential differences between CDD and other forms of ASD. WES and CNV analyses identified one or more rare de novo, homozygous, and/or hemizygous (mother-to-son transmission on chrX) variants for most probands that were not shared by unaffected sibling controls. There were no clearly deleterious variants or highly recurrent candidate genes. Candidate genes that were found to be most conserved at variant position and most intolerant of variation, such as TRRAP, ZNF236, and KIAA2018, play a role or may be involved in transcription. Using the human BrainSpan transcriptome dataset, CDD candidate genes were found to be more highly expressed in non-neocortical regions than neocortical regions. This expression profile was similar to that of an independent cohort of ASD probands with regression. The non-neocortical regions overlapped with those identified by fMRI as abnormally hyperactive in response to viewing faces, such as the thalamus, cerebellum, caudate, and hippocampus. Eye-tracking analysis showed that, among individuals with ASD, subjects with CDD focused on eyes the most when shown pictures of faces. CONCLUSIONS: Given that cohort sizes were limited by the rarity of CDD, and the challenges of conducting non-sedated fMRI and eye tracking in subjects with ASD and significant ID, this is an exploratory study designed to investigate the neurobiological features of CDD. In addition to reporting the first multimodal analysis of CDD, a combination of fMRI and eye-tracking analyses are being presented for the first time for low-functioning individuals with ASD. Our results suggest differences between CDD and other forms of ASD on the neurobiological as well as clinical level.

Automated quantitative analysis of tissue microarray of 443 patients with colorectal adenocarcinoma: low expression of Bcl-2 predicts poor survival.

BACKGROUND: Bcl-2 has been implicated in the development and progression of a number of cancers including colorectal cancer. Reports of Bcl-2 expression in colorectal cancer and patient outcomes have been inconsistent due to small cohorts and semi-quantitative grading methods. STUDY DESIGN: We used a high throughput tissue microarray system (automated quantitative analysis [AQUA]), analyzing colorectal adenocarcinoma samples from 443 patients resected during the period of 1967 to 1986. This system uses fully quantitative, automated fluorescent microscopy to accurately assess Bcl-2 expression in colorectal cancer samples. Clinicopathologic variables were collected prospectively and were assessed using log-rank tests and Cox proportional hazards models. RESULTS: At a median follow-up of 54 months, the 5- and 10-year disease-specific survivals for all patients were 59.2% and 52.1%, respectively. Loss of Bcl-2 expression was seen in 70.4% of tumors and was associated with a decreased 5-year disease-specific survival (55.8% vs 75.6%, p = 0.001 and relative risk [RR] 1.8) and decreased 5-year overall survival (45.8% vs 56.5%, p = 0.046 and RR 1.2). On univariate analysis, T stage, N stage, and loss of Bcl-2 expression predicted poor disease-specific survival. On multivariate analysis, Bcl-2 expression was an independent prognostic factor for disease-specific survival (p = 0.034). CONCLUSIONS: Our results indicate that loss of Bcl-2 expression in colorectal cancer is associated with decreased disease-specific and overall survival. This finding may help identify a subset of patients with a more aggressive phenotype and guide adjuvant chemotherapy choices.

A germline mutation in the BRCA1 3'UTR predicts Stage IV breast cancer.

BACKGROUND: A germline, variant in the BRCA1 3'UTR (rs8176318) was previously shown to predict breast and ovarian cancer risk in women from high-risk families, as well as increased risk of triple negative breast cancer. Here, we tested the hypothesis that this variant predicts tumor biology, like other 3'UTR mutations in cancer. METHODS: The impact of the BRCA1-3'UTR-variant on BRCA1 gene expression, and altered response to external stimuli was tested in vitro using a luciferase reporter assay. Gene expression was further tested in vivo by immunoflourescence staining on breast tumor tissue, comparing triple negative patient samples with the variant (TG or TT) or non-variant (GG) BRCA1 3'UTR. To determine the significance of the variant on clinically relevant endpoints, a comprehensive collection of West-Irish breast cancer patients were tested for the variant. Finally, an association of the variant with breast screening clinical phenotypes was evaluated using a cohort of women from the High Risk Breast Program at the University of Vermont. RESULTS: Luciferase reporters with the BRCA1-3'UTR-variant (T allele) displayed significantly lower gene expression, as well as altered response to external hormonal stimuli, compared to the non-variant 3'UTR (G allele) in breast cancer cell lines. This was confirmed clinically by the finding of reduced BRCA1 gene expression in triple negative samples from patients carrying the homozygous TT variant, compared to non-variant patients. The BRCA1-3'UTR-variant (TG or TT) also associated with a modest increased risk for developing breast cancer in the West-Irish cohort (OR=1.4, 95% CI 1.1-1.8, p=0.033). More importantly, patients with the BRCA1-3'UTR-variant had a 4-fold increased risk of presenting with Stage IV disease (p=0.018, OR=3.37, 95% CI 1.3-11.0). Supporting that this finding is due to tumor biology, and not difficulty screening, obese women with the BRCA1-3'UTR-variant had significantly less dense breasts (p=0.0398) in the Vermont cohort. CONCLUSION: A variant in the 3'UTR of BRCA1 is functional, leading to decreased BRCA1 expression, modest increased breast cancer risk, and most importantly, presentation with stage IV breast cancer, likely due to aggressive tumor biology.

Tissue microarray analysis of 560 patients with colorectal adenocarcinoma: high expression of HuR predicts poor survival.

The purpose of this study is to characterize the expression of HuR in colorectal carcinoma and determine its correlation with clinical outcome. Differential expression of HuR has been suggested to be of prognostic significance in carcinomas of the ovaries, stomach, and breast. HuR regulates the expression of a variety of proteins critical to carcinogenesis via the pathways of cell-cycle progress, invasion, and metastasis. Increasing evidence suggests that angiogenic pathways are involved. A tissue microarray consisting of tumors from 560 patients with colorectal adenocarcinoma was analyzed for HuR protein expression using a quantitative, automated immunofluorescent microscopy system (AQUA). Clinical data corresponding to each examined specimen collected through an institutional review board (IRB)-approved protocol were analyzed using chi-squared test, Cox regression, and Kaplan-Meier analysis. Median follow-up was 54 months. Along with tumor stage and overall tumor-node-metastasis (TNM) stage, HuR expression was found to be an independent predictor of survival. In patients in the highest quartile of total HuR expression, survival was 22.8 months less than those in the lower quartiles (40.6 versus 63.4 months, p = 0.04). Furthermore, HuR levels correlate positively with expression of vascular endothelial growth factor (VEGF) and CD31, a marker for vascular endothelium. We conclude that expression of high levels of HuR correlates with features of advanced disease and portends poorer survival in patients with colorectal adenocarcinoma. These results further suggest that HuR exerts its tumorigenic effects through VEGF-mediated angiogenesis and may be a novel therapeutic target in colorectal cancer.

High levels of vascular endothelial growth factor and its receptors (VEGFR-1, VEGFR-2, neuropilin-1) are associated with worse outcome in breast cancer.

Vascular endothelial growth factor has been shown to be up-regulated in breast cancers. Vascular endothelial growth factor receptors, VEGFR-1 and VEGFR-2, are the principal mediators of its effects. Together with VEGFR-1 and VEGFR-2, neuropilin-1 may act as a coreceptor for vascular endothelial growth factor. Although vascular endothelial growth factor exerts important effects on endothelial cells, VEGFRs are likely present on tumor cells as well. We used AQUA to analyze tumor-specific expression of vascular endothelial growth factor, VEGFR-1, VEGFR-2, and neuropilin-1 on a large cohort of breast cancer tissue microarray. Two-fold redundant arrays were constructed from 642 cases of primary breast adenocarcinomas. Automated image analysis with AQUA (Automated Quantitative Analysis) was then performed to determine a quantitative expression score. Scores from redundant arrays were normalized and averaged. Kaplan-Meier survival analysis showed that high levels of vascular endothelial growth factor, VEGFR-1, VEGFR-2, and neuropilin-1 were all significantly associated with survival (Miller Siegmeund corrected P = .0020, .0160, and .0320, respectively). In addition, vascular endothelial growth factor and neuropilin-1 retained a significant association with survival independent of other standard prognostic factors. Vascular endothelial growth factor, VEGFR-1 and -2, and neuropilin-1 are expressed to varying degrees in primary breast cancers and have prognostic significance. Further study of the functional significance of this finding is warranted as well as the prognostic value of these biomarkers in other tumor microenvironment-specific compartments (eg, vessels).

Biomarker validation: in situ analysis of protein expression using semiquantitative immunohistochemistry-based techniques.

Biomarker-driven cancer research is common in the current literature. Much of this research is a result of the increase in genomic and proteomic high-throughput technologies, which have increased our knowledge and also produced an abundance of data with unclear clinical significance. Immunohistochemistry-based assessment of protein expression is a natural validation method of expression-profiling data that is easily performed on tissue samples collected prospectively or from archived samples. Coupled with tissue microarray technology and the increasing number of available automated, quantitative systems to read these arrays, we now have an efficient method of validating biomarkers for prognostic and predictive capabilities and for the identification of drug development targets.