Genetic diversity of 1,845 rhesus macaques improves genetic variation interpretation and identifies disease models.

Understanding and treating human diseases require valid animal models. Leveraging the genetic diversity in rhesus macaque populations across eight primate centers in the United States, we conduct targeted-sequencing on 1845 individuals for 374 genes linked to inherited human retinal and neurodevelopmental diseases. We identify over 47,000 single nucleotide variants, a substantial proportion of which are shared with human populations. By combining rhesus and human allele frequencies with established variant prediction methods, we develop a machine learning-based score that outperforms established methods in predicting missense variant pathogenicity. Remarkably, we find a marked number of loss-of-function variants and putative deleterious variants, which may lead to the development of rhesus disease models. Through phenotyping of macaques carrying a pathogenic OPA1:p.A8S variant, we identify a genetic model of autosomal dominant optic atrophy. Finally, we present a public website housing variant and genotype data from over two thousand rhesus macaques.

Factorized visual representations in the primate visual system and deep neural networks.

Object classification has been proposed as a principal objective of the primate ventral visual stream and has been used as an optimization target for deep neural network models (DNNs) of the visual system. However, visual brain areas represent many different types of information, and optimizing for classification of object identity alone does not constrain how other information may be encoded in visual representations. Information about different scene parameters may be discarded altogether ('invariance'), represented in non-interfering subspaces of population activity ('factorization') or encoded in an entangled fashion. In this work, we provide evidence that factorization is a normative principle of biological visual representations. In the monkey ventral visual hierarchy, we found that factorization of object pose and background information from object identity increased in higher-level regions and strongly contributed to improving object identity decoding performance. We then conducted a large-scale analysis of factorization of individual scene parameters - lighting, background, camera viewpoint, and object pose - in a diverse library of DNN models of the visual system. Models which best matched neural, fMRI, and behavioral data from both monkeys and humans across 12 datasets tended to be those which factorized scene parameters most strongly. Notably, invariance to these parameters was not as consistently associated with matches to neural and behavioral data, suggesting that maintaining non-class information in factorized activity subspaces is often preferred to dropping it altogether. Thus, we propose that factorization of visual scene information is a widely used strategy in brains and DNN models thereof.

When looking at a picture, we can quickly identify a recognizable object, such as an apple, applying a single word label to it. Although extensive neuroscience research has focused on how human and monkey brains achieve this recognition, our understanding of how the brain and brain-like computer models interpret other complex aspects of a visual scene ­ such as object position and environmental context ­ remains incomplete. In particular, it was not clear to what extent object recognition comes at the expense of other important scene details. For example, various aspects of the scene might be processed simultaneously. On the other hand, general object recognition may interfere with processing of such details. To investigate this, Lindsey and Issa analyzed 12 monkey and human brain datasets, as well as numerous computer models, to explore how different aspects of a scene are encoded in neurons and how these aspects are represented by computational models. The analysis revealed that preventing effective separation and retention of information about object pose and environmental context worsened object identification in monkey cortex neurons. In addition, the computer models that were the most brain-like could independently preserve the other scene details without interfering with object identification. The findings suggest that human and monkey high level ventral visual processing systems are capable of representing the environment in a more complex way than previously appreciated. In the future, studying more brain activity data could help to identify how rich the encoded information is and how it might support other functions like spatial navigation. This knowledge could help to build computational models that process the information in the same way, potentially improving their understanding of real-world scenes.

Monkey see, monkey do: Tracking iPS-cardiomyocyte survival and maturation in autografts.

Induced pluripotent stem cell-derived cardiomyocyte (iPSC-CM) therapy has emerged as a highly promising field of heart repair. Lin et al.1 presented compelling evidence on the long-term engraftment and maturation of autologous iPSC-CMs in two rhesus macaques, demonstrating unprecedented cardiac autografting data in large animal models without the need of immunosuppressants.

An exploratory deep learning approach to investigate tuberculosis pathogenesis in nonhuman primate model: Combining automated radiological analysis with clinical and biomarkers data.

BACKGROUND: Tuberculosis (TB) kills approximately 1.6 million people yearly despite the fact anti-TB drugs are generally curative. Therefore, TB-case detection and monitoring of therapy, need a comprehensive approach. Automated radiological analysis, combined with clinical, microbiological, and immunological data, by machine learning (ML), can help achieve it. METHODS: Six rhesus macaques were experimentally inoculated with pathogenic Mycobacterium tuberculosis in the lung. Data, including Computed Tomography (CT), were collected at 0, 2, 4, 8, 12, 16, and 20 weeks. RESULTS: Our ML-based CT analysis (TB-Net) efficiently and accurately analyzed disease progression, performing better than standard deep learning model (LLM OpenAI's CLIP Vi4). TB-Net based results were more consistent than, and confirmed independently by, blinded manual disease scoring by two radiologists and exhibited strong correlations with blood biomarkers, TB-lesion volumes, and disease-signs during disease pathogenesis. CONCLUSION: The proposed approach is valuable in early disease detection, monitoring efficacy of therapy, and clinical decision making.

Intensive Surveillance of Porcine-Rhesus Kidney Xenotransplant Using Different Ultrasound Techniques.

BACKGROUND: Significant progress has been made in kidney xenotransplantation in the past few years, and this field is accelerating towards clinical translation. Therefore, surveillance of the xenograft with appropriate tools is of great importance. Ultrasonography has been widely used in kidney allotransplantation and served as an economical and non-invasive method to monitor the allograft. However, questions remain whether the ultrasonographic criteria established for human kidney allograft could also be applied in xenotransplantation. METHODS: In the current study, we established a porcine-rhesus life sustaining kidney xenotransplantation model. The xenograft underwent intensive surveillance using gray-scale, colorful Doppler ultrasound as well as 2D shear wave elastography. The kidney growth, blood perfusion, and cortical stiffness were measured twice a day. These parameters were compared with the clinical data including urine output, chemistry, and pathological findings. RESULTS: The observation continued for 16 days after transplantation. Decline of urine output and elevated serum creatinine were observed on POD9 and biopsy proven antibody-mediated rejection was seen on the same day. The xenograft underwent substantial growth, with the long axis length increased by 32% and the volume increased by threefold at the end of observation. The resistive index of the xenograft arteries elevated in response to rejection, together with impaired cortical perfusion, while the peak systolic velocity (PSV) was not compromised. The cortical stiffness also increased along with rejection. CONCLUSION: In summary, the ultrasound findings of kidney xenograft shared similarities with those in allograft but possessed some unique features. A modified criteria needs to be established for further application of ultrasound in kidney xenotransplantation.

Amylase and lipase levels in the metabolic syndrome and type 2 diabetes: A longitudinal study in rhesus monkeys.

Latent associations between low serum amylase and reduced plasma insulin levels and increased adiposity have been described previously in a small study of asymptomatic middle-aged humans. In the present study, we sought to determine the nature of such changes during the longitudinal progression from metabolically normal to overt type 2 diabetes mellitus (T2DM) in nonhuman primates (NHPs), a disease that appears to be the same in both pathophysiology and underlying mechanisms as that which most commonly develops in middle-aged adult humans. Amylase and lipase levels were characterized in 157 unrelated adult rhesus monkeys (Macaca mulatta); 38% developed T2DM while under study. In all monkeys, multivariable linear regression analysis revealed that amylase could be negatively predicted by % body fat (ß -0.29; p = 0.002), age (ß -0.27; p = 0.005), and HbA1c (ß -0.18; p = 0.037). Amylase levels were positively predicted by lipase levels (ß = 0.19; p = -0.024) in all NHPs included in the study. Amylase was significantly lower in NHPs with metabolic syndrome (p < 0.001), prediabetes (PreDM) (p < 0.001), and T2DM (p < 0.001) compared to metabolically normal adult NHPs. Lipase increased in NHPs with PreDM (p = 0.005) and T2DM (p = 0.04) compared to normal NHPs. This is the first longitudinal study of any species, including humans, to show the dynamics of amylase and lipase during the metabolic progression from normal to metabolic syndrome, to PreDM and then to overt T2DM. The extraordinary similarity between humans and monkeys in T2DM, in pancreatic pathophysiology and in metabolic functions give these findings high translational value.

Metagenomic comparison of intestinal microbiota between normal and liver fibrotic rhesus macaques (Macaca mulatta).

Liver fibrosis is an important pathological process in chronic liver disease and cirrhosis. Recent studies have found a close association between intestinal microbiota and the development of liver fibrosis. To determine whether there are differences in the intestinal microbiota between rhesus macaques with liver fibrosis (MG) and normal rhesus macaques (MN), fecal samples were collected from 8 male MG and 12 male MN. The biological composition of the intestinal microbiota was then detected using 16S rRNA gene sequencing. The results revealed statistically significant differences in ASVs and Chao1 in the alpha-diversity and the beta-diversity of intestinal microbiota between MG and MN. Both groups shared Prevotella and Lactobacillus as common dominant microbiota. However, beneficial bacteria such as Lactobacillus were significantly less abundant in MG (P = 0.02). Predictive functional analysis using PICRUSt2 gene prediction revealed that MG exhibited a higher relative abundance of functions related to substance transport and metabolic pathways. This study may provide insight into further exploration of the mechanisms by which intestinal microbiota affect liver fibrosis and its potential future use in treating liver fibrosis.

Multiscale chemogenetic dissection of fronto-temporal top-down regulation for object memory in primates.

Visual object memory is a fundamental element of various cognitive abilities, and the underlying neural mechanisms have been extensively examined especially in the anterior temporal cortex of primates. However, both macroscopic large-scale functional network in which this region is embedded and microscopic neuron-level dynamics of top-down regulation it receives for object memory remains elusive. Here, we identified the orbitofrontal node as a critical partner of the anterior temporal node for object memory by combining whole-brain functional imaging during rest and a short-term object memory task in male macaques. Focal chemogenetic silencing of the identified orbitofrontal node downregulated both the local orbitofrontal and remote anterior temporal nodes during the task, in association with deteriorated mnemonic, but not perceptual, performance. Furthermore, imaging-guided neuronal recordings in the same monkeys during the same task causally revealed that orbitofrontal top-down modulation enhanced stimulus-selective mnemonic signal in individual anterior temporal neurons while leaving bottom-up perceptual signal unchanged. Furthermore, similar activity difference was also observed between correct and mnemonic error trials before silencing, suggesting its behavioral relevance. These multifaceted but convergent results provide a multiscale causal understanding of dynamic top-down regulation of the anterior temporal cortex along the ventral fronto-temporal network underpinning short-term object memory in primates.

Comprehensive Flow Cytometric, Immunohistologic, and Molecular Assessment of Thymus Function in Rhesus Macaques.

The critical importance of the thymus for generating new naive T cells that protect against novel infections and are tolerant to self-antigens has led to a recent revival of interest in monitoring thymic function in species other than humans and mice. Nonhuman primates such as rhesus macaques (Macaca mulatta) provide particularly useful animal models for translational research in immunology. In this study, we tested the performance of a 15-marker multicolor Ab panel for flow cytometric phenotyping of lymphocyte subsets directly from rhesus whole blood, with validation by thymectomy and T cell depletion. Immunohistochemical and multiplex RNA expression analysis of thymus tissue biopsies and molecular assays on PBMCs were used to further validate thymus function. Results identify Ab panels that can accurately classify rhesus naive T cells (CD3+CD45RA+CD197+ or CD3+CD28+CD95-) and recent thymic emigrants (CD8+CD28+CD95-CD103+CD197+) using just 100 µl of whole blood and commercially available fluorescent Abs. An immunohistochemical panel reactive with pan-cytokeratin (CK), CK14, CD3, Ki-67, CCL21, and TdT provides histologic evidence of thymopoiesis from formalin-fixed, paraffin-embedded thymus tissues. Identification of mRNAs characteristic of both functioning thymic epithelial cells and developing thymocytes and/or molecular detection of products of TCR gene rearrangement provide additional complementary methods to evaluate thymopoiesis, without requiring specific Abs. Combinations of multiparameter flow cytometry, immunohistochemistry, multiplex gene expression, and TCR excision circle assays can comprehensively evaluate thymus function in rhesus macaques while requiring only minimal amounts of peripheral blood or biopsied thymus tissue.

Facial expressivity in dominant macaques is linked to group cohesion.

Social living affords primates (including humans) many benefits. Communication has been proposed to be the key mechanism used to bond social connections, which could explain why primates have evolved such expressive faces. We assessed whether the facial expressivity of the dominant male (quantified from the coding of anatomically based facial movement) was related to social network properties (based on social proximity and grooming) in nine groups of captive rhesus macaques (Macaca mulatta) housed in uniform physical and social environments. More facially expressive dominant male macaques were more socially connected and had more cohesive social groups. These findings show that inter-individual differences in facial expressivity are related to differential social outcomes at both an individual and group level. More expressive individuals occupy more beneficial social positions, which could help explain the selection for complex facial communication in primates.

Hippocampal connectivity patterns echo macroscale cortical evolution in the primate brain.

While the hippocampus is key for human cognitive abilities, it is also a phylogenetically old cortex and paradoxically considered evolutionarily preserved. Here, we introduce a comparative framework to quantify preservation and reconfiguration of hippocampal organisation in primate evolution, by analysing the hippocampus as an unfolded cortical surface that is geometrically matched across species. Our findings revealed an overall conservation of hippocampal macro- and micro-structure, which shows anterior-posterior and, perpendicularly, subfield-related organisational axes in both humans and macaques. However, while functional organisation in both species followed an anterior-posterior axis, we observed a marked reconfiguration in the latter across species, which mirrors a rudimentary integration of the default-mode-network in non-human primates. Here we show that microstructurally preserved regions like the hippocampus may still undergo functional reconfiguration in primate evolution, due to their embedding within heteromodal association networks.

Development of a hydrogen peroxide-inactivated vaccine that protects against viscerotropic yellow fever in a non-human primate model.

Yellow fever virus (YFV) is endemic in >40 countries and causes viscerotropic disease with up to 20%-60% mortality. Successful live-attenuated yellow fever (YF) vaccines were developed in the mid-1930s, but their use is restricted or formally contraindicated in vulnerable populations including infants, the elderly, and people with compromised immune systems. In these studies, we describe the development of a next-generation hydrogen peroxide-inactivated YF vaccine and determine immune correlates of protection based on log neutralizing index (LNI) and neutralizing titer-50% (NT50) studies. In addition, we compare neutralizing antibody responses and protective efficacy of hydrogen peroxide-inactivated YF vaccine candidates to live-attenuated YFV-17D (YF-VAX) in a rhesus macaque model of viscerotropic YF. Our results indicate that an optimized, inactivated YF vaccine elicits protective antibody responses that prevent viral dissemination and lethal infection in rhesus macaques and may be a suitable alternative for vaccinating vulnerable populations who are not eligible to receive replicating live-attenuated YF vaccines.

Genotypic variation in the promoter region of the CRH-248 gene interacts with early rearing experiences to disrupt the development of the HPA axis in infant rhesus macaques (Macaca mulatta).

Aberrant functioning of the hypothalamic-pituitary-adrenal (HPA) axis is a hallmark of conditions such as depression, anxiety disorders, and post-traumatic stress disorder. Early-life adversity and genetic variation can interaction to disrupt HPA axis regulation, potentially contributing to certain forms of psychopathology. This study employs a rhesus macaque model to investigate how early parental neglect interacts with a single nucleotide polymorphism within the promoter region of the corticotropin-releasing hormone (CRH-248) gene, impacting the development of the HPA axis. For the initial six months of life, 307 rhesus monkey infants (n = 146 females, n = 161 males) were either reared with their mothers (MR) in conditions emulating the natural environment (control group) or raised without maternal care in groups with constant or 3-hours daily access to same-aged peers (NR). Blood samples collected on days 30, 60, 90, and 120 of life under stressful conditions were assayed for plasma cortisol and adrenocorticotropic hormone (ACTH) concentrations. Findings revealed that NR subjects exhibited a significant blunting of both ACTH and cortisol concentrations. Notably, there was a gene-by-environment interaction observed for ACTH and cortisol levels, with NR subjects with the polymorphism displaying higher ACTH concentrations and lower cortisol concentrations. To the extent that these results generalize to humans, they suggest that early parental neglect may render individuals vulnerable to HPA axis dysfunction, a susceptibility that is modulated by CRH-248 genotype-a gene-by-environment interaction that leaves a lasting developmental signature.

Correlation between the number of interstitial neurons of the white matter and number of neurons within cortical layers: Histological analyses in postnatal macaque.

We have examined the number and distribution of NeuN-immunoreactive cortical white matter interstitial cells (WMICs) and compared them to the neurons in layers 1-6 across the overlying cortex in coronal sections from postnatal macaques. The data have been gathered from over 300 selected regions at gyral crowns, at sulci, and at linear regions of the cortex where we also determined cortical layer thicknesses: standard thicknesses and tangential thicknesses. Cortical thicknesses and cell numbers showed variability according to gyral, linear, or sulcal regions. In spite of these variations, our standardized cell numbers in layers 1 to 6b and interstitial cells underlying layer 6b-white matter boundary have shown a consistent correlation between the number of WMICs and the number of layer 5 and 6a cortical neurons on all cortical regions studied: for each WMIC, there are on the order of five cortical neurons in layer 5 and approximately three cortical neurons in layer 6a, irrespective of the origins of the selected cortical area or whether they are from gyral, linear, or sulcal regions. We propose that the number of interstitial neurons in the postnatal macaque cortex is correlated to the density of neurons within layers 5 and 6a and, from a clinical perspective, the change in density or distribution of interstitial neurons in schizophrenia or epilepsy may in fact be linked to the number of layers 5 and 6a neurons.

Lingual neoplasia in nonhuman primates: Description of five cases and a literature review.

BACKGROUND: Documentation of lingual tumors is scarce in nonhuman primates. METHODS: Through a multi-institutional retrospective study we compile cases of primary and metastatic neoplasia in non-human primates. RESULTS: We describe five cases of lingual neoplasia. Three cases are primary lingual tumors: chondro-osteoblastic lipoma in a howler monkey, squamous cell carcinoma, and fibroma in two baboons. We describe two cases of metastatic lymphoma in the tongue in rhesus macaques. A literature review of published lingual neoplasia in nonhuman primates is included in this manuscript. CONCLUSION: Lingual neoplasia is seldom reported in non-human primates.

Recombinant OC43 SARS-CoV-2 spike replacement virus: An improved BSL-2 proxy virus for SARS-CoV-2 neutralization assays.

We generated a replication-competent OC43 human seasonal coronavirus (CoV) expressing the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike in place of the native spike (rOC43-CoV2 S). This virus is highly attenuated relative to OC43 and SARS-CoV-2 in cultured cells and animals and is classified as a biosafety level 2 (BSL-2) agent by the NIH biosafety committee. Neutralization of rOC43-CoV2 S and SARS-CoV-2 by S-specific monoclonal antibodies and human sera is highly correlated, unlike recombinant vesicular stomatitis virus-CoV2 S. Single-dose immunization with rOC43-CoV2 S generates high levels of neutralizing antibodies against SARS-CoV-2 and fully protects human ACE2 transgenic mice from SARS-CoV-2 lethal challenge, despite nondetectable replication in respiratory and nonrespiratory organs. rOC43-CoV2 S induces S-specific serum and airway mucosal immunoglobulin A and IgG responses in rhesus macaques. rOC43-CoV2 S has enormous value as a BSL-2 agent to measure S-specific antibodies in the context of a bona fide CoV and is a candidate live attenuated SARS-CoV-2 mucosal vaccine that preferentially replicates in the upper airway.

Context-invariant beliefs are supported by dynamic reconfiguration of single unit functional connectivity in prefrontal cortex of male macaques.

Natural behaviors occur in closed action-perception loops and are supported by dynamic and flexible beliefs abstracted away from our immediate sensory milieu. How this real-world flexibility is instantiated in neural circuits remains unknown. Here, we have male macaques navigate in a virtual environment by primarily leveraging sensory (optic flow) signals, or by more heavily relying on acquired internal models. We record single-unit spiking activity simultaneously from the dorsomedial superior temporal area (MSTd), parietal area 7a, and the dorso-lateral prefrontal cortex (dlPFC). Results show that while animals were able to maintain adaptive task-relevant beliefs regardless of sensory context, the fine-grain statistical dependencies between neurons, particularly in 7a and dlPFC, dynamically remapped with the changing computational demands. In dlPFC, but not 7a, destroying these statistical dependencies abolished the area's ability for cross-context decoding. Lastly, correlational analyses suggested that the more unit-to-unit couplings remapped in dlPFC, and the less they did so in MSTd, the less were population codes and behavior impacted by the loss of sensory evidence. We conclude that dynamic functional connectivity between neurons in prefrontal cortex maintain a stable population code and context-invariant beliefs during naturalistic behavior.

Orthogonality of sensory and contextual categorical dynamics embedded in a continuum of responses from the second somatosensory cortex.

How does the brain simultaneously process signals that bring complementary information, like raw sensory signals and their transformed counterparts, without any disruptive interference? Contemporary research underscores the brain's adeptness in using decorrelated responses to reduce such interference. Both neurophysiological findings and artificial neural networks support the notion of orthogonal representation for signal differentiation and parallel processing. Yet, where, and how raw sensory signals are transformed into more abstract representations remains unclear. Using a temporal pattern discrimination task in trained monkeys, we revealed that the second somatosensory cortex (S2) efficiently segregates faithful and transformed neural responses into orthogonal subspaces. Importantly, S2 population encoding for transformed signals, but not for faithful ones, disappeared during a nondemanding version of this task, which suggests that signal transformation and their decoding from downstream areas are only active on-demand. A mechanistic computation model points to gain modulation as a possible biological mechanism for the observed context-dependent computation. Furthermore, individual neural activities that underlie the orthogonal population representations exhibited a continuum of responses, with no well-determined clusters. These findings advocate that the brain, while employing a continuum of heterogeneous neural responses, splits population signals into orthogonal subspaces in a context-dependent fashion to enhance robustness, performance, and improve coding efficiency.

Myelin dystrophy impairs signal transmission and working memory in a multiscale model of the aging prefrontal cortex.

Normal aging leads to myelin alterations in the rhesus monkey dorsolateral prefrontal cortex (dlPFC), which are positively correlated with degree of cognitive impairment. It is hypothesized that remyelination with shorter and thinner myelin sheaths partially compensates for myelin degradation, but computational modeling has not yet explored these two phenomena together systematically. Here, we used a two-pronged modeling approach to determine how age-related myelin changes affect a core cognitive function: spatial working memory. First, we built a multicompartment pyramidal neuron model fit to monkey dlPFC empirical data, with an axon including myelinated segments having paranodes, juxtaparanodes, internodes, and tight junctions. This model was used to quantify conduction velocity (CV) changes and action potential (AP) failures after demyelination and subsequent remyelination. Next, we incorporated the single neuron results into a spiking neural network model of working memory. While complete remyelination nearly recovered axonal transmission and network function to unperturbed levels, our models predict that biologically plausible levels of myelin dystrophy, if uncompensated by other factors, can account for substantial working memory impairment with aging. The present computational study unites empirical data from ultrastructure up to behavior during normal aging, and has broader implications for many demyelinating conditions, such as multiple sclerosis or schizophrenia.

Recapitulation of HIV-1 Neutralization Breadth in Plasma by the Combination of Two Broadly Neutralizing Antibodies from Different Lineages in the Same SHIV-Infected Rhesus Macaque.

Viral infection generally induces polyclonal neutralizing antibody responses. However, how many lineages of antibody responses can fully represent the neutralization activities in sera has not been well studied. Using the newly designed stable HIV-1 Env trimer as hook, we isolated two distinct broadly neutralizing antibodies (bnAbs) from Chinese rhesus macaques infected with SHIV1157ipd3N4 for 5 years. One lineage of neutralizing antibodies (JT15 and JT16) targeted the V2-apex in the Env trimers, similar to the J038 lineage bnAbs identified in our previous study. The other lineage neutralizing antibody (JT18) targeted the V3 crown region in the Env, which strongly competed with human 447-52D. Each lineage antibody neutralized a different set of viruses. Interestingly, when the two neutralizing antibodies from different lineages isolated from the same macaque were combined, the mixture had a neutralization breath very similar to that from the cognate sera. Our study demonstrated that a minimum of two different neutralizing antibodies can fully recapitulate the serum neutralization breadth. This observation can have important implications in AIDS vaccine design.