Monoallelic de novo AJAP1 loss-of-function variants disrupt trans-synaptic control of neurotransmitter release.

Adherens junction-associated protein 1 (AJAP1) has been implicated in brain diseases; however, a pathogenic mechanism has not been identified. AJAP1 is widely expressed in neurons and binds to γ-aminobutyric acid type B receptors (GBRs), which inhibit neurotransmitter release at most synapses in the brain. Here, we show that AJAP1 is selectively expressed in dendrites and trans-synaptically recruits GBRs to presynaptic sites of neurons expressing AJAP1. We have identified several monoallelic AJAP1 variants in individuals with epilepsy and/or neurodevelopmental disorders. Specifically, we show that the variant p.(W183C) lacks binding to GBRs, resulting in the inability to recruit them. Ultrastructural analysis revealed significantly decreased presynaptic GBR levels in Ajap1-/- and Ajap1W183C/+ mice. Consequently, these mice exhibited reduced GBR-mediated presynaptic inhibition at excitatory and inhibitory synapses, along with impaired synaptic plasticity. Our study reveals that AJAP1 enables the postsynaptic neuron to regulate the level of presynaptic GBR-mediated inhibition, supporting the clinical relevance of loss-of-function AJAP1 variants.

Noise induces intercellular Ca 2+ signaling waves and the unfolded protein response in the hearing cochlea.

Exposure to loud noise is a common cause of acquired hearing loss. Disruption of subcellular calcium homeostasis and downstream stress pathways in the endoplasmic reticulum and mitochondria, including the unfolded protein response, have been implicated in the pathophysiology of noise-induced hearing loss. However, studies on the association between calcium homeostasis and stress pathways has been limited due to limited ability to measure calcium dynamics in mature-hearing, noise-exposed mice. We used a genetically encoded calcium indicator mouse model in which GcAMP is expressed specifically in hair cells or supporting cells under control of Myo15Cre or Sox2Cre, respectively. We performed live calcium imaging and UPR gene expression analysis in 8-week-old mice exposed to levels of noise that cause cochlear synaptopathy (98 db SPL) or permanent hearing loss (106 dB SPL). UPR activation occurred immediately after noise exposure and was noise dose-dependent, with the pro-apoptotic pathway upregulated only after 106 dB noise exposure. Spontaneous calcium transients in hair cells and intercellular calcium waves in supporting cells, which are present in neonatal cochleae, were quiescent in mature-hearing cochleae, but re-activated upon noise exposure. 106 dB noise exposure was associated with more persistent and expansive ICS wave activity. These findings demonstrate a strong and dose-dependent association between noise exposure, UPR activation, and changes in calcium homeostasis in hair cells and supporting cells, suggesting that targeting these pathways may be effective to develop treatments for noise-induced hearing loss.

Discovering the gene-brain-behavior link in autism via generative machine learning.

Autism is traditionally diagnosed behaviorally but has a strong genetic basis. A genetics-first approach could transform understanding and treatment of autism. However, isolating the gene-brain-behavior relationship from confounding sources of variability is a challenge. We demonstrate a novel technique, 3D transport-based morphometry (TBM), to extract the structural brain changes linked to genetic copy number variation (CNV) at the 16p11.2 region. We identified two distinct endophenotypes. In data from the Simons Variation in Individuals Project, detection of these endophenotypes enabled 89 to 95% test accuracy in predicting 16p11.2 CNV from brain images alone. Then, TBM enabled direct visualization of the endophenotypes driving accurate prediction, revealing dose-dependent brain changes among deletion and duplication carriers. These endophenotypes are sensitive to articulation disorders and explain a portion of the intelligence quotient variability. Genetic stratification combined with TBM could reveal new brain endophenotypes in many neurodevelopmental disorders, accelerating precision medicine, and understanding of human neurodiversity.

Exome copy number variant detection, analysis, and classification in a large cohort of families with undiagnosed rare genetic disease.

Copy number variants (CNVs) are significant contributors to the pathogenicity of rare genetic diseases and, with new innovative methods, can now reliably be identified from exome sequencing. Challenges still remain in accurate classification of CNV pathogenicity. CNV calling using GATK-gCNV was performed on exomes from a cohort of 6,633 families (15,759 individuals) with heterogeneous phenotypes and variable prior genetic testing collected at the Broad Institute Center for Mendelian Genomics of the Genomics Research to Elucidate the Genetics of Rare Diseases consortium and analyzed using the seqr platform. The addition of CNV detection to exome analysis identified causal CNVs for 171 families (2.6%). The estimated sizes of CNVs ranged from 293 bp to 80 Mb. The causal CNVs consisted of 140 deletions, 15 duplications, 3 suspected complex structural variants (SVs), 3 insertions, and 10 complex SVs, the latter two groups being identified by orthogonal confirmation methods. To classify CNV variant pathogenicity, we used the 2020 American College of Medical Genetics and Genomics/ClinGen CNV interpretation standards and developed additional criteria to evaluate allelic and functional data as well as variants on the X chromosome to further advance the framework. We interpreted 151 CNVs as likely pathogenic/pathogenic and 20 CNVs as high-interest variants of uncertain significance. Calling CNVs from existing exome data increases the diagnostic yield for individuals undiagnosed after standard testing approaches, providing a higher-resolution alternative to arrays at a fraction of the cost of genome sequencing. Our improvements to the classification approach advances the systematic framework to assess the pathogenicity of CNVs.

TMTC4 is a hair cell-specific human deafness gene.

Transmembrane and tetratricopeptide repeat 4 (Tmtc4) is a deafness gene in mice. Tmtc4-KO mice have rapidly progressive postnatal hearing loss due to overactivation of the unfolded protein response (UPR); however, the cellular basis and human relevance of Tmtc4-associated hearing loss in the cochlea was not heretofore appreciated. We created a hair cell-specific conditional KO mouse that phenocopies the constitutive KO with postnatal onset deafness, demonstrating that Tmtc4 is a hair cell-specific deafness gene. Furthermore, we identified a human family in which Tmtc4 variants segregate with adult-onset progressive hearing loss. Lymphoblastoid cells derived from multiple affected and unaffected family members, as well as human embryonic kidney cells engineered to harbor each of the variants, demonstrated that the human Tmtc4 variants confer hypersensitivity of the UPR toward apoptosis. These findings provide evidence that TMTC4 is a deafness gene in humans and further implicate the UPR in progressive hearing loss.

Modulating the Unfolded Protein Response with ISRIB Mitigates Cisplatin Ototoxicity.

Cisplatin is a common chemotherapy drug with a nearly universal side effect of ototoxicity. The cellular mechanisms underlying cisplatin ototoxicity are poorly understood. Efforts in drug development to prevent or reverse cisplatin ototoxicity have largely focused on pathways of oxidative stress and apoptosis. An effective treatment for cisplatin ototoxicity, sodium thiosulfate, is associated with reduced survival in disseminated hepatoblastoma, highlighting the need for more specific drugs. The unfolded protein response (UPR) and endoplasmic reticulum (ER) stress pathways have been shown to be involved in the pathogenesis of noise-induced hearing loss and cochlear synaptopathy in vivo , and these pathways have been implicated broadly in cisplatin cytotoxicity. This study sought to determine whether the UPR can be targeted to prevent cisplatin ototoxicity. Neonatal cochlear cultures and HEK cells were exposed to cisplatin and UPR-modulating drugs, and UPR marker gene expression and cell death measured. Treatment with ISRIB, a drug that activates eif2B and downregulates the pro-apoptotic PERK/CHOP pathway of the UPR, was tested in an in vivo mouse model of cisplatin ototoxicity and well as a head and neck squamous cell carcinoma (HNSCC) cell-based assay of cisplatin cytotoxicity. Cisplatin exhibited a biphasic, non-linear dose-response of cell death and apoptosis that correlated with different patterns of UPR marker gene expression in HEK cells and cochlear cultures. ISRIB treatment protected against cisplatin-induced hearing loss and hair-cell death, but did not impact cisplatin's cytotoxic effects on HNSCC cell viability. These findings demonstrate that targeting the pro-apoptotic PERK/CHOP pathway with ISRIB can mitigate cisplatin ototoxicity without reducing anti-cancer cell effects, suggesting that this may be a viable strategy for drug development.

Exome copy number variant detection, analysis and classification in a large cohort of families with undiagnosed rare genetic disease.

Copy number variants (CNVs) are significant contributors to the pathogenicity of rare genetic diseases and with new innovative methods can now reliably be identified from exome sequencing. Challenges still remain in accurate classification of CNV pathogenicity. CNV calling using GATK-gCNV was performed on exomes from a cohort of 6,633 families (15,759 individuals) with heterogeneous phenotypes and variable prior genetic testing collected at the Broad Institute Center for Mendelian Genomics of the GREGoR consortium. Each family's CNV data was analyzed using the seqr platform and candidate CNVs classified using the 2020 ACMG/ClinGen CNV interpretation standards. We developed additional evidence criteria to address situations not covered by the current standards. The addition of CNV calling to exome analysis identified causal CNVs for 173 families (2.6%). The estimated sizes of CNVs ranged from 293 bp to 80 Mb with estimates that 44% would not have been detected by standard chromosomal microarrays. The causal CNVs consisted of 141 deletions, 15 duplications, 4 suspected complex structural variants (SVs), 3 insertions and 10 complex SVs, the latter two groups being identified by orthogonal validation methods. We interpreted 153 CNVs as likely pathogenic/pathogenic and 20 CNVs as high interest variants of uncertain significance. Calling CNVs from existing exome data increases the diagnostic yield for individuals undiagnosed after standard testing approaches, providing a higher resolution alternative to arrays at a fraction of the cost of genome sequencing. Our improvements to the classification approach advances the systematic framework to assess the pathogenicity of CNVs.

The variant landscape and function of DDX3X in cancer and neurodevelopmental disorders.

RNA molecules rely on proteins across their life cycle. DDX3X encodes an X-linked DEAD-box RNA helicase with a Y-linked paralog, DDX3Y. DDX3X is central to the RNA life cycle and is implicated in many conditions, including cancer and the neurodevelopmental disorder DDX3X syndrome. DDX3X-linked conditions often exhibit sex differences, possibly due to differences between expression or function of the X- and Y-linked paralogs DDX3X and DDX3Y. DDX3X-related diseases have different mutational landscapes, indicating different roles of DDX3X. Understanding the role of DDX3X in normal and disease states will inform the understanding of DDX3X in disease. We review the function of DDX3X and DDX3Y, discuss how mutation type and sex bias contribute to human diseases involving DDX3X, and review possible DDX3X-targeting treatments.

De novo CLCN3 variants affecting Gly327 cause severe neurodevelopmental syndrome with brain structural abnormalities.

A recent study revealed that monoallelic missense or biallelic loss-of-function variants in the chloride voltage-gated channel 3 (CLCN3) cause neurodevelopmental disorders resulting in brain abnormalities. Functional studies suggested that some missense variants had varying gain-of-function effects on channel activity. Meanwhile, two patients with homozygous frameshift variants showed severe neuropsychiatric disorders and a range of brain structural abnormalities. Here we describe two patients with de novo CLCN3 variants affecting the same amino acid, Gly327 (p.(Gly327Ser) and p.(Gly327Asp)). They showed severe neurological phenotypes including global developmental delay, intellectual disability, hypotonia, failure to thrive, and various brain abnormalities. They also presented with characteristic brain and ophthalmological abnormalities, hippocampal and retinal degradation, which were observed in patients harboring homozygous loss-of-function variants. These findings were also observed in CLCN3-deficient mice, indicating that the monoallelic missense variant may also have a dominant negative effect. This study will expand the phenotypic spectrum of CLCN3-related disorders.

Clinical Characteristics of Seizures and Epilepsy in Individuals With Recurrent Deletions and Duplications in the 16p11.2 Region.

Background and Objectives: Deletions and duplications at 16p11.2 (BP4 to BP5; 29.5-30.1 Mb) have been associated with several neurodevelopmental and neuropsychiatric disorders including autism spectrum disorder, intellectual disability (ID), and schizophrenia. Seizures have also been reported in individuals with these particular copy number variants, but the epilepsy phenotypes have not been well-delineated. We aimed to systematically characterize the seizure types, epilepsy syndromes, and epilepsy severity in a large cohort of individuals with these 16p11.2 deletions and duplications. Methods: The cohort of ascertained participants with the recurrent 16p11.2 copy number variant was assembled through the multicenter Simons Variation in Individuals Project. Detailed data on individuals identified as having a history of seizures were obtained using a semistructured phone interview and review of medical records, EEG, and MRI studies obtained clinically or as part of the Simons Variation in Individuals Project. Results: Among 129 individuals with the 16p11.2 deletion, 31 (24%) had at least one seizure, including 23 (18%) who met criteria for epilepsy; 42% of them fit the phenotype of classic or atypical Self-limited (Familial) Infantile Epilepsy (Se(F)IE). Among 106 individuals with 16p11.2 duplications, 16 (15%) had at least one seizure, including 11 (10%) who met criteria for epilepsy. The seizure types and epilepsy syndromes were heterogeneous in this group. Most of the individuals in both the deletion and duplication groups had well-controlled seizures with subsequent remission. Pharmacoresistant epilepsy was uncommon. Seizures responded favorably to phenobarbital, carbamazepine, and oxcarbazepine in the deletion group, specifically in the Se(F)IE, and to various antiseizure medications in the duplication group. Discussion: These findings delineate the spectrum of seizures and epilepsies in the recurrent 16p11.2 deletions and duplications and provide potential diagnostic, therapeutic, and prognostic information.

Systematic analysis and prediction of genes associated with monogenic disorders on human chromosome X.

Disease gene discovery on chromosome (chr) X is challenging owing to its unique modes of inheritance. We undertook a systematic analysis of human chrX genes. We observe a higher proportion of disorder-associated genes and an enrichment of genes involved in cognition, language, and seizures on chrX compared to autosomes. We analyze gene constraints, exon and promoter conservation, expression, and paralogues, and report 127 genes sharing one or more attributes with known chrX disorder genes. Using machine learning classifiers trained to distinguish disease-associated from dispensable genes, we classify 247 genes, including 115 of the 127, as having high probability of being disease-associated. We provide evidence of an excess of variants in predicted genes in existing databases. Finally, we report damaging variants in CDK16 and TRPC5 in patients with intellectual disability or autism spectrum disorders. This study predicts large-scale gene-disease associations that could be used for prioritization of X-linked pathogenic variants.

Loss of non-motor kinesin KIF26A causes congenital brain malformations via dysregulated neuronal migration and axonal growth as well as apoptosis.

Kinesins are canonical molecular motors but can also function as modulators of intracellular signaling. KIF26A, an unconventional kinesin that lacks motor activity, inhibits growth-factor-receptor-bound protein 2 (GRB2)- and focal adhesion kinase (FAK)-dependent signal transduction, but its functions in the brain have not been characterized. We report a patient cohort with biallelic loss-of-function variants in KIF26A, exhibiting a spectrum of congenital brain malformations. In the developing brain, KIF26A is preferentially expressed during early- and mid-gestation in excitatory neurons. Combining mice and human iPSC-derived organoid models, we discovered that loss of KIF26A causes excitatory neuron-specific defects in radial migration, localization, dendritic and axonal growth, and apoptosis, offering a convincing explanation of the disease etiology in patients. Single-cell RNA sequencing in KIF26A knockout organoids revealed transcriptional changes in MAPK, MYC, and E2F pathways. Our findings illustrate the pathogenesis of KIF26A loss-of-function variants and identify the surprising versatility of this non-motor kinesin.

Heterozygous variants in MYH10 associated with neurodevelopmental disorders and congenital anomalies with evidence for primary cilia-dependent defects in Hedgehog signaling.

PURPOSE: Nonmuscle myosin II complexes are master regulators of actin dynamics that play essential roles during embryogenesis with vertebrates possessing 3 nonmuscle myosin II heavy chain genes, MYH9, MYH10, and MYH14. As opposed to MYH9 and MYH14, no recognizable disorder has been associated with MYH10. We sought to define the clinical characteristics and molecular mechanism of a novel autosomal dominant disorder related to MYH10. METHODS: An international collaboration identified the patient cohort. CAS9-mediated knockout cell models were used to explore the mechanism of disease pathogenesis. RESULTS: We identified a cohort of 16 individuals with heterozygous MYH10 variants presenting with a broad spectrum of neurodevelopmental disorders and variable congenital anomalies that affect most organ systems and were recapitulated in animal models of altered MYH10 activity. Variants were typically de novo missense changes with clustering observed in the motor domain. MYH10 knockout cells showed defects in primary ciliogenesis and reduced ciliary length with impaired Hedgehog signaling. MYH10 variant overexpression produced a dominant-negative effect on ciliary length. CONCLUSION: These data presented a novel genetic cause of isolated and syndromic neurodevelopmental disorders related to heterozygous variants in the MYH10 gene with implications for disrupted primary cilia length control and altered Hedgehog signaling in disease pathogenesis.

Newborn screening for neurodevelopmental diseases: Are we there yet?

In the US, newborn screening (NBS) is a unique health program that supports health equity and screens virtually every baby after birth, and has brought timely treatments to babies since the 1960's. With the decreasing cost of sequencing and the improving methods to interpret genetic data, there is an opportunity to add DNA sequencing as a screening method to facilitate the identification of babies with treatable conditions that cannot be identified in any other scalable way, including highly penetrant genetic neurodevelopmental disorders (NDD). However, the lack of effective dietary or drug-based treatments has made it nearly impossible to consider NDDs in the current NBS framework, yet it is anticipated that any treatment will be maximally effective if started early. Hence there is a critical need for large scale pilot studies to assess if and how NDDs can be effectively screened at birth, if parents desire that information, and what impact early diagnosis may have. Here we attempt to provide an overview of the recent advances in NDD treatments, explore the possible framework of setting up a pilot study to genetically screen for NDDs, highlight key technical, practical, and ethical considerations and challenges, and examine the policy and health system implications.

ANKLE2-related microcephaly: A variable microcephaly syndrome resembling Zika infection.

OBJECTIVE: This study delineates the clinical and molecular spectrum of ANKLE2-related microcephaly (MIC), as well as highlights shared pathological mechanisms between ANKLE2 and the Zika virus. METHODS: We identified 12 individuals with MIC and variants in ANKLE2 with a broad range of features. Probands underwent thorough phenotypic evaluations, developmental assessments, and anthropometric measurements. Brain imaging studies were systematically reviewed for developmental abnormalities. We functionally interrogated a subset of identified ANKLE2 variants in Drosophila melanogaster. RESULTS: All individuals had MIC (z-score ≤ -3), including nine with congenital MIC. We identified a broad range of brain abnormalities including simplified cortical gyral pattern, full or partial callosal agenesis, increased extra-axial spaces, hypomyelination, cerebellar vermis hypoplasia, and enlarged cisterna magna. All probands had developmental delays in at least one domain, with speech and language delays being the most common. Six probands had skin findings characteristic of ANKLE2 including hyper- and hypopigmented macules. Only one individual had scalp rugae. Functional characterization in Drosophila recapitulated the human MIC phenotype. Of the four variants tested, p.Val229Gly, p.Arg236*, and p.Arg536Cys acted as partial-loss-of-function variants, whereas the c.1421-1G>C splicing variant demonstrated a strong loss-of-function effect. INTERPRETATION: Deleterious variants in the ANKLE2 gene cause a unique MIC syndrome characterized by congenital or postnatal MIC, a broad range of structural brain abnormalities, and skin pigmentary changes. Thorough functional characterization has identified shared pathogenic mechanisms between ANKLE2-related MIC and congenital Zika virus infection. This study further highlights the importance of a thorough diagnostic evaluation including molecular diagnostic testing in individuals with MIC.

De Novo GLI3 Pathogenic Variants May Cause Hypotonia and a Range of Brain Malformations Without Skeletal Abnormalities.

BACKGROUND: GLI3 encodes a zinc finger transcription factor that plays a role in the sonic hedgehog pathway. Germline pathogenic GLI3 variants are associated with Greig cephalopolysyndactyly and Pallister-Hall syndromes, two syndromes involving brain malformation and polydactyly. METHODS: We identified patients with pathogenic GLI3 variants and brain malformations in the absence of polydactyly or other skeletal malformation. RESULTS: Two patients were identified. Patient #1 is a 4-year-old boy with hypotonia and global developmental delay. Brain MRI showed a focal cortical dysplasia, but he had no history of seizures. Genetic testing identified a de novo likely pathogenic GLI3 variant: c.4453A>T, p.Asn1485Tyr. Patient #2 is a 4-year-old boy with hypotonia, macrocephaly, and global developmental delay. His brain MRI showed partial agenesis of the corpus callosum, dilatation of the right lateral ventricle, and absent hippocampal commissure. Genetic testing identified a de novo pathogenic GLI3 variant: c.4236_4237del, p.Gln1414AspfsTer21. Neither patient had polydactyly or any apparent skeletal abnormality. CONCLUSIONS: These patients widen the spectrum of clinical features that may be associated with GLI3 pathogenic variants to include hypotonia, focal cortical dysplasia, and other brain malformations, in the absence of apparent skeletal malformation. Further study is needed to determine if GLI3 pathogenic variants are a more common cause of focal cortical dysplasia or corpus callosum agenesis than presently recognized.

Whole-exome sequencing with targeted analysis and epilepsy after acute symptomatic neonatal seizures.

BACKGROUND: The contribution of pathogenic gene variants with development of epilepsy after acute symptomatic neonatal seizures is not known. METHODS: Case-control study of 20 trios in children with a history of acute symptomatic neonatal seizures: 10 with and 10 without post-neonatal epilepsy. We performed whole-exome sequencing (WES) and identified pathogenic de novo, transmitted, and non-transmitted variants from established and candidate epilepsy association genes and correlated prevalence of these variants with epilepsy outcomes. We performed a sensitivity analysis with genes associated with coronary artery disease (CAD). We analyzed variants throughout the exome to evaluate for differential enrichment of functional properties using exploratory KEGG searches. RESULTS: Querying 200 established and candidate epilepsy genes, pathogenic variants were identified in 5 children with post-neonatal epilepsy yet in only 1 child without subsequent epilepsy. There was no difference in the number of trios with non-transmitted pathogenic variants in epilepsy or CAD genes. An exploratory KEGG analysis demonstrated a relative enrichment in cell death pathways in children without subsequent epilepsy. CONCLUSIONS: In this pilot study, children with epilepsy after acute symptomatic neonatal seizures had a higher prevalence of coding variants with a targeted epilepsy gene sequencing analysis compared to those patients without subsequent epilepsy. IMPACT: We performed whole-exome sequencing (WES) in 20 trios, including 10 children with epilepsy and 10 without epilepsy, both after acute symptomatic neonatal seizures. Children with post-neonatal epilepsy had a higher burden of pathogenic variants in epilepsy-associated genes compared to those without post-neonatal epilepsy. Future studies evaluating this association may lead to a better understanding of the risk of epilepsy after acute symptomatic neonatal seizures and elucidate molecular pathways that are dysregulated after brain injury and implicated in epileptogenesis.

Major brain malformations: corpus callosum dysgenesis, agenesis of septum pellucidum and polymicrogyria in patients with BCORL1-related disorders.

OBJECTIVE: BCORL1, a transcriptional co-repressor, has a role in cortical migration, neuronal differentiation, maturation, and cerebellar development. We describe BCORL1 as a new genetic cause for major brain malformations. METHODS AND RESULTS: We report three patients from two unrelated families with neonatal onset intractable epilepsy and profound global developmental delay. Brain MRI of two siblings from the first family depicted hypoplastic corpus callosum and septal agenesis (ASP) in the older brother and unilateral perisylvian polymicrogyria (PMG) in the younger one. MRI of the patient from the second family demonstrated complete agenesis of corpus callosum (CC). Whole Exome Sequencing revealed a novel hemizygous variant in NM_021946.5 (BCORL1):c.796C>T (p.Pro266Ser) in the two siblings from the first family and the NM_021946.5 (BCORL1): c.3376G>A; p.Asp1126Asn variant in the patient from the second family, both variants inherited from healthy mothers. We reviewed the patients' charts and MRIs and compared the phenotype to the other published BCORL1-related cases. Brain malformations have not been previously described in association with the BCORL1 phenotype. We discuss the potential influence of BCORL1 on brain development. CONCLUSIONS: We suggest that BCORL1 variants present with a spectrum of neurodevelopmental disorders and can lead to major brain malformations originating at different stages of fetal development. We suggest adding BCORL1 to the genetic causes of PMG, ASP, and CC dysgenesis.

Autism-associated biomarkers: test-retest reliability and relationship to quantitative social trait variation in rhesus monkeys.

BACKGROUND: Rhesus monkeys (Macaca mulatta) exhibit pronounced individual differences in social traits as measured by the macaque Social Responsiveness Scale-Revised. The macaque Social Responsiveness Scale was previously adapted from the Social Responsiveness Scale, an instrument designed to assess social and autistic trait variation in humans. To better understand potential biological underpinnings of this behavioral variation, we evaluated the trait-like consistency of several biological measures previously implicated in autism (e.g., arginine vasopressin, oxytocin, and their receptors, as well as ERK1/2, PTEN, and AKT(1-3) from the RAS-MAPK and PI3K-AKT pathways). We also tested which biological measures predicted macaque Social Responsiveness Scale-Revised scores. METHODS: Cerebrospinal fluid and blood samples were collected from N = 76 male monkeys, which, as a sample, showed a continuous distribution on the macaque Social Responsiveness Scale-Revised. In a subset of these subjects (n = 43), samples were collected thrice over a 10-month period. The following statistical tests were used: "Case 2A" intra-class correlation coefficients of consistency, principal component analysis, and general linear modeling. RESULTS: All biological measures (except AKT) showed significant test-retest reliability within individuals across time points. We next performed principal component analysis on data from monkeys with complete biological measurement sets at the first time point (n = 57), to explore potential correlations between the reliable biological measures and their relationship to macaque Social Responsiveness Scale-Revised score; a three-component solution was found. Follow-up analyses revealed that cerebrospinal fluid arginine vasopressin concentration, but no other biological measure, robustly predicted individual differences in macaque Social Responsiveness Scale-Revised scores, such that monkeys with the lowest cerebrospinal fluid arginine vasopressin concentration exhibited the greatest social impairment. Finally, we confirmed that this result held in the larger study sample (in which cerebrospinal fluid arginine vasopressin values were available from n = 75 of the subjects). CONCLUSIONS: These findings indicate that cerebrospinal fluid arginine vasopressin concentration is a stable trait-like measure and that it is linked to quantitative social trait variation in male rhesus monkeys.