Transcriptomic and open chromatin atlas of high-resolution anatomical regions in the rhesus macaque brain

Xenopus tropicalis osteoblast-specific open chromatin regions reveal promoters and enhancers involved in human skeletal phenotypes and shed light on early vertebrate evolution

While understanding the genetic underpinnings of osteogenesis has far-reaching implications for skeletal diseases and evolution, a comprehensive characterization of the osteoblastic regulatory landscape in non-mammalian vertebrates is still lacking. Here, we compared the ATAC-Seq profile of Xenopus tropicalis (Xt) osteoblasts to a variety of non mineralizing control tissues, and identified osteoblast-specific nucleosome free regions (NFRs) at 527 promoters and 6747 distal regions. Sequence analyses, Gene Ontology, RNA-Seq and ChIP-Seq against four key histone marks confirmed that the distal regions correspond to bona fide osteogenic transcriptional enhancers exhibiting a shared regulatory logic with mammals. We report 425 regulatory regions conserved with human and globally associated to skeletogenic genes. Of these, 35 regions have been shown to impact human skeletal phenotypes by GWAS, including one trps1 enhancer and the runx2 promoter, two genes which are respectively involved in trichorhinophalangeal syndrome type I and cleidocranial dysplasia. Intriguingly, 60 osteoblastic NFRs also align to the genome of the elephant shark, a species lacking osteoblasts and bone tissue. To tackle this paradox, we chose to focus on dlx5 because its conserved promoter, known to integrate regulatory inputs during mammalian osteogenesis, harbours an osteoblast-specific NFR in both frog and human. Hence, we show that dlx5 is expressed in Xt and elephant shark odontoblasts, supporting a common cellular and genetic origin of bone and dentine. Taken together, our work (i) unravels the Xt osteogenic regulatory landscape, (ii) illustrates how cross-species comparisons harvest data relevant to human biology and (iii) reveals that a set of genes including bnc2, dlx5, ebf3, mir199a, nfia, runx2 and zfhx4 drove the development of a primitive form of mineralized skeletal tissue deep in the vertebrate lineage.

Reciprocal negative feedback between Prrx1 and miR-140-3p regulates rapid chondrogenesis in the regenerating antler

During growth phase, antlers exhibit a very rapid rate of chondrogenesis. The antler is formed from its growth center reserve mesenchyme (RM) cells, which have been found to be the derivatives of paired related homeobox 1 (Prrx1)-positive periosteal cells. However, the underlying mechanism that drives rapid chondrogenesis is not known. Herein, the miRNA expression profiles and chromatin states of three tissue layers (RM, precartilage, and cartilage) at different stages of differentiation within the antler growth center were analyzed by RNA-sequencing and ATAC-sequencing. We found that miR-140-3p was the miRNA that exhibited the greatest degree of upregulation in the rapidly growing antler, increasing from the RM to the cartilage layer. We also showed that Prrx1 was a key upstream regulator of miR-140-3p, which firmly confirmed by Prrx1 CUT&Tag sequencing of RM cells. Through multiple approaches (three-dimensional chondrogenic culture and xenogeneic antler model), we demonstrated that Prrx1 and miR-140-3p functioned as reciprocal negative feedback in the antler growth center, and downregulating PRRX1/upregulating miR-140-3p promoted rapid chondrogenesis of RM cells and xenogeneic antler. Thus, we conclude that the reciprocal negative feedback between Prrx1 and miR-140-3p is essential for balancing mesenchymal proliferation and chondrogenic differentiation in the regenerating antler. We further propose that the mechanism underlying chondrogenesis in the regenerating antler would provide a reference for helping understand the regulation of human cartilage regeneration and repair.

Multimodal Nature of the Single-cell Primate Brain Atlas: Morphology, Transcriptome, Electrophysiology, and Connectivity

Primates exhibit complex brain structures that augment cognitive function. The neocortex fulfills high-cognitive functions through billions of connected neurons. These neurons have distinct transcriptomic, morphological, and electrophysiological properties, and their connectivity principles vary. These features endow the primate brain atlas with a multimodal nature. The recent integration of next-generation sequencing with modified patch-clamp techniques is revolutionizing the way to census the primate neocortex, enabling a multimodal neuronal atlas to be established in great detail: (1) single-cell/single-nucleus RNA-seq technology establishes high-throughput transcriptomic references, covering all major transcriptomic cell types; (2) patch-seq links the morphological and electrophysiological features to the transcriptomic reference; (3) multicell patch-clamp delineates the principles of local connectivity. Here, we review the applications of these technologies in the primate neocortex and discuss the current advances and tentative gaps for a comprehensive understanding of the primate neocortex.

Integrated multi-omics analysis of brain aging in female nonhuman primates reveals altered signaling pathways relevant to age-related disorders

The prefrontal cortex (PFC) has been implicated as a key brain region responsible for age-related cognitive decline. Little is known about aging-related molecular changes in PFC that may mediate these effects. To date, no studies have used untargeted discovery methods with integrated analyses to determine PFC molecular changes in healthy female primates. We quantified PFC changes associated with healthy aging in female baboons by integrating multiple omics data types (transcriptomics, proteomics, metabolomics) from samples across the adult age span. Our integrated omics approach using unbiased weighted gene co-expression network analysis to integrate data and treat age as a continuous variable, revealed highly interconnected known and novel pathways associated with PFC aging. We found Gamma-aminobutyric acid (GABA) tissue content associated with these signaling pathways, providing 1 potential biomarker to assess PFC changes with age. These highly coordinated pathway changes during aging may represent early steps for aging-related decline in PFC functions, such as learning and memory, and provide potential biomarkers to assess cognitive status in humans.

The molecular genetic landscape of human brain size variation

Human brain size changes dynamically through early development, peaks in adolescence, and varies up to 2-fold among adults. However, the molecular genetic underpinnings of interindividual variation in brain size remain unknown. Here, we leveraged postmortem brain RNA sequencing and measurements of brain weight (BW) in 2,531 individuals across three independent datasets to identify 928 genome-wide significant associations with BW. Genes associated with higher or lower BW showed distinct neurodevelopmental trajectories and spatial patterns that mapped onto functional and cellular axes of brain organization. Expression of BW genes was predictive of interspecies differences in brain size, and bioinformatic annotation revealed enrichment for neurogenesis and cell-cell communication. Genome-wide, transcriptome-wide, and phenome-wide association analyses linked BW gene sets to neuroimaging measurements of brain size and brain-related clinical traits. Cumulatively, these results represent a major step toward delineating the molecular pathways underlying human brain size variation in health and disease.

Sex differences in brain transcriptomes of juvenile Cynomolgus macaques

Background: Behavioral, social, and physical characteristics are posited to distinguish the sexes, yet research on transcription-level sexual differences in the brain is limited. Here, we investigated sexually divergent brain transcriptomics in prepubertal cynomolgus macaques, a commonly used surrogate species to humans. Methods: A transcriptomic profile using RNA sequencing was generated for the temporal lobe, ventral midbrain, and cerebellum of 3 female and 3 male cynomolgus macaques previously treated with an Adeno-associated virus vector mix. Statistical analyses to determine differentially expressed protein-coding genes in all three lobes were conducted using DeSeq2 with a false discovery rate corrected P value of .05. Results: We identified target genes in the temporal lobe, ventral midbrain, and cerebellum with functions in translation, immunity, behavior, and neurological disorders that exhibited statistically significant sexually divergent expression. Conclusions: We provide potential mechanistic insights to the epidemiological differences observed between the sexes with regards to mental health and infectious diseases, such as COVID19. Our results provide pre-pubertal information on sexual differences in non-human primate brain transcriptomics and may provide insight to health disparities between the biological sexes in humans.

Genetic Approaches for Neural Circuits Dissection in Non-human Primates

Genetic tools, which can be used for the morphology study of specific neurons, pathway-selective connectome mapping, neuronal activity monitoring, and manipulation with a spatiotemporal resolution, have been widely applied to the understanding of complex neural circuit formation, interactions, and functions in rodents. Recently, similar genetic approaches have been tried in non-human primates (NHPs) in neuroscience studies for dissecting the neural circuits involved in sophisticated behaviors and clinical brain disorders, although they are still very preliminary. In this review, we introduce the progress made in the development and application of genetic tools for brain studies on NHPs. We also discuss the advantages and limitations of each approach and provide a perspective for using genetic tools to study the neural circuits of NHPs.

Spatiotemporal proteomic atlas of multiple brain regions across early fetal to neonatal stages in cynomolgus monkey

Fetal stages are critical periods for brain development. However, the protein molecular signature and dynamics of the human brain remain unclear due to sampling difficulty and ethical limitations. Non-human primates present similar developmental and neuropathological features to humans. This study constructed a spatiotemporal proteomic atlas of cynomolgus macaque brain development from early fetal to neonatal stages. Here we showed that (1) the variability across stages was greater than that among brain regions, and comparisons of cerebellum vs. cerebrum and cortical vs. subcortical regions revealed region-specific dynamics across early fetal to neonatal stages; (2) fluctuations in abundance of proteins associated with neural disease suggest the risk of nervous disorder at early fetal stages; (3) cross-species analysis (human, monkey, and mouse) and comparison between proteomic and transcriptomic data reveal the proteomic specificity and genes with mRNA/protein discrepancy. This study provides insight into fetal brain development in primates.

Cynomolgus-rhesus hybrid macaques serve as a platform for imprinting studies

Genomic imprinting can lead to allele-specific expression (ASE), where one allele is preferentially expressed more than the other. Perturbations in genomic imprinting or ASE genes have been widely observed across various neurological disorders, notably autism spectrum disorder (ASD). In this study, we crossed rhesus cynomolgus monkeys to produce hybrid monkeys and established a framework to evaluate their allele-specific gene expression patterns using the parental genomes as a reference. Our proof-of-concept analysis of the hybrid monkeys identified 353 genes with allele-biased expression in the brain, enabling us to determine the chromosomal locations of ASE clusters. Importantly, we confirmed a significant enrichment of ASE genes associated with neuropsychiatric disorders, including ASD, highlighting the potential of hybrid monkey models in advancing our understanding of genomic imprinting.

Neuroinflammatory transcriptional programs induced in rhesus pre-frontal cortex white matter during acute SHIV infection

Background

Immunosurveillance of the central nervous system (CNS) is vital to resolve infection and injury. However, immune activation within the CNS in the setting of chronic viral infections, such as HIV-1, is strongly linked to progressive neurodegeneration and cognitive decline. Establishment of HIV-1 in the CNS early following infection underscores the need to delineate features of acute CNS immune activation, as these early inflammatory events may mediate neurodegenerative processes. Here, we focused on elucidating molecular programs of neuroinflammation in brain regions based on vulnerability to neuroAIDS and/or neurocognitive decline. To this end, we assessed transcriptional profiles within the subcortical white matter of the pre-frontal cortex (PFCw), as well as synapse dense regions from hippocampus, superior temporal cortex, and caudate nucleus, in rhesus macaques following infection with Simian/Human Immunodeficiency Virus (SHIV.C.CH505).

Methods

We performed RNA extraction and sequenced RNA isolated from 3 mm brain punches. Viral RNA was quantified in the brain and cerebrospinal fluid by RT-qPCR assays targeting SIV Gag. Neuroinflammation was assessed by flow cytometry and multiplex ELISA assays.

Results

RNA sequencing and flow cytometry data demonstrated immune surveillance of the rhesus CNS by innate and adaptive immune cells during homeostasis. Following SHIV infection, viral entry and integration within multiple brain regions demonstrated vulnerabilities of key cognitive and motor function brain regions to HIV-1 during the acute phase of infection. SHIV-induced transcriptional alterations were concentrated to the PFCw and STS with upregulation of gene expression pathways controlling innate and T-cell inflammatory responses. Within the PFCw, gene modules regulating microglial activation and T cell differentiation were induced at 28 days post-SHIV infection, with evidence for stimulation of immune effector programs characteristic of neuroinflammation. Furthermore, enrichment of pathways regulating mitochondrial respiratory capacity, synapse assembly, and oxidative and endoplasmic reticulum stress were observed. These acute neuroinflammatory features were substantiated by increased influx of activated T cells into the CNS.

Conclusions

Our data show pervasive immune surveillance of the rhesus CNS at homeostasis and reveal perturbations of important immune, neuronal, and synaptic pathways within key anatomic regions controlling cognition and motor function during acute HIV infection. These findings provide a valuable framework to understand early molecular features of HIV associated neurodegeneration.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12974-022-02610-y.

Somatic retrotransposition in the developing rhesus macaque brain

The retrotransposon LINE-1 (L1) is central to the recent evolutionary history of the human genome and continues to drive genetic diversity and germline pathogenesis. However, the spatiotemporal extent and biological significance of somatic L1 activity are poorly defined and are virtually unexplored in other primates. From a single L1 lineage active at the divergence of apes and Old World monkeys, successive L1 subfamilies have emerged in each descendant primate germline. As revealed by case studies, the presently active human L1 subfamily can also mobilize during embryonic and brain development in vivo. It is unknown whether nonhuman primate L1s can similarly generate somatic insertions in the brain. Here we applied approximately 40× single-cell whole-genome sequencing (scWGS), as well as retrotransposon capture sequencing (RC-seq), to 20 hippocampal neurons from two rhesus macaques (Macaca mulatta). In one animal, we detected and PCR-validated a somatic L1 insertion that generated target site duplications, carried a short 5′ transduction, and was present in ∼7% of hippocampal neurons but absent from cerebellum and nonbrain tissues. The corresponding donor L1 allele was exceptionally mobile in vitro and was embedded in PRDM4, a gene expressed throughout development and in neural stem cells. Nanopore long-read methylome and RNA-seq transcriptome analyses indicated young retrotransposon subfamily activation in the early embryo, followed by repression in adult tissues. These data highlight endogenous macaque L1 retrotransposition potential, provide prototypical evidence of L1-mediated somatic mosaicism in a nonhuman primate, and allude to L1 mobility in the brain over the past 30 million years of human evolution.

Layer-specific pyramidal neuron properties underlie diverse anterior cingulate cortical motor and limbic networks

The laminar cellular and circuit mechanisms by which the anterior cingulate cortex (ACC) exerts flexible control of motor and affective information for goal-directed behavior have not been elucidated. Using multimodal tract-tracing, in vitro patch-clamp recording and computational approaches in rhesus monkeys (M. mulatta), we provide evidence that specialized motor and affective network dynamics can be conferred by layer-specific biophysical and structural properties of ACC pyramidal neurons targeting two key downstream structures -the dorsal premotor cortex (PMd) and the amygdala (AMY). AMY-targeting neurons exhibited significant laminar differences, with L5 more excitable (higher input resistance and action potential firing rates) than L3 neurons. Between-pathway differences were found within L5, with AMY-targeting neurons exhibiting greater excitability, apical dendritic complexity, spine densities, and diversity of inhibitory inputs than PMd-targeting neurons. Simulations using a pyramidal-interneuron network model predict that these layer- and pathway-specific single-cell differences contribute to distinct network oscillatory dynamics. L5 AMY-targeting networks are more tuned to slow oscillations well-suited for affective and contextual processing timescales, while PMd-targeting networks showed strong beta/gamma synchrony implicated in rapid sensorimotor processing. These findings are fundamental to our broad understanding of how layer-specific cellular and circuit properties can drive diverse laminar activity found in flexible behavior.

LncRNA LINC00281/Annexin A2 Regulates Vascular Smooth Muscle Cell Phenotype Switching via the Nuclear Factor-Kappa B Signaling Pathway

Abnormal phenotype switch in vascular smooth muscle cells (VSMCs) plays an important role in the initiation and progression of vascular proliferative diseases. Annexin A2 (ANXA2), related to the pro-inflammatory response, contributes to the proliferation and migration of VSMCs. This study explored the mechanisms involved in the regulation of VSMC phenotype modulation via ANXA2. The results revealed that ANXA2 promotes the phosphorylation of p65 and co-translocates with p65 into the nucleus, resulting in VSMC proliferation, migration, and dedifferentiation. Based on bioinformatics predictions and RNA immunoprecipitation assays, LINC00281 was confirmed to be an upstream regulator of ANXA2. Taken together, ANXA2, which is negatively regulated by the long noncoding RNA (lncRNA) LINC00281, has significant importance in the regulation of VSMC proliferation, migration, and phenotype switching via the nuclear factor-kappa B (NF-кB) p65 signaling pathway. This indicates that the lncRNA LINC00281/ANXA2/NF-кB p65 signaling pathway might be a new therapeutic target for vascular proliferative diseases.

Comprehensive enhancer-target gene assignments improve gene set level interpretation of genome-wide regulatory data

BACKGROUND

Revealing the gene targets of distal regulatory elements is challenging yet critical for interpreting regulome data. Experiment-derived enhancer-gene links are restricted to a small set of enhancers and/or cell types, while the accuracy of genome-wide approaches remains elusive due to the lack of a systematic evaluation. We combined multiple spatial and in silico approaches for defining enhancer locations and linking them to their target genes aggregated across >500 cell types, generating 1860 human genome-wide distal enhancer-to-target gene definitions (EnTDefs). To evaluate performance, we used gene set enrichment (GSE) testing on 87 independent ENCODE ChIP-seq datasets of 34 transcription factors (TFs) and assessed concordance of results with known TF Gene Ontology annotations, and other benchmarks.

RESULTS

The top ranked 741 (40%) EnTDefs significantly outperform the common, naïve approach of linking distal regions to the nearest genes, and the top 10 EnTDefs perform well when applied to ChIP-seq data of other cell types. The GSE-based ranking of EnTDefs is highly concordant with ranking based on overlap with curated benchmarks of enhancer-gene interactions. Both our top general EnTDef and cell-type-specific EnTDefs significantly outperform seven independent computational and experiment-based enhancer-gene pair datasets. We show that using our top EnTDefs for GSE with either genome-wide DNA methylation or ATAC-seq data is able to better recapitulate the biological processes changed in gene expression data performed in parallel for the same experiment than our lower-ranked EnTDefs.

CONCLUSIONS

Our findings illustrate the power of our approach to provide genome-wide interpretation regardless of cell type.

A triple-network organization for the mouse brain

The triple-network model of psychopathology is a framework to explain the functional and structural neuroimaging phenotypes of psychiatric and neurological disorders. It describes the interactions within and between three distributed networks: the salience, default-mode, and central executive networks. These have been associated with brain disorder traits in patients. Homologous networks have been proposed in animal models, but their integration into a triple-network organization has not yet been determined. Using resting-state datasets, we demonstrate conserved spatio-temporal properties between triple-network elements in human, macaque, and mouse. The model predictions were also shown to apply in a mouse model for depression. To validate spatial homologies, we developed a data-driven approach to convert mouse brain maps into human standard coordinates. Finally, using high-resolution viral tracers in the mouse, we refined an anatomical model for these networks and validated this using optogenetics in mice and tractography in humans. Unexpectedly, we find serotonin involvement within the salience rather than the default-mode network. Our results support the existence of a triple-network system in the mouse that shares properties with that of humans along several dimensions, including a disease condition. Finally, we demonstrate a method to humanize mouse brain networks that opens doors to fully data-driven trans-species comparisons.

Transcriptional and anatomical diversity of medium spiny neurons in the primate striatum

Medium spiny neurons (MSNs) constitute the vast majority of striatal neurons and the principal interface between dopamine reward signals and functionally diverse cortico-basal ganglia circuits. Information processing in these circuits is dependent on distinct MSN types: cell types that are traditionally defined according to their projection targets or dopamine receptor expression. Single-cell transcriptional studies have revealed greater MSN heterogeneity than predicted by traditional circuit models, but the transcriptional landscape in the primate striatum remains unknown. Here, we set out to establish molecular definitions for MSN subtypes in Rhesus monkeys and to explore the relationships between transcriptionally defined subtypes and anatomical subdivisions of the striatum. Our results suggest at least nine MSN subtypes, including dorsal striatum subtypes associated with striosome and matrix compartments, ventral striatum subtypes associated with the nucleus accumbens shell and olfactory tubercle, and an MSN-like cell type restricted to μ-opioid receptor rich islands in the ventral striatum. Although each subtype was demarcated by discontinuities in gene expression, continuous variation within subtypes defined gradients corresponding to anatomical locations and, potentially, functional specializations. These results lay the foundation for achieving cell-type-specific transgenesis in the primate striatum and provide a blueprint for investigating circuit-specific information processing.

clusterProfiler 4.0: A universal enrichment tool for interpreting omics data

Functional enrichment analysis is pivotal for interpreting high-throughput omics data in life science. It is crucial for this type of tool to use the latest annotation databases for as many organisms as possible. To meet these requirements, we present here an updated version of our popular Bioconductor package, clusterProfiler 4.0. This package has been enhanced considerably compared with its original version published 9 years ago. The new version provides a universal interface for functional enrichment analysis in thousands of organisms based on internally supported ontologies and pathways as well as annotation data provided by users or derived from online databases. It also extends the dplyr and ggplot2 packages to offer tidy interfaces for data operation and visualization. Other new features include gene set enrichment analysis and comparison of enrichment results from multiple gene lists. We anticipate that clusterProfiler 4.0 will be applied to a wide range of scenarios across diverse organisms.

•

clusterProfiler supports exploring functional characteristics of both coding and non-coding genomics data for thousands of species with up-to-date gene annotation

•

It provides a universal interface for gene functional annotation from a variety of sources and thus can be applied in diverse scenarios

•

It provides a tidy interface to access, manipulate, and visualize enrichment results to help users achieve efficient data interpretation

•

Datasets obtained from multiple treatments and time points can be analyzed and compared in a single run, easily revealing functional consensus and differences among distinct conditions

Evolutionary Changes in Pathways and Networks of Genes Expressed in the Brains of Humans and Macaques

As the key organ that separates humans from nonhuman primates, the brain has continuously evolved to adapt to environmental and climatic changes. Although humans share most genetic, molecular, and cellular features with other primates such as macaques, there are significant differences in the structure and function of the brain between humans and these species. Thus, exploring the differences between the brains of human and nonhuman primates in the context of evolution will provide insights into the development, functionality, and diseases of the human central nervous system (CNS). Since the genes involved in many aspects of the human brain are under common pressures of natural selection, their evolutionary features can be analyzed collectively at the pathway level. In this study, the molecular mechanisms underlying human brain capabilities were explored by comparing the evolution features of pathways enriched in genes expressed in the human brain and the macaque brain. We identified 31 pathways with differential evolutionary properties, including those related to neurological diseases, signal transduction, immunological response, and metabolic processes. By analyzing genes differentially expressed in brain regions or development stages between humans and macaques,  9 and 4 pathways with differential evolutionary properties were detected, respectively. We further performed crosstalk analysis on the pathways to obtain an intuitive correlation between the pathways, which is helpful in understanding the mechanisms of interaction between pathways. Our results provide on a comprehensive view of the evolutionary pathways of the human CNS and can serve as a reference for the study of human brain development.

Different Flavors of Astrocytes: Revising the Origins of Astrocyte Diversity and Epigenetic Signatures to Understand Heterogeneity after Injury

Astrocytes are a specific type of neuroglial cells that confer metabolic and structural support to neurons. Astrocytes populate all regions of the nervous system and adopt a variety of phenotypes depending on their location and their respective functions, which are also pleiotropic in nature. For example, astrocytes adapt to pathological conditions with a specific cellular response known as reactive astrogliosis, which includes extensive phenotypic and transcriptional changes. Reactive astrocytes may lose some of their homeostatic functions and gain protective or detrimental properties with great impact on damage propagation. Different astrocyte subpopulations seemingly coexist in reactive astrogliosis, however, the source of such heterogeneity is not completely understood. Altered cellular signaling in pathological compared to healthy conditions might be one source fueling astrocyte heterogeneity. Moreover, diversity might also be encoded cell-autonomously, for example as a result of astrocyte subtype specification during development. We hypothesize and propose here that elucidating the epigenetic signature underlying the phenotype of each astrocyte subtype is of high relevance to understand another regulative layer of astrocyte heterogeneity, in general as well as after injury or as a result of other pathological conditions. High resolution methods should allow enlightening diverse cell states and subtypes of astrocyte, their adaptation to pathological conditions and ultimately allow controlling and manipulating astrocyte functions in disease states. Here, we review novel literature reporting on astrocyte diversity from a developmental perspective and we focus on epigenetic signatures that might account for cell type specification.

Chromatin accessibility in neuropsychiatric disorders

Epigenetic mechanisms have recently emerged as critical regulators of brain function in health and disease. By controlling the accessibility and the expression of specific genes, these pathways can mediate transient and long-lasting changes in neuronal function in both physiological and pathological contexts. Due to the extreme complexity of the epigenetic regulatory landscape, emerging methods that directly assay chromatin accessibility are of particular interest. Here, I review recent insights gained on open and closed chromatin states in the brain, with emphasis on neuropsychiatric disorders. These advances generated an invaluable wealth of information that can help us better understand gene regulation in the brain and its impairments that contribute to the development of disease.

Brain transcriptomics of nonhuman primates: A review

The brain is one of the most important and intricate organs in our bodies. Interpreting brain function and illustrating the changes and molecular mechanisms during physiological or pathological processes are essential but sometimes difficult to achieve. In addition to histology, ethology and pharmacology, the development of transcriptomics alleviates this condition by enabling high-throughput observation of the brain at various levels of anatomical specificity. Moreover, because human brain samples are scarce, the brains of nonhuman primates are important alternative models. Here in this review, we summarize the applications of transcriptomics in nonhuman primate brain studies, including investigations of brain development, aging, toxic effects and diseases. Overall, as a powerful tool with developmental potential, transcriptomics has been widely utilized in neuroscience.