Structurally divergent and recurrently mutated regions of primate genomes.

We sequenced and assembled using multiple long-read sequencing technologies the genomes of chimpanzee, bonobo, gorilla, orangutan, gibbon, macaque, owl monkey, and marmoset. We identified 1,338,997 lineage-specific fixed structural variants (SVs) disrupting 1,561 protein-coding genes and 136,932 regulatory elements, including the most complete set of human-specific fixed differences. We estimate that 819.47 Mbp or ∼27% of the genome has been affected by SVs across primate evolution. We identify 1,607 structurally divergent regions wherein recurrent structural variation contributes to creating SV hotspots where genes are recurrently lost (e.g., CARD, C4, and OLAH gene families) and additional lineage-specific genes are generated (e.g., CKAP2, VPS36, ACBD7, and NEK5 paralogs), becoming targets of rapid chromosomal diversification and positive selection (e.g., RGPD gene family). High-fidelity long-read sequencing has made these dynamic regions of the genome accessible for sequence-level analyses within and between primate species.

Efficient embryonic homozygous gene conversion via RAD51-enhanced interhomolog repair.

Searching for factors to improve knockin efficiency for therapeutic applications, biotechnology, and generation of non-human primate models of disease, we found that the strand exchange protein RAD51 can significantly increase Cas9-mediated homozygous knockin in mouse embryos through an interhomolog repair (IHR) mechanism. IHR is a hallmark of meiosis but only occurs at low frequencies in somatic cells, and its occurrence in zygotes is controversial. Using multiple approaches, we provide evidence for an endogenous IHR mechanism in the early embryo that can be enhanced by RAD51. This process can be harnessed to generate homozygotes from wild-type zygotes using exogenous donors and to convert heterozygous alleles into homozygous alleles without exogenous templates. Furthermore, we identify additional IHR-promoting factors and describe features of IHR events. Together, our findings show conclusive evidence for IHR in mouse embryos and describe an efficient method for enhanced gene conversion.

Targeting Peripheral Somatosensory Neurons to Improve Tactile-Related Phenotypes in ASD Models.

Somatosensory over-reactivity is common among patients with autism spectrum disorders (ASDs) and is hypothesized to contribute to core ASD behaviors. However, effective treatments for sensory over-reactivity and ASDs are lacking. We found distinct somatosensory neuron pathophysiological mechanisms underlie tactile abnormalities in different ASD mouse models and contribute to some ASD-related behaviors. Developmental loss of ASD-associated genes Shank3 or Mecp2 in peripheral mechanosensory neurons leads to region-specific brain abnormalities, revealing links between developmental somatosensory over-reactivity and the genesis of aberrant behaviors. Moreover, acute treatment with a peripherally restricted GABAA receptor agonist that acts directly on mechanosensory neurons reduced tactile over-reactivity in six distinct ASD models. Chronic treatment of Mecp2 and Shank3 mutant mice improved body condition, some brain abnormalities, anxiety-like behaviors, and some social impairments but not memory impairments, motor deficits, or overgrooming. Our findings reveal a potential therapeutic strategy targeting peripheral mechanosensory neurons to treat tactile over-reactivity and select ASD-related behaviors.

CRISPR-Cas9 knockin mice for genome editing and cancer modeling.

CRISPR-Cas9 is a versatile genome editing technology for studying the functions of genetic elements. To broadly enable the application of Cas9 in vivo, we established a Cre-dependent Cas9 knockin mouse. We demonstrated in vivo as well as ex vivo genome editing using adeno-associated virus (AAV)-, lentivirus-, or particle-mediated delivery of guide RNA in neurons, immune cells, and endothelial cells. Using these mice, we simultaneously modeled the dynamics of KRAS, p53, and LKB1, the top three significantly mutated genes in lung adenocarcinoma. Delivery of a single AAV vector in the lung generated loss-of-function mutations in p53 and Lkb1, as well as homology-directed repair-mediated Kras(G12D) mutations, leading to macroscopic tumors of adenocarcinoma pathology. Together, these results suggest that Cas9 mice empower a wide range of biological and disease modeling applications.

A transcriptomic taxonomy of mouse brain-wide spinal projecting neurons.

The brain controls nearly all bodily functions via spinal projecting neurons (SPNs) that carry command signals from the brain to the spinal cord. However, a comprehensive molecular characterization of brain-wide SPNs is still lacking. Here we transcriptionally profiled a total of 65,002 SPNs, identified 76 region-specific SPN types, and mapped these types into a companion atlas of the whole mouse brain1. This taxonomy reveals a three-component organization of SPNs: (1) molecularly homogeneous excitatory SPNs from the cortex, red nucleus and cerebellum with somatotopic spinal terminations suitable for point-to-point communication; (2) heterogeneous populations in the reticular formation with broad spinal termination patterns, suitable for relaying commands related to the activities of the entire spinal cord; and (3) modulatory neurons expressing slow-acting neurotransmitters and/or neuropeptides in the hypothalamus, midbrain and reticular formation for 'gain setting' of brain-spinal signals. In addition, this atlas revealed a LIM homeobox transcription factor code that parcellates the reticulospinal neurons into five molecularly distinct and spatially segregated populations. Finally, we found transcriptional signatures of a subset of SPNs with large soma size and correlated these with fast-firing electrophysiological properties. Together, this study establishes a comprehensive taxonomy of brain-wide SPNs and provides insight into the functional organization of SPNs in mediating brain control of bodily functions.

Conserved and divergent gene regulatory programs of the mammalian neocortex.

Divergence of cis-regulatory elements drives species-specific traits1, but how this manifests in the evolution of the neocortex at the molecular and cellular level remains unclear. Here we investigated the gene regulatory programs in the primary motor cortex of human, macaque, marmoset and mouse using single-cell multiomics assays, generating gene expression, chromatin accessibility, DNA methylome and chromosomal conformation profiles from a total of over 200,000 cells. From these data, we show evidence that divergence of transcription factor expression corresponds to species-specific epigenome landscapes. We find that conserved and divergent gene regulatory features are reflected in the evolution of the three-dimensional genome. Transposable elements contribute to nearly 80% of the human-specific candidate cis-regulatory elements in cortical cells. Through machine learning, we develop sequence-based predictors of candidate cis-regulatory elements in different species and demonstrate that the genomic regulatory syntax is highly preserved from rodents to primates. Finally, we show that epigenetic conservation combined with sequence similarity helps to uncover functional cis-regulatory elements and enhances our ability to interpret genetic variants contributing to neurological disease and traits.

Targeting thalamic circuits rescues motor and mood deficits in PD mice.

Although bradykinesia, tremor and rigidity are the hallmark motor defects in patients with Parkinson's disease (PD), patients also experience motor learning impairments and non-motor symptoms such as depression1. The neural circuit basis for these different symptoms of PD are not well understood. Although current treatments are effective for locomotion deficits in PD2,3, therapeutic strategies targeting motor learning deficits and non-motor symptoms are lacking4-6. Here we found that distinct parafascicular (PF) thalamic subpopulations project to caudate putamen (CPu), subthalamic nucleus (STN) and nucleus accumbens (NAc). Whereas PF→CPu and PF→STN circuits are critical for locomotion and motor learning, respectively, inhibition of the PF→NAc circuit induced a depression-like state. Whereas chemogenetically manipulating CPu-projecting PF neurons led to a long-term restoration of locomotion, optogenetic long-term potentiation (LTP) at PF→STN synapses restored motor learning behaviour in an acute mouse model of PD. Furthermore, activation of NAc-projecting PF neurons rescued depression-like phenotypes. Further, we identified nicotinic acetylcholine receptors capable of modulating PF circuits to rescue different PD phenotypes. Thus, targeting PF thalamic circuits may be an effective strategy for treating motor and non-motor deficits in PD.

The NIH Somatic Cell Genome Editing program.

The move from reading to writing the human genome offers new opportunities to improve human health. The United States National Institutes of Health (NIH) Somatic Cell Genome Editing (SCGE) Consortium aims to accelerate the development of safer and more-effective methods to edit the genomes of disease-relevant somatic cells in patients, even in tissues that are difficult to reach. Here we discuss the consortium's plans to develop and benchmark approaches to induce and measure genome modifications, and to define downstream functional consequences of genome editing within human cells. Central to this effort is a rigorous and innovative approach that requires validation of the technology through third-party testing in small and large animals. New genome editors, delivery technologies and methods for tracking edited cells in vivo, as well as newly developed animal models and human biological systems, will be assembled-along with validated datasets-into an SCGE Toolkit, which will be disseminated widely to the biomedical research community. We visualize this toolkit-and the knowledge generated by its applications-as a means to accelerate the clinical development of new therapies for a wide range of conditions.

Evolutionary and developmental specialization of foveal cell types in the marmoset.

In primates, high-acuity vision is mediated by the fovea, a small specialized central region of the retina. The fovea, unique to the anthropoid lineage among mammals, undergoes notable neuronal morphological changes during postnatal maturation. However, the extent of cellular similarity across anthropoid foveas and the molecular underpinnings of foveal maturation remain unclear. Here, we used high-throughput single-cell RNA sequencing to profile retinal cells of the common marmoset (Callithrix jacchus), an early divergent in anthropoid evolution from humans, apes, and macaques. We generated atlases of the marmoset fovea and peripheral retina for both neonates and adults. Our comparative analysis revealed that marmosets share almost all their foveal types with both humans and macaques, highlighting a conserved cellular structure among primate foveas. Furthermore, by tracing the developmental trajectory of cell types in the foveal and peripheral retina, we found distinct maturation paths for each. In-depth analysis of gene expression differences demonstrated that cone photoreceptors and Müller glia (MG), among others, show the greatest molecular divergence between these two regions. Utilizing single-cell ATAC-seq and gene-regulatory network inference, we uncovered distinct transcriptional regulations differentiating foveal cones from their peripheral counterparts. Further analysis of predicted ligand-receptor interactions suggested a potential role for MG in supporting the maturation of foveal cones. Together, these results provide valuable insights into foveal development, structure, and evolution.

Innovations present in the primate interneuron repertoire.

Primates and rodents, which descended from a common ancestor around 90 million years ago1, exhibit profound differences in behaviour and cognitive capacity; the cellular basis for these differences is unknown. Here we use single-nucleus RNA sequencing to profile RNA expression in 188,776 individual interneurons across homologous brain regions from three primates (human, macaque and marmoset), a rodent (mouse) and a weasel (ferret). Homologous interneuron types-which were readily identified by their RNA-expression patterns-varied in abundance and RNA expression among ferrets, mice and primates, but varied less among primates. Only a modest fraction of the genes identified as 'markers' of specific interneuron subtypes in any one species had this property in another species. In the primate neocortex, dozens of genes showed spatial expression gradients among interneurons of the same type, which suggests that regional variation in cortical contexts shapes the RNA expression patterns of adult neocortical interneurons. We found that an interneuron type that was previously associated with the mouse hippocampus-the 'ivy cell', which has neurogliaform characteristics-has become abundant across the neocortex of humans, macaques and marmosets but not mice or ferrets. We also found a notable subcortical innovation: an abundant striatal interneuron type in primates that had no molecularly homologous counterpart in mice or ferrets. These interneurons expressed a unique combination of genes that encode transcription factors, receptors and neuropeptides and constituted around 30% of striatal interneurons in marmosets and humans.

Distinct subnetworks of the thalamic reticular nucleus.

The thalamic reticular nucleus (TRN), the major source of thalamic inhibition, regulates thalamocortical interactions that are critical for sensory processing, attention and cognition1-5. TRN dysfunction has been linked to sensory abnormality, attention deficit and sleep disturbance across multiple neurodevelopmental disorders6-9. However, little is known about the organizational principles that underlie its divergent functions. Here we performed an integrative study linking single-cell molecular and electrophysiological features of the mouse TRN to connectivity and systems-level function. We found that cellular heterogeneity in the TRN is characterized by a transcriptomic gradient of two negatively correlated gene-expression profiles, each containing hundreds of genes. Neurons in the extremes of this transcriptomic gradient express mutually exclusive markers, exhibit core or shell-like anatomical structure and have distinct electrophysiological properties. The two TRN subpopulations make differential connections with the functionally distinct first-order and higher-order thalamic nuclei to form molecularly defined TRN-thalamus subnetworks. Selective perturbation of the two subnetworks in vivo revealed their differential role in regulating sleep. In sum, our study provides a comprehensive atlas of TRN neurons at single-cell resolution and links molecularly defined subnetworks to the functional organization of thalamocortical circuits.

Author Correction: Innovations present in the primate interneuron repertoire.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Multi-animal pose estimation, identification and tracking with DeepLabCut.

Estimating the pose of multiple animals is a challenging computer vision problem: frequent interactions cause occlusions and complicate the association of detected keypoints to the correct individuals, as well as having highly similar looking animals that interact more closely than in typical multi-human scenarios. To take up this challenge, we build on DeepLabCut, an open-source pose estimation toolbox, and provide high-performance animal assembly and tracking-features required for multi-animal scenarios. Furthermore, we integrate the ability to predict an animal's identity to assist tracking (in case of occlusions). We illustrate the power of this framework with four datasets varying in complexity, which we release to serve as a benchmark for future algorithm development.

Atypical behaviour and connectivity in SHANK3-mutant macaques.

Mutation or disruption of the SH3 and ankyrin repeat domains 3 (SHANK3) gene represents a highly penetrant, monogenic risk factor for autism spectrum disorder, and is a cause of Phelan-McDermid syndrome. Recent advances in gene editing have enabled the creation of genetically engineered non-human-primate models, which might better approximate the behavioural and neural phenotypes of autism spectrum disorder than do rodent models, and may lead to more effective treatments. Here we report CRISPR-Cas9-mediated generation of germline-transmissible mutations of SHANK3 in cynomolgus macaques (Macaca fascicularis) and their F1 offspring. Genotyping of somatic cells as well as brain biopsies confirmed mutations in the SHANK3 gene and reduced levels of SHANK3 protein in these macaques. Analysis of data from functional magnetic resonance imaging revealed altered local and global connectivity patterns that were indicative of circuit abnormalities. The founder mutants exhibited sleep disturbances, motor deficits and increased repetitive behaviours, as well as social and learning impairments. Together, these results parallel some aspects of the dysfunctions in the SHANK3 gene and circuits, as well as the behavioural phenotypes, that characterize autism spectrum disorder and Phelan-McDermid syndrome.

Anterior thalamic circuits crucial for working memory.

Alterations in the structure and functional connectivity of anterior thalamic nuclei (ATN) have been linked to reduced cognition during aging. However, ATN circuits that contribute to higher cognitive functions remain understudied. We found that the anteroventral (AV) subdivision of ATN is necessary specifically during the maintenance phase of a spatial working memory task. This function engages the AV→parasubiculum (PaS)→entorhinal cortex (EC) circuit. Aged mice showed a deficit in spatial working memory, which was associated with a decrease in the excitability of AV neurons. Activation of AV neurons or the AV→PaS circuit in aged mice was sufficient to rescue their working memory performance. Furthermore, rescued aged mice showed improved behavior-induced neuronal activity in prefrontal cortex (PFC), a critical site for working memory processes. Although the direct activation of PFC neurons in aged mice also rescued their working memory performance, we found that these animals exhibited increased levels of anxiety, which was not the case for AV→PaS circuit manipulations in aged mice. These results suggest that targeting AV thalamus in aging may not only be beneficial for cognitive functions but that this approach may have fewer unintended effects compared to direct PFC manipulations.

Sapap4 deficiency leads to postsynaptic defects and abnormal behaviors relevant to hyperkinetic neuropsychiatric disorder in mice.

Postsynaptic proteins play critical roles in synaptic development, function, and plasticity. Dysfunction of postsynaptic proteins is strongly linked to neurodevelopmental and psychiatric disorders. SAP90/PSD95-associated protein 4 (SAPAP4; also known as DLGAP4) is a key component of the PSD95-SAPAP-SHANK excitatory postsynaptic scaffolding complex, which plays important roles at synapses. However, the exact function of the SAPAP4 protein in the brain is poorly understood. Here, we report that Sapap4 knockout (KO) mice have reduced spine density in the prefrontal cortex and abnormal compositions of key postsynaptic proteins in the postsynaptic density (PSD) including reduced PSD95, GluR1, and GluR2 as well as increased SHANK3. These synaptic defects are accompanied by a cluster of abnormal behaviors including hyperactivity, impulsivity, reduced despair/depression-like behavior, hypersensitivity to low dose of amphetamine, memory deficits, and decreased prepulse inhibition, which are reminiscent of mania. Furthermore, the hyperactivity of Sapap4 KO mice could be partially rescued by valproate, a mood stabilizer used for mania treatment in humans. Together, our findings provide evidence that SAPAP4 plays an important role at synapses and reinforce the view that dysfunction of the postsynaptic scaffolding protein SAPAP4 may contribute to the pathogenesis of hyperkinetic neuropsychiatric disorder.

Opportunities and limitations of genetically modified nonhuman primate models for neuroscience research.

The recently developed new genome-editing technologies, such as the CRISPR/Cas system, have opened the door for generating genetically modified nonhuman primate (NHP) models for basic neuroscience and brain disorders research. The complex circuit formation and experience-dependent refinement of the human brain are very difficult to model in vitro, and thus require use of in vivo whole-animal models. For many neurodevelopmental and psychiatric disorders, abnormal circuit formation and refinement might be at the center of their pathophysiology. Importantly, many of the critical circuits and regional cell populations implicated in higher human cognitive function and in many psychiatric disorders are not present in lower mammalian brains, while these analogous areas are replicated in NHP brains. Indeed, neuropsychiatric disorders represent a tremendous health and economic burden globally. The emerging field of genetically modified NHP models has the potential to transform our study of higher brain function and dramatically facilitate the development of effective treatment for human brain disorders. In this paper, we discuss the importance of developing such models, the infrastructure and training needed to maximize the impact of such models, and ethical standards required for using these models.

SHANK proteins: roles at the synapse and in autism spectrum disorder.

Several large-scale genomic studies have supported an association between cases of autism spectrum disorder and mutations in the genes SH3 and multiple ankyrin repeat domains protein 1 (SHANK1), SHANK2 and SHANK3, which encode a family of postsynaptic scaffolding proteins that are present at glutamatergic synapses in the CNS. An evaluation of human genetic data, as well as of in vitro and in vivo animal model data, may allow us to understand how disruption of SHANK scaffolding proteins affects the structure and function of neural circuits and alters behaviour.

Thalamic reticular impairment underlies attention deficit in Ptchd1(Y/-) mice.

Developmental disabilities, including attention-deficit hyperactivity disorder (ADHD), intellectual disability (ID), and autism spectrum disorders (ASD), affect one in six children in the USA. Recently, gene mutations in patched domain containing 1 (PTCHD1) have been found in ~1% of patients with ID and ASD. Individuals with PTCHD1 deletion show symptoms of ADHD, sleep disruption, hypotonia, aggression, ASD, and ID. Although PTCHD1 is probably critical for normal development, the connection between its deletion and the ensuing behavioural defects is poorly understood. Here we report that during early post-natal development, mouse Ptchd1 is selectively expressed in the thalamic reticular nucleus (TRN), a group of GABAergic neurons that regulate thalamocortical transmission, sleep rhythms, and attention. Ptchd1 deletion attenuates TRN activity through mechanisms involving small conductance calcium-dependent potassium currents (SK). TRN-restricted deletion of Ptchd1 leads to attention deficits and hyperactivity, both of which are rescued by pharmacological augmentation of SK channel activity. Global Ptchd1 deletion recapitulates learning impairment, hyper-aggression, and motor defects, all of which are insensitive to SK pharmacological targeting and not found in the TRN-restricted deletion mouse. This study maps clinically relevant behavioural phenotypes onto TRN dysfunction in a human disease model, while also identifying molecular and circuit targets for intervention.