Comparative transcriptome analysis between rhesus macaques ( Macaca mulatta) and crab-eating macaques ( M. fascicularis)

Understanding gene expression variations between species is pivotal for deciphering the evolutionary diversity in phenotypes. Rhesus macaques ( Macaca mulatta, MMU) and crab-eating macaques ( M. fascicularis, MFA) serve as crucial nonhuman primate biomedical models with different phenotypes. To date, however, large-scale comparative transcriptome research between these two species has not yet been fully explored. Here, we conducted systematic comparisons utilizing newly sequenced RNA-seq data from 84 samples (41 MFA samples and 43 MMU samples) encompassing 14 common tissues. Our findings revealed a small fraction of genes (3.7%) with differential expression between the two species, as well as 36.5% of genes with tissue-specific expression in both macaques. Comparison of gene expression between macaques and humans indicated that 22.6% of orthologous genes displayed differential expression in at least two tissues. Moreover, 19.41% of genes that overlapped with macaque-specific structural variants showed differential expression between humans and macaques. Of these, the FAM220A gene exhibited elevated expression in humans compared to macaques due to lineage-specific duplication. In summary, this study presents a large-scale transcriptomic comparison between MMU and MFA and between macaques and humans. The discovery of gene expression variations not only enhances the biomedical utility of macaque models but also contributes to the wider field of primate genomics.

Full-length 16S rDNA sequencing based on Oxford Nanopore Technologies revealed the association between gut-pharyngeal microbiota and tuberculosis in cynomolgus macaques

Tuberculosis (TB) is an infectious disease caused by the Mycobacterium tuberculosis complex (Mtbc), which develops from asymptomatic latent TB to active stages. The microbiome was purposed as a potential factor affecting TB pathogenesis, but the study was limited. The present study explored the association between gut-pharyngeal microbiome and TB stages in cynomolgus macaques using the full-length 16S rDNA amplicon sequencing based on Oxford Nanopore Technologies. The total of 71 macaques was divided into TB (-) control, TB (+) latent and TB (+) active groups. The differential abundance analysis showed that Haemophilus hemolyticus was decreased, while Prevotella species were increased in the pharyngeal microbiome of TB (+) macaques. In addition, Eubacterium coprostanoligenes in the gut was enriched in TB (+) macaques. Alteration of these bacteria might affect immune regulation and TB severity, but details of mechanisms should be further explored and validated. In summary, microbiota may be associated with host immune regulation and affect TB progression. The findings suggested the potential mechanisms of host-microbes interaction, which may improve the understanding of the role of microbiota and help develop therapeutics for TB in the future.

The multilevel extensive diversity across the cynomolgus macaque captured by ultra-deep adaptive immune receptor repertoire sequencing

The extent to which AIRRs differ among and within individuals remains elusive. Via ultra-deep repertoire sequencing of 22 and 25 tissues in three cynomolgus macaques, respectively, we identified 84 and 114 novel IGHV and TRBV alleles, confirming 72 (85.71%) and 100 (87.72%) of them. The heterogeneous V gene usage patterns were influenced, in turn, by genetics, isotype (for BCRs only), tissue group, and tissue. A higher proportion of intragroup shared clones in the intestinal tissues than those in other tissues suggests a close intra-intestinal adaptive immunity network. Significantly higher mutation burdens in the public clones and the inter-tissue shared IgM and IgD clones indicate that they might target the shared antigens. This study reveals the extensive heterogeneity of the AIRRs at various levels and has broad fundamental and clinical implications. The data generated here will serve as an invaluable resource for future studies on adaptive immunity in health and diseases.

The HER2-directed antibody-drug conjugate DHES0815A in advanced and/or metastatic breast cancer: preclinical characterization and phase 1 trial results

Approved antibody-drug conjugates (ADCs) for HER2-positive breast cancer include trastuzumab emtansine and trastuzumab deruxtecan. To develop a differentiated HER2 ADC, we chose an antibody that does not compete with trastuzumab or pertuzumab for binding, conjugated to a reduced potency PBD (pyrrolobenzodiazepine) dimer payload. PBDs are potent cytotoxic agents that alkylate and cross-link DNA. In our study, the PBD dimer is modified to alkylate, but not cross-link DNA. This HER2 ADC, DHES0815A, demonstrates in vivo efficacy in models of HER2-positive and HER2-low cancers and is well-tolerated in cynomolgus monkey safety studies. Mechanisms of action include induction of DNA damage and apoptosis, activity in non-dividing cells, and bystander activity. A dose-escalation study (ClinicalTrials.gov: NCT03451162) in patients with HER2-positive metastatic breast cancer, with the primary objective of evaluating the safety and tolerability of DHES0815A and secondary objectives of characterizing the pharmacokinetics, objective response rate, duration of response, and formation of anti-DHES0815A antibodies, is reported herein. Despite early signs of anti-tumor activity, patients at higher doses develop persistent, non-resolvable dermal, ocular, and pulmonary toxicities, which led to early termination of the phase 1 trial.

Adaptation of HLA testing to characterize the cynomolgus macaque MHC polymorphisms and alloantibody signatures

Nonhuman primates are the closest animal models to humans with respect to genetics and physiology. Consequently, a critical component of immunogenetics research relies on drawing inferences from the cynomolgus macaque to inform human trials. Despite the conserved organization of the Major Histocompatibility Complex (MHC) between cynomolgus macaques and humans, MHC genotyping of cynomolgus macaques is challenging due to high rates of copy number variants, duplications, and rearrangements, particularly at the MHC class I loci. Furthermore, the limited availability of commercial reagents specific to cynomolgus macaques that can be used to characterize anti-MHC class I and class II antibody (Ab) specificities in cynomolgus macaques presents a major bottleneck in translational research. Here we successfully characterized cynomolgus macaque Mafa class I and class II serologic specificities in 86 animals originating from various geographical regions using the complement dependent cytotoxicity (CDC) assay with human HLA class I and class II monoclonal antibody (mAb) typing trays. Further, we successfully induced and characterized anti-Mafa class I and class II alloantibody specificity using HLA single antigen bead assays. We also subsequently tracked the alloAb burden in the animals during treatment with anti-B lymphocyte stimulator (BLyS) treatment. Altogether, these methods can be easily used in translational research to serotype MHC class I and class II specificity in macaques, characterize their alloAb specificity, and evaluate the efficacy of novel therapeutic modalities in depleting circulating alloAbs in these animals.

MHC-DRB alleles with amino acids Val11, Phe13, and the shared epitopes promote collagen-induced arthritis and a rapid IgG1 response in Filipino cynomolgus macaques

Macaques are useful animal models for studying the pathogenesis of rheumatoid arthritis (RA) and the development of anti-rheumatic drugs. The purpose of this study was to identify the major histocompatibility complex (MHC) polymorphisms associated with the pathology of collagen-induced arthritis (CIA) and anti-collagen IgG induction in a cynomolgus macaque model, as MHC polymorphisms affect the onset of CIA in other animal models. Nine female Filipino cynomolgus macaques were immunized with bovine type II collagen (b-CII) to induce CIA, which was diagnosed clinically by scoring the symptoms of joint swelling over 9 weeks. MHC polymorphisms and anti-b-CII antibody titers were compared between symptomatic and asymptomatic macaques. Four of 9 (44%) macaques were defined as the CIA-affected group. Anti-b-CII IgG in the affected group increased in titer approximately 3 weeks earlier compared with the asymptomatic group. The mean plasma IgG1 titer in the CIA-affected group was significantly higher (p < 0.05) than that of the asymptomatic group. Furthermore, the cynomolgus macaque MHC (Mafa)-DRB1*10:05 or Mafa-DRB1*10:07 alleles, which contain the well-documented RA-susceptibility five amino acid sequence known as the shared epitope (SE) in positions 70 to 74, with valine at position 11 (Val11, V11) and phenylalanine at position 13 (Phe13, F13), were detected in the affected group. In contrast, no MHC polymorphisms specific to the asymptomatic group were identified. In conclusion, the presence of V11 and F13 along with SE in the MHC-DRB1 alleles seems essential for the production of IgG1 and the rapid induction of severe CIA in female Filipino cynomolgus macaques.

Characterizing bacterial and fungal communities along the longitudinal axis of the intestine in cynomolgus monkeys

IMPORTANCE

Gut microbiota varies along the gastrointestinal (GI) tract and exerts profound influences on the host's physiology, immunity, and nutrition. Given that gut microbes interact with the host closely and the gastrointestinal function differed from the small to the large intestine, it is essential to characterize the gut biogeography of the microbial community. Here, we focused on intestinal bacteria and fungi in cynomolgus monkeys and determined their spatial distribution along the GI tract by performing 16S and 18S rRNA gene sequencing. The composition and function of bacterial and fungal communities differed significantly at different biogeographic sites of the intestine, and the site-specific correlations between intestinal bacteria and fungi were revealed. Thus, our studies characterized the gut biogeography of bacteria and fungi in NHPs and revealed their site-specific correlations along the GI tract.

Liver injury in cynomolgus monkeys following intravenous and intrathecal scAAV9 gene therapy delivery

Hepatotoxicity associated with intravenous/intrathecal adeno-associated virus (AAV) gene therapy has been observed in preclinical species and patients. In nonhuman primates, hepatotoxicity following self-complementary AAV9 administration varies from asymptomatic transaminase elevation with minimal to mild microscopic changes to symptomatic elevations of liver function and thromboinflammatory markers with microscopic changes consistent with marked hepatocellular necrosis and deteriorating clinical condition. These transient acute liver injury marker elevations occur from 3-4 days post intravenous administration to ∼2 weeks post intrathecal administration. No transaminase elevation or microscopic changes were observed with intrathecal administration of empty capsids or a "promoterless genome" vector, suggesting that liver injury after cerebrospinal fluid dosing in nonhuman primates is driven by viral transduction and transgene expression. Co-administration of prednisolone after intravenous or intrathecal dosing did not prevent liver enzyme or microscopic changes despite a reduction of T lymphocyte infiltration in liver tissue. Similarly, co-administration of rituximab/everolimus with intrathecal dosing failed to block AAV-driven hepatotoxicity. Self-complementary AAV-induced acute liver injury appears to correlate with high hepatocellular vector load, macrophage activation, and type 1 interferon innate virus-sensing pathway responses. The current work characterizes key aspects pertaining to early AAV-driven hepatotoxicity in cynomolgus macaques, highlighting the usefulness of this nonclinical species in that context.

Comparative analysis of gut microbiota between common (Macaca fascicularis fascicularis) and Burmese (M. f. aurea) long-tailed macaques in different habitats

The environment has an important effect on the gut microbiota-an essential part of the host's health-and is strongly influenced by the dietary pattern of the host as these together shape the composition and functionality of the gut microbiota in humans and other animals. This study compared the gut microbiota of Macaca fascicularis fascicularis and M. f. aurea in mangrove and island populations using 16S rRNA gene sequencing on a nanopore platform to investigate the effect of the environment and/or diet. The results revealed that the M. f. fascicularis populations that received anthropogenic food exhibited a higher richness and evenness of gut microbiota than the M. f. aurea populations in different habitats. Firmicutes and Bacteroidetes were the two most abundant bacterial phyla in the gut microbiota of both these subspecies; however, the relative abundance of these phyla was significantly higher in M. f. aurea than in M. f. fascicularis. This variation in the gut microbiota between the two subspecies in different habitats mostly resulted from the differences in their diets. Moreover, the specific adaptation of M. f. aurea to different environments with a different food availability had a significant effect on their microbial composition.

C9orf72 poly(PR) aggregation in nucleus induces ALS/FTD-related neurodegeneration in cynomolgus monkeys

Poly(PR) is a dipeptide repeat protein comprising proline and arginine residues. It is one of the translational product of expanded G4C2 repeats in the C9orf72 gene, and its accumulation is contributing to the neuropathogenesis of C9orf72-associated amyotrophic lateral sclerosis and/or frontotemporal dementia (C9-ALS/FTD). In this study, we demonstrate that poly(PR) protein alone is sufficient to induce neurodegeneration related to ALS/FTD in cynomolgus monkeys. By delivering poly(PR) via AAV, we observed that the PR proteins were located within the nucleus of infected cells. The expression of (PR)50 protein, consisting of 50 PR repeats, led to increased loss of cortical neurons, cytoplasmic lipofuscin, and gliosis in the brain, as well as demyelination and loss of ChAT positive neurons in the spinal cord of monkeys. While, these pathologies were not observed in monkeys expressing (PR)5, a protein comprising only 5 PR repeats. Furthermore, the (PR)50-expressing monkeys exhibited progressive motor deficits, cognitive impairment, muscle atrophy, and abnormal electromyography (EMG) potentials, which closely resemble clinical symptoms seen in C9-ALS/FTD patients. By longitudinally tracking these monkeys, we found that changes in cystatin C and chitinase-1 (CHIT1) levels in the cerebrospinal fluid (CSF) corresponded to the phenotypic progression of (PR)50-induced disease. Proteomic analysis revealed that the major clusters of dysregulated proteins were nuclear-localized, and downregulation of the MECP2 protein was implicated in the toxic process of poly(PR). This research indicates that poly(PR) expression alone induces neurodegeneration and core phenotypes associated with C9-ALS/FTD in monkeys, which may provide insights into the mechanisms of disease pathogenesis.

Genome sequencing and application of Taiwanese macaque Macaca cyclopis

Formosan macaque (Macaca cyclopis) is the only non-human primate in Taiwan Island. We performed de novo hybrid assembly for M. cyclopis using Illumina paired-end short reads, mate-pair reads and Nanopore long reads and obtained 5065 contigs with a N50 of 2.66 megabases. M. cyclopis contigs >  = 10 kb were assigned to chromosomes using Indian rhesus macaque (Macaca mulatta mulatta) genome assembly Mmul_10 as reference, resulting in a draft of M. cyclopis genome of 2,846,042,475 bases, distributed in 21 chromosomes. The draft genome contains 23,462 transcriptional origins (genes), capable of expressing 716,231 exons in 59,484 transcripts. Genome-based phylogenetic study using the assembled M. cyclopis genome together with genomes of four other macaque species, human, orangutan and chimpanzee showed similar result as previously reported. However, the M. cyclopis species was found to diverge from Chinese M. mulatta lasiota about 1.8 million years ago. Fossil gene analysis detected the presence of gap and pol endogenous viral elements of simian retrovirus in all macaques tested, including M. fascicularis, M. m. mulatta and M. cyclopis. However, M. cyclopis showed ~ 2 times less in number and more uniform in chromosomal locations. The constrain in foreign genome disturbance, presumably due to geographical isolation, should be able to simplify genomics-related investigations, making M. cyclopis an ideal primate species for medical research.

Genotyping of Macaque Population on the Brain-Derived Neurotrophic Factor Gene Polymorphism by Mismatch Amplification Mutation Assay (MAMA)-PCR

Human rs6265 (196G>A) polymorphism in the BDNF gene is associated with many clinically significant phenotypic manifestations. Rhesus monkey (Macaca mulatta) has a functionally significant rs309950446 ( 136G>A) polymorphism. To determine this polymorphism in macaques, we used mismatch amplification mutation assay (MAMA)-PCR method with non-complementary nucleotide to the template chain at the 3rd position from the 3'-end of the allele-specific primers (mismatch primers), which allowed the best discrimination of the alleles. Genotyping of male rhesus monkeys (n=178) and cynomolgus monkeys (Macaca fascicularis) (n=90) was carried out. The A/A, G/G, and G/A genotypes were found in 16, 34, and 50% rhesus macaques, respectively. In the cynomolgus macaques, the mutant polymorphic allele was not detected. The study results allow considering rhesus macaques as a potential biological model for assessment of the gen-environment interaction of the BDNF gene polymorphism.

Brain-wide and cell-specific transcriptomic insights into MRI-derived cortical morphology in macaque monkeys

Integrative analyses of transcriptomic and neuroimaging data have generated a wealth of information about biological pathways underlying regional variability in imaging-derived brain phenotypes in humans, but rarely in nonhuman primates due to the lack of a comprehensive anatomically-defined atlas of brain transcriptomics. Here we generate complementary bulk RNA-sequencing dataset of 819 samples from 110 brain regions and single-nucleus RNA-sequencing dataset, and neuroimaging data from 162 cynomolgus macaques, to examine the link between brain-wide gene expression and regional variation in morphometry. We not only observe global/regional expression profiles of macaque brain comparable to human but unravel a dorsolateral-ventromedial gradient of gene assemblies within the primate frontal lobe. Furthermore, we identify a set of 971 protein-coding and 34 non-coding genes consistently associated with cortical thickness, specially enriched for neurons and oligodendrocytes. These data provide a unique resource to investigate nonhuman primate models of human diseases and probe cross-species evolutionary mechanisms.

Complete sequencing of a cynomolgus macaque major histocompatibility complex haplotype

Macaques provide the most widely used nonhuman primate models for studying the immunology and pathogenesis of human diseases. Although the macaque major histocompatibility complex (MHC) region shares most features with the human leukocyte antigen (HLA) region, macaques have an expanded repertoire of MHC class I genes. Although a chimera of two rhesus macaque MHC haplotypes was first published in 2004, the structural diversity of MHC genomic organization in macaques remains poorly understood owing to a lack of adequate genomic reference sequences. We used ultralong Oxford Nanopore and high-accuracy Pacific Biosciences (PacBio) HiFi sequences to fully assemble the ∼5.2-Mb M3 haplotype of an MHC-homozygous, Mauritian-origin cynomolgus macaque (Macaca fascicularis). The MHC homozygosity allowed us to assemble a single MHC haplotype unambiguously and avoid chimeric assemblies that hampered previous efforts to characterize this exceptionally complex genomic region in macaques. The high quality of this new assembly is exemplified by the identification of an extended cluster of six Mafa-AG genes that contains a recent duplication with a highly similar ∼48.5-kb block of sequence. The MHC class II region of this M3 haplotype is similar to the previously sequenced rhesus macaque haplotype and HLA class II haplotypes. The MHC class I region, in contrast, contains 13 MHC-B genes, four MHC-A genes, and three MHC-E genes (vs. 19 MHC-B, two MHC-A, and one MHC-E in the previously sequenced haplotype). These results provide an unambiguously assembled single contiguous cynomolgus macaque MHC haplotype with fully curated gene annotations that will inform infectious disease and transplantation research.

Selection of a PD-1 blocking antibody from a novel fully human phage display library

Programmed cell death protein 1 (PD-1) is an immunoregulatory target which is recognized by different monoclonal antibodies, approved for the therapy of multiple types of cancer. Different anti-PD-1 antibodies display different therapeutic properties and there is a pharmaceutical interest to generate and characterize novel anti-PD-1 antibodies. We screened multiple human antibody phage display libraries to target novel epitopes on the PD-1 surface and we discovered a unique and previously undescribed binding specificity (termed D12) from a new antibody library (termed AMG). The library featured antibody fragments in single-chain fragment variable (scFv) format, based on the IGHV3-23*03 (VH ) and IGKV1-39*01 (Vκ) genes. The D12 antibody was characterized by surface plasmon resonance (SPR), cross-reacted with the Cynomolgus monkey antigen and bound to primary human T cells, as shown by flow cytometry. The antibody blocked the PD-1/PD-L1 interaction in vitro with an EC50 value which was comparable to the one of nivolumab, a clinically approved antibody. The fine details of the interaction between D12 and PD-1 were elucidated by x-ray crystallography of the complex at a 3.5 Å resolution, revealing an unprecedented conformational change at the N-terminus of PD-1 following D12 binding, as well as partial overlap with the binding site for the cognate PD-L1 and PD-L2 ligands which prevents their binding. The results of the study suggest that the expansion of antibody library repertoires may facilitate the discovery of novel binding specificities with unique properties that hold promises for the modulation of PD-1 activity in vitro and in vivo.

Dynamic evolution of Mhc haplotypes in cynomolgus macaques of different geographic origins

The major histocompatibility complex (MHC) plays a key role in immune defense, and the Mhc genes of cynomolgus macaque display a high degree of polymorphism. Based on their geographic distribution, different populations of cynomolgus macaques are recognized. Here we present the characterization of the Mhc class I and II repertoire of a large pedigreed group of cynomolgus macaques originating from the mainland north of the isthmus of Kra (N = 42). Segregation analyses resulted in the definition of 81 unreported Mafa-A/B/DRB/DQ/DP haplotypes, which include 32 previously unknown DRB regions. In addition, we report 13 newly defined Mafa-A/B/DRB/DQ/DP haplotypes in a group of cynomolgus macaques originating from the mainland south of the isthmus of Kra/Maritime Southeast Asia (N = 16). A relatively high level of sharing of Mafa-A (51%) and Mafa-B (40%) lineage groups is observed between the populations native to the north and the south of isthmus of Kra. At the allelic level, however, the Mafa-A/B haplotypes seem to be characteristic of a population. An overall comparison of all currently known data revealed that each geographic population has its own specific combinations of Mhc class I and II haplotypes. This illustrates the dynamic evolution of the cynomolgus macaque Mhc region, which was most likely generated by recombination and maintained by selection due to the differential pathogenic pressures encountered in different geographic areas.

PINK1 gene mutation by pair truncated sgRNA/Cas9-D10A in cynomolgus monkeys

Mutations of PTEN-induced kinase I (PINK1) cause early-onset Parkinson's disease (PD) with selective neurodegeneration in humans. However, current PINK1 knockout mouse and pig models are unable to recapitulate the typical neurodegenerative phenotypes observed in PD patients. This suggests that generating PINK1 disease models in non-human primates (NHPs) that are close to humans is essential to investigate the unique function of PINK1 in primate brains. Paired single guide RNA (sgRNA)/Cas9-D10A nickases and truncated sgRNA/Cas9, both of which can reduce off-target effects without compromising on-target editing, are two optimized strategies in the CRISPR/Cas9 system for establishing disease animal models. Here, we combined the two strategies and injected Cas9-D10A mRNA and two truncated sgRNAs into one-cell-stage cynomolgus zygotes to target the PINK1 gene. We achieved precise and efficient gene editing of the target site in three newborn cynomolgus monkeys. The frame shift mutations of PINK1 in mutant fibroblasts led to a reduction in mRNA. However, western blotting and immunofluorescence staining confirmed the PINK1 protein levels were comparable to that in wild-type fibroblasts. We further reprogramed mutant fibroblasts into induced pluripotent stem cells (iPSCs), which showed similar ability to differentiate into dopamine (DA) neurons. Taken together, our results showed that co-injection of Cas9-D10A nickase mRNA and sgRNA into one-cell-stage cynomolgus embryos enabled the generation of human disease models in NHPs and target editing by pair truncated sgRNA/Cas9-D10A in PINK1 gene exon 2 did not impact protein expression.

Chromosomal-scale de novo genome assemblies of Cynomolgus Macaque and Common Marmoset

Cynomolgus macaque (Macaca fascicularis) and common marmoset (Callithrix jacchus) have been widely used in human biomedical research. Long-standing primate genome assemblies used the human genome as a reference for ordering and orienting the assembled fragments into chromosomes. Here we performed de novo genome assembly of these two species without any human genome-based bias observed in the genome assemblies released earlier. We assembled PacBio long reads, and the resultant contigs were scaffolded with Hi-C data, which were further refined based on Hi-C contact maps and alternate de novo assemblies. The assemblies achieved scaffold N50 lengths of 149 Mb and 137 Mb for cynomolgus macaque and common marmoset, respectively. The high fidelity of our assembly is also ascertained by BAC-end concordance in common marmoset. Our assembly of cynomolgus macaque outperformed all the available assemblies of this species in terms of contiguity. The chromosome-scale genome assemblies produced in this study are valuable resources for non-human primate models and provide an important baseline in human biomedical research.

Testing Sex-Biased Admixture Origin of Macaque Species Using Autosomal and X-Chromosomal Genomic Sequences

The role of sex-specific demography in hybridization and admixture of genetically diverged species and populations is essential to understand the origins of the genomic diversity of sexually reproducing organisms. In order to infer how sex-linked loci have been differentiated undergoing frequent hybridization and admixture, we examined 17 whole-genome sequences of seven species representing the genus Macaca, which shows frequent inter-specific hybridization and predominantly female philopatry. We found that hybridization and admixture were prevalent within these species. For three cases of suggested hybrid origin of species/subspecies, Macaca arctoides, Macaca fascicularis ssp. aurea, and Chinese Macaca mulatta, we examined the level of admixture of X chromosomes, which is less affected by male-biased migration than that of autosomes. In one case, we found that Macaca cyclopis and Macaca fuscata was genetically closer to Chinese M. mulatta than to the Indian M. mulatta, and the admixture level of Chinese M. mulatta and M. fuscata/cyclopis was more pronounced on the X chromosome than on autosomes. Since the mitochondrial genomes of Chinese M. mulatta, M. cyclopis, and M. fuscata were found to cluster together, and the mitochondrial genome of Indian M. mulatta is more distantly related, the observed pattern of genetic differentiation on X-chromosomal loci is consistent with the nuclear swamping hypothesis, in which strong, continuous male-biased introgression from the ancestral Chinese M. mulatta population to a population related to M. fuscata and M. cyclopis generated incongruencies between the genealogies of the mitochondrial and nuclear genomes.

Monkeys mutant for PKD1 recapitulate human autosomal dominant polycystic kidney disease

Autosomal dominant polycystic kidney disease (ADPKD) caused by PKD1 mutations is one of the most common hereditary disorders. However, the key pathological processes underlying cyst development and exacerbation in pre-symptomatic stages remain unknown, because rodent models do not recapitulate critical disease phenotypes, including disease onset in heterozygotes. Here, using CRISPR/Cas9, we generate ADPKD models with PKD1 mutations in cynomolgus monkeys. As in humans and mice, near-complete PKD1 depletion induces severe cyst formation mainly in collecting ducts. Importantly, unlike in mice, PKD1 heterozygote monkeys exhibit cyst formation perinatally in distal tubules, possibly reflecting the initial pathology in humans. Many monkeys in these models survive after cyst formation, and cysts progress with age. Furthermore, we succeed in generating selective heterozygous mutations using allele-specific targeting. We propose that our models elucidate the onset and progression of ADPKD, which will serve as a critical basis for establishing new therapeutic strategies, including drug treatments.

The Cynomolgus Macaque MHC Polymorphism in Experimental Medicine

Among the non-human primates used in experimental medicine, cynomolgus macaques (Macaca fascicularis hereafter referred to as Mafa) are increasingly selected for the ease with which they are maintained and bred in captivity. Macaques belong to Old World monkeys and are phylogenetically much closer to humans than rodents, which are still the most frequently used animal model. Our understanding of the Mafa genome has progressed rapidly in recent years and has greatly benefited from the latest technical advances in molecular genetics. Cynomolgus macaques are widespread in Southeast Asia and numerous studies have shown a distinct genetic differentiation of continental and island populations. The major histocompatibility complex of cynomolgus macaque (Mafa MHC) is organized in the same way as that of human, but it differs from the latter by its high degree of classical class I gene duplication. Human polymorphic MHC regions play a pivotal role in allograft transplantation and have been associated with more than 100 diseases and/or phenotypes. The Mafa MHC polymorphism similarly plays a crucial role in experimental allografts of organs and stem cells. Experimental results show that the Mafa MHC class I and II regions influence the ability to mount an immune response against infectious pathogens and vaccines. MHC also affects cynomolgus macaque reproduction and impacts on numerous biological parameters. This review describes the Mafa MHC polymorphism and the methods currently used to characterize it. We discuss some of the major areas of experimental medicine where an effect induced by MHC polymorphism has been demonstrated.

CRISPR/Cas9-mediated disruption of SHANK3 in monkey leads to drug-treatable autism-like symptoms

Monogenic mutations in the SHANK3 gene, which encodes a postsynaptic scaffold protein, play a causative role in autism spectrum disorder (ASD). Although a number of mouse models with Shank3 mutations have been valuable for investigating the pathogenesis of ASD, species-dependent differences in behaviors and brain structures post considerable challenges to use small animals to model ASD and to translate experimental therapeutics to the clinic. We have used clustered regularly interspersed short palindromic repeat/CRISPR-associated nuclease 9 to generate a cynomolgus monkey model by disrupting SHANK3 at exons 6 and 12. Analysis of the live mutant monkey revealed the core behavioral abnormalities of ASD, including impaired social interaction and repetitive behaviors, and reduced brain network activities detected by positron-emission computed tomography (PET). Importantly, these abnormal behaviors and brain activities were alleviated by the antidepressant fluoxetine treatment. Our findings provide the first demonstration that the genetically modified non-human primate can be used for translational research of therapeutics for ASD.

DNA-based Determination of Ancestry in Cynomolgus Macaques (Macaca fascicularis)

Interest in the genetic composition of cynomolgus macaques (Macaca fascicularis) has increased due to the rising demand for NHP models in human biomedical research. Significant genetic differences among regional populations of cynomolgus macaques can confound interpretations of research results because they do not solely reflect differences in experimental treatment effects. Therefore, the common origin of cynomolgus macaques used as research subjects should be verified by using region-specific genetic markers to minimize the influence of underlying genetic variation among animals selected as research subjects on phenotypes under study. We compared the effectiveness of 18 short tandem repeat (STR) markers with that of 83 single-nucleotide polymorphism (SNP) markers to differentiate the ancestry of cynomolgus macaques from 6 different populations (Cambodia, Sumatra, Mauritius, Singapore, and the islands of Luzon and Zamboanga in the Philippines). Genetic diversity indices such as allele numbers and expected heterozygosity based on SNP were lower and exhibited lower standard errors than those provided by STR, probably because, unlike STR, most SNP are biallelic and consequently exhibit maximal expected heterozygosity values of 0.50. However, the standard error of estimates of observed heterozygosity based on SNP was higher than that for STR, perhaps reflecting sampling errors. Only 27 SNP were required to match the resolving power of 17 STR to detect population structure, that is, 1.6 SNP:1 STR. Whereas STR only differentiated the Mauritian population from all other populations, SNP detected 4 genetically distinct groups (Cambodia, Singapore-Sumatra, Mauritius, and Zamboanga). SNP are poised to become as valuable as STR for understanding and detecting genetic structure among cynomolgus macaques. Although STR will remain an important tool for cynomolgus macaque population studies, SNP have the potential to become the mainstream marker type.

Characterisation of the hepatitis B virus cross-species transmission pattern via Na+/taurocholate co-transporting polypeptides from 11 New World and Old World primate species

The hepatic Na+/taurocholate co-transporting polypeptide (NTCP in man, Ntcp in animals) is the high-affinity receptor for the hepatitis B (HBV) and hepatitis D (HDV) viruses. Species barriers for human HBV/HDV within the order Primates were previously attributed to Ntcp sequence variations that disable virus-receptor interaction. However, only a limited number of primate Ntcps have been analysed so far. In the present study, a total of 11 Ntcps from apes, Old and New World monkeys were cloned and expressed in vitro to characterise their interaction with HBV and HDV. All Ntcps showed intact bile salt transport. Human NTCP as well as the Ntcps from the great apes chimpanzee and orangutan showed transport-competing binding of HBV derived myr-preS1-peptides. In contrast, all six Ntcps from the group of Old World monkeys were insensitive to HBV myr-preS1-peptide binding and HBV/HDV infection. This is basically predetermined by the amino acid arginine at position 158 of all studied Old World monkey Ntcps. An exchange from arginine to glycine (as present in humans and great apes) at this position (R158G) alone was sufficient to achieve full transport-competing HBV myr-preS1-peptide binding and susceptibility for HBV/HDV infection. New World monkey Ntcps showed higher sequence heterogeneity, but in two cases with 158G showed transport-competing HBV myr-preS1-peptide binding, and in one case (Saimiri sciureus) even susceptibility for HBV/HDV infection. In conclusion, amino acid position 158 of NTCP/Ntcp is sufficient to discriminate between the HBV/HDV susceptible group of humans and great apes (158G) and the non-susceptible group of Old World monkeys (158R). In the case of the phylogenetically more distant New World monkey Ntcps amino acid 158 plays a significant, but not exclusive role.

Analysis of nuclear and organellar genomes of Plasmodium knowlesi in humans reveals ancient population structure and recent recombination among host-specific subpopulations

The macaque parasite Plasmodium knowlesi is a significant concern in Malaysia where cases of human infection are increasing. Parasites infecting humans originate from genetically distinct subpopulations associated with the long-tailed (Macaca fascicularis (Mf)) or pig-tailed macaques (Macaca nemestrina (Mn)). We used a new high-quality reference genome to re-evaluate previously described subpopulations among human and macaque isolates from Malaysian-Borneo and Peninsular-Malaysia. Nuclear genomes were dimorphic, as expected, but new evidence of chromosomal-segment exchanges between subpopulations was found. A large segment on chromosome 8 originating from the Mn subpopulation and containing genes encoding proteins expressed in mosquito-borne parasite stages, was found in Mf genotypes. By contrast, non-recombining organelle genomes partitioned into 3 deeply branched lineages, unlinked with nuclear genomic dimorphism. Subpopulations which diverged in isolation have re-connected, possibly due to deforestation and disruption of wild macaque habitats. The resulting genomic mosaics reveal traits selected by host-vector-parasite interactions in a setting of ecological transition.

Transgenic Mice Overexpressing the Divalent Metal Transporter 1 Exhibit Iron Accumulation and Enhanced Parkin Expression in the Brain

Exposure to divalent metals such as iron and manganese is thought to increase the risk for Parkinson's disease (PD). Under normal circumstances, cellular iron and manganese uptake is regulated by the divalent metal transporter 1 (DMT1). Accordingly, alterations in DMT1 levels may underlie the abnormal accumulation of metal ions and thereby disease pathogenesis. Here, we have generated transgenic mice overexpressing DMT1 under the direction of a mouse prion promoter and demonstrated its robust expression in several regions of the brain. When fed with iron-supplemented diet, DMT1-expressing mice exhibit rather selective accumulation of iron in the substantia nigra, which is the principal region affected in human PD cases, but otherwise appear normal. Alongside this, the expression of Parkin is also enhanced, likely as a neuroprotective response, which may explain the lack of phenotype in these mice. When DMT1 is overexpressed against a Parkin null background, the double-mutant mice similarly resisted a disease phenotype even when fed with iron- or manganese-supplemented diet. However, these mice exhibit greater vulnerability toward 6-hydroxydopamine-induced neurotoxicity. Taken together, our results suggest that iron accumulation alone is not sufficient to cause neurodegeneration and that multiple hits are required to promote PD.

Major histocompatibility complex haplotyping and long-amplicon allele discovery in cynomolgus macaques from Chinese breeding facilities

Very little is currently known about the major histocompatibility complex (MHC) region of cynomolgus macaques (Macaca fascicularis; Mafa) from Chinese breeding centers. We performed comprehensive MHC class I haplotype analysis of 100 cynomolgus macaques from two different centers, with animals from different reported original geographic origins (Vietnamese, Cambodian, and Cambodian/Indonesian mixed-origin). Many of the samples were of known relation to each other (sire, dam, and progeny sets), making it possible to characterize lineage-level haplotypes in these animals. We identified 52 Mafa-A and 74 Mafa-B haplotypes in this cohort, many of which were restricted to specific sample origins. We also characterized full-length MHC class I transcripts using Pacific Biosciences (PacBio) RS II single-molecule real-time (SMRT) sequencing. This technology allows for complete read-through of unfragmented MHC class I transcripts (~1100 bp in length), so no assembly is required to unambiguously resolve novel full-length sequences. Overall, we identified 311 total full-length transcripts in a subset of 72 cynomolgus macaques from these Chinese breeding facilities; 130 of these sequences were novel and an additional 115 extended existing short database sequences to span the complete open reading frame. This significantly expands the number of Mafa-A, Mafa-B, and Mafa-I full-length alleles in the official cynomolgus macaque MHC class I database. The PacBio technique described here represents a general method for full-length allele discovery and genotyping that can be extended to other complex immune loci such as MHC class II, killer immunoglobulin-like receptors, and Fc gamma receptors.

Ancestry, Plasmodium cynomolgi prevalence and rhesus macaque admixture in cynomolgus macaques (Macaca fascicularis) bred for export in Chinese breeding farms

BACKGROUND

Most cynomolgus macaques (Macaca fascicularis) used in the United States as animal models are imported from Chinese breeding farms without documented ancestry. Cynomolgus macaques with varying rhesus macaque ancestry proportions may exhibit differences, such as susceptibility to malaria, that affect their suitability as a research model.

METHODS

DNA of 400 cynomolgus macaques from 10 Chinese breeding farms was genotyped to characterize their regional origin and rhesus ancestry proportion. A nested PCR assay was used to detect Plasmodium cynomolgi infection in sampled individuals.

RESULTS

All populations exhibited high levels of genetic heterogeneity and low levels of inbreeding and genetic subdivision. Almost all individuals exhibited an Indochinese origin and a rhesus ancestry proportion of 5%-48%. The incidence of P. cynomolgi infection in cynomolgus macaques is strongly associated with proportion of rhesus ancestry.

CONCLUSIONS

The varying amount of rhesus ancestry in cynomolgus macaques underscores the importance of monitoring their genetic similarity in malaria research.

Alu-Derived Alternative Splicing Events Specific to Macaca Lineages in CTSF Gene

Cathepsin F, which is encoded by CTSF, is a cysteine proteinase ubiquitously expressed in several tissues. In a previous study, novel transcripts of the CTSF gene were identified in the crab-eating monkey deriving from the integration of an Alu element-AluYRa1. The occurrence of AluYRa1-derived alternative transcripts and the mechanism of exonization events in the CTSF gene of human, rhesus monkey, and crab-eating monkey were investigated using PCR and reverse transcription PCR on the genomic DNA and cDNA isolated from several tissues. Results demonstrated that AluYRa1 was only integrated into the genome of Macaca species and this lineage-specific integration led to exonization events by producing a conserved 3' splice site. Six transcript variants (V1-V6) were generated by alternative splicing (AS) events, including intron retention and alternative 5' splice sites in the 5' and 3' flanking regions of CTSF_AluYRa1. Among them, V3-V5 transcripts were ubiquitously expressed in all tissues of rhesus monkey and crab-eating monkey, whereas AluYRa1-exonized V1 was dominantly expressed in the testis of the crab-eating monkey, and V2 was only expressed in the testis of the two monkeys. These five transcript variants also had different amino acid sequences in the C-terminal region of CTSF, as compared to reference sequences. Thus, species-specific Alu-derived exonization by lineage-specific integration of Alu elements and AS events seems to have played an important role during primate evolution by producing transcript variants and gene diversification.

Catalysis and Structure of Zebrafish Urate Oxidase Provide Insights into the Origin of Hyperuricemia in Hominoids

Urate oxidase (Uox) catalyses the first reaction of oxidative uricolysis, a three-step enzymatic pathway that allows some animals to eliminate purine nitrogen through a water-soluble compound. Inactivation of the pathway in hominoids leads to elevated levels of sparingly soluble urate and puts humans at risk of hyperuricemia and gout. The uricolytic activities lost during evolution can be replaced by enzyme therapy. Here we report on the functional and structural characterization of Uox from zebrafish and the effects on the enzyme of the missense mutation (F216S) that preceded Uox pseudogenization in hominoids. Using a kinetic assay based on the enzymatic suppression of the spectroscopic interference of the Uox reaction product, we found that the F216S mutant has the same turnover number of the wild-type enzyme but a much-reduced affinity for the urate substrate and xanthine inhibitor. Our results indicate that the last functioning Uox in hominoid evolution had an increased Michaelis constant, possibly near to upper end of the normal range of urate in the human serum (~300 μM). Changes in the renal handling of urate during primate evolution can explain the genetic modification of uricolytic activities in the hominoid lineage without the need of assuming fixation of deleterious mutations.

Gain of a New Exon by a Lineage-Specific Alu Element-Integration Event in the BCS1L Gene during Primate Evolution

BCS1L gene encodes mitochondrial protein and is a member of conserved AAA protein family. This gene is involved in the incorporation of Rieske FeS and Qcr10p into complex III of respiratory chain. In our previous study, AluYRa2-derived alternative transcript in rhesus monkey genome was identified. However, this transcript has not been reported in human genome. In present study, we conducted evolutionary analysis of AluYRa2-exonized transcript with various primate genomic DNAs and cDNAs from humans, rhesus monkeys, and crab-eating monkeys. Remarkably, our results show that AluYRa2 element has only been integrated into genomes of Macaca species. This Macaca lineage-specific integration of AluYRa2 element led to exonization event in the first intron region of BCS1L gene by producing a conserved 3' splice site. Intriguingly, in rhesus and crab-eating monkeys, more diverse transcript variants by alternative splicing (AS) events, including exon skipping and different 5' splice sites from humans, were identified. Alignment of amino acid sequences revealed that AluYRa2-exonized transcript has short N-terminal peptides. Therefore, AS events play a major role in the generation of various transcripts and proteins during primate evolution. In particular, lineage-specific integration of Alu elements and species-specific Alu-derived exonization events could be important sources of gene diversification in primates.

Analysis of full-length mitochondrial DNA D-loop sequences from Macaca fascicularis of different geographical origin reveals novel haplotypes

BACKGROUND

Cynomolgus macaques are indigenous to Asia occupying a range of geographical areas. A non-indigenous population established on Mauritius approximately 500 years ago. Mauritian cynomolgus macaques are recognised as having low genetic diversity compared to Indonesian macaques, from which they originated. As cynomolgus macaques are widely used as a biomedical model, there have been many studies of their genetic relationships. However, population diversity and relationships have only been assessed through analysis of either the hypervariable region I or II separately within the D-loop region of the mitochondrial genome in these macaques.

METHODS

Using sequencing, we defined haplotypes encompassing the full D-loop sequence for Mauritian and Indonesian cynomolgus macaques.

RESULTS

We evaluated the haplotype relationships by constructing a median-joining network based on full-length D-loop sequences, which has not been reported previously.

CONCLUSION

Our data allow a complete D-loop haplotype, including a hereto unreported polymorphic region, to be defined to aid the resolution of populations of cynomolgus macaques and which highlights the value in analysing both D-loop hypervariable regions in concert.

Comparison of Methods for Determining ABO Blood Type in Cynomolgus Macaques (Macaca fascicularis)

Thorough examination of ABO blood type in cynomolgus monkeys is an essential experimental step to prevent humoral rejection during transplantation research. In the present study, we evaluated current methods of ABO blood-antigen typing in cynomolgus monkeys by comparing the outcomes obtained by reverse hemagglutination, single-nucleotide polymorphism (SNP) analysis, and buccal mucosal immunohistochemistry. Among 21 animals, 5 were type A regardless of the method. However, of 8 serologically type B animals, 3 had a heterozygous type AB SNP profile, among which 2 failed to express A antigen, as shown by immunohistochemical analysis. Among 8 serologically type AB animals, 2 appeared to be type A by SNP analysis and immunohistochemistry. None of the methods identified any type O subjects. We conclude that the expression of ABO blood-group antigens is regulated by an incompletely understood process and that using both SNP and immunohistochemistry might minimize the risk of incorrect results obtained from the conventional hemagglutination assay.

Use of genome-wide heterospecific single-nucleotide polymorphisms to estimate linkage disequilibrium in rhesus and cynomolgus macaques

Rhesus and cynomolgus macaques are frequently used in biomedical research, and the availability of their reference genomes now provides for their use in genome-wide association studies. However, little is known about linkage disequilibrium (LD) in their genomes, which can affect the design and success of such studies. Here we studied LD by using 1781 conserved single-nucleotide polymorphisms (SNPs) in 183 rhesus macaques (Macaca mulatta), including 97 purebred Chinese and 86 purebred Indian animals, and 96 cynomolgus macaques (M. fascicularis fascicularis). Correlation between loci pairs decayed to 0.02 at 1146.83, 2197.92, and 3955.83 kb for Chinese rhesus, Indian rhesus, and cynomolgus macaques, respectively. Differences between the observed heterozygosity and minor allele frequency (MAF) of pairs of these 3 taxa were highly statistically significant. These 3 nonhuman primate taxa have significantly different genetic diversities (heterozygosity and MAF) and rates of LD decay. Our study confirms a much lower rate of LD decay in Indian than in Chinese rhesus macaques relative to that previously reported. In contrast, the especially low rate of LD decay in cynomolgus macaques suggests the particular usefulness of this species in genome-wide association studies. Although conserved markers, such as those used here, are required for valid LD comparisons among taxa, LD can be assessed with less bias by using species-specific markers, because conserved SNPs may be ancestral and therefore not informative for LD.

[Small hairpin RNA targeting inhibition of NF-κB gene in endometriosis therapy of Macaca fascicularis]

OBJECTIVE

To observe the therapeutic effect of NF-κB gene short hairpin RNA (shRNA) on endometriosis and identify the function of NF-κB on the maintenance and development of endometriosis in Macaca fascicularis.

METHODS

The Macaca fascicularis model of endometriosis was developed, which divided into experimental group, negative control group and simple model group. The high specificity adenovirus vector mediated shRNA targeting NF-κB gene and negative control shRNA adenovirus with no-load NF-κB gene were synthesised. The experimental group injected the adenovirus which carried the NF-κB shRNA into the endometriosis lesions under laparoscopy surgery, the negative control group with no-load shRNA adenovirus and the simple models group injected with normal saline. Four weeks later after the injection, an observed operation was performed through laparoscopy and some lesions were collected. The CD34 immunohistochemistry of these lesions were done to detect the microvessel density, then the variation of the microvessel density among each group were observed. The expression of the NF-κB and proliferating cell nuclear antigen (PCNA) were detected through western blot.

RESULTS

First, the Macaca fascicularis model of endometriosis was successful developed, and the experimental group has an evident atrophy in ectopic lesions compared with the previous. The lesions' microvessel density in experimental group decreased evidently compared with the negative control group and simple model group (0.002 0±0.000 3 versus 0.021 9±0.002 6 versus 0.024 5±0.003 3), and the differences was statistically significant (P < 0.01). The expression of PCNA (0.37±0.17 versus 0.57±0.26 versus 0.57±0.28) and NF-κB (0.338±0.174 versus 0.678±0.021 versus 0.645±0.098) in experiment group was lower than the negative control group and simple model group, the differences were statistically significant (all P < 0.01).

CONCLUSION

Through targeting suppressed the NF-κB gene expression by NF-κB shRNA, we can inhibit the development of endometriosis through reducing the ability of angiogenesis and cell proliferation of ectopic endometrial cells.

Genetic substructure in cynomolgus macaques (Macaca fascicularis) on the island of Mauritius

BACKGROUND

Nonhuman primates are commonly used in biomedical research as animal models of human disease and behavior. Compared to common rodent models, nonhuman primates are genetically, physiologically, behaviorally and neurologically more similar to humans owing to more recent shared ancestry and therefore provide the advantage of greater translational validity in preclinical studies. The cynomolgus macaque (Macaca fascicularis) is one of the most commonly used nonhuman primates in academic and industry settings, yet population genetic research has revealed significant substructure throughout the species distribution that may confound studies. Cynomolgus monkeys introduced to Mauritius specifically have previously been thought to maintain the least genetic heterogeneity of all cynomolgus monkeys, although recent work, including work from our lab, suggests macaques from Mauritius too may harbor cryptic substructure.

RESULTS

To evaluate putative substructure in Mauritian cynomolgus macaques, we designed a panel of 96 single nucleotide polymorphisms based on preliminary findings from previous work to screen 246 of cynomolgus monkeys from two primary suppliers. Results from this study support substructure in Mauritian macaques and suggest a minimum of two populations and maybe three on Mauritius, with moderate admixture.

CONCLUSION

These findings inform the natural history of these monkeys suggesting either a previously unrecognized physical or ecological barrier to gene flow on Mauritius and/or the breakdown of historic substructure resulting from the history of macaque introduction to the island. These findings are relevant to ongoing research using these models in part because of increased appreciation of segregating common variation with functional effects and may be used to better inform animal selection in preclinical research.

Effect of MHC haplotype on immune response upon experimental SHIVSF162P4cy infection of Mauritian cynomolgus macaques

Little is known about the effects of Major Histocompatibility Complex (MHC) haplotypes on immunity to primate lentiviruses involving both acquired and innate immune responses. We present statistical evidence of the influence of MHC polymorphism on antiviral immunity of Mauritian cynomolgus macaques (MCM) following simian/human immunodeficiency virus SHIV_SF162P4cy infection, involving the production of pro- and anti-inflammatory cytokines and α-defensins, which may modulate acquired immune responses. During the acute phase of infection, IL-10 correlated positively with viral load and negatively with CD4+T cell counts. Furthermore, α-defensins production was directly correlated with plasma viral RNA, particularly at peak of viral load. When the effects of the MHC were analyzed, a significant association between lower anti-Env binding and neutralizing antibody levels with class IB M4 haplotype and with class IA, IB M4 haplotype, respectively, was observed in the post-acute phase. Lower antibody responses may have resulted into a poor control of infection thus explaining the previously reported lower CD4 T cell counts in these monkeys. Class II M3 haplotype displayed significantly lower acute and post-acute IL-10 levels. In addition, significantly lower levels of α-defensins were detected in class IA M3 haplotype monkeys than in non-M3 macaques, in the post-acute phase of infection. These data indicate that the MHC could contribute to the delicate balance of pro-inflammatory mechanisms, particularly with regard to the association between IL-10 and α-defensins in lentivirus infection. Our results show that host genetic background, virological and immunological parameters should be considered for the design and interpretation of HIV-1 vaccine efficacy studies.

Evidence for a role for interleukin-17, Th17 cells and iron homeostasis in protective immunity against tuberculosis in cynomolgus macaques

Tuberculosis (TB) remains a major global public health problem. The only vaccine, BCG, gives variable protection, especially in adults, so several new vaccines are in clinical trials. There are no correlates of protective immunity to TB; therefore vaccines progress through lengthy and expensive pre-clinical assessments and human trials. Correlates of protection could act as early end-points during clinical trials, accelerating vaccine development and reducing costs. A genome-wide microarray was utilised to identify potential correlates of protection and biomarkers of disease induced post-BCG vaccination and post-Mycobacterium tuberculosis challenge in PPD-stimulated peripheral blood mononuclear cells from cynomolgus macaques where the outcome of infection was known. Gene expression post BCG-vaccination and post challenge was compared with gene expression when the animals were naïve. Differentially expressed genes were identified using a moderated T test with Benjamini Hochberg multiple testing correction. After BCG vaccination and six weeks post-M. tuberculosis challenge, up-regulation of genes related to a Th1 and Th17 response was observed in disease controllers. At post-mortem, RT-PCR revealed an up-regulation of iron regulatory genes in animals that developed TB and down-regulation of these genes in disease controllers, indicating the ability to successfully withhold iron may be important in the control of TB disease. The induction of a balanced Th1 and Th17 response, together with expression of effector cytokines, such as IFNG, IL2, IL17, IL21 and IL22, could be used as correlates of a protective host response.

Identification and characterization of microRNAs in the crab-eating macaque (Macaca fascicularis) using transcriptome analysis

MicroRNAs (miRNAs), with an average length between 16 nt and 26 nt, are small non-coding RNAs that can repress gene expression on the post-transcriptional level. Macaca fascicularis (M. fascicularis), one of the most important nonhuman primate animal models, is widely used in basic and applied preclinical research, especially in studies that involve neuroscience and disease. However, due to the lack of a complete genome sequence, the miRNAs in M. fascicularis have not been completely characterized. In this study, 86 putative M. fascicularis miRNAs were identified using a strategy of our design. The expression of some of these miRNAs in the tissue was confirmed by qRT-PCR. The function and pathway of their targeted genes were analyzed to reveal the potential relevance of miRNA regulation on diseases and physiological processes. The current study provides insight into potential miRNAs and forms a useful knowledge base for the future understanding of the function of miRNAs in M. fascicularis.

Finding the factors of reduced genetic diversity on X chromosomes of Macaca fascicularis: male-driven evolution, demography, and natural selection

The ratio of genetic diversity on X chromosomes relative to autosomes in organisms with XX/XY sex chromosomes could provide fundamental insight into the process of genome evolution. Here we report this ratio for 24 cynomolgus monkeys (Macaca fascicularis) originating in Indonesia, Malaysia, and the Philippines. The average X/A diversity ratios in these samples was 0.34 and 0.20 in the Indonesian-Malaysian and Philippine populations, respectively, considerably lower than the null expectation of 0.75. A Philippine population supposed to derive from an ancestral population by founding events showed a significantly lower ratio than the parental population, suggesting a demographic effect for the reduction. Taking sex-specific mutation rate bias and demographic effect into account, expected X/A diversity ratios generated by computer simulations roughly agreed with the observed data in the intergenic regions. In contrast, silent sites in genic regions on X chromosomes showed strong reduction in genetic diversity and the observed X/A diversity ratio in the genic regions cannot be explained by mutation rate bias and demography, indicating that natural selection also reduces the level of polymorphism near genes. Whole-genome analysis of a female cynomolgus monkey also supported the notion of stronger reduction of genetic diversity near genes on the X chromosome.

TRIM5α polymorphism identification in cynomolgus macaques of Vietnamese origin and Chinese rhesus macaques

TRIM5α is a retroviral restriction factor, in which the B30.2 (SPRY) and coiled-coil domains cooperate to determine the specificity of TRIM5α-mediated capture of retroviral capsids. Here, all exons of TRIM5α were analyzed in 39 Vietnamese cynomolgus macaques (VCE) and 29 Chinese rhesus macaques (CR). Our results revealed the presence of 22 alleles using the PHASE 2.0 software package (PHylogenetics And Sequence Evolution), including two novel species-specific alleles with a low frequency in VCE in exons 4 and 8, which encoded the coiled-coil and B30.2 (SPRY) domains, respectively. Nine alleles were detected in both VCE and CR, while four alleles were likely shared between the species. Of these alleles, the highest frequencies of 38% and 26% occurred in VCE and CR, respectively. Importantly, we found that some alleles encoded the same coiled-coil domain, but not the SPRY domain. In contrast, other alleles encoded the same SPRY domain, but not the coiled-coil domain. Our findings will contribute to the understanding of the interaction between the two domains and the determination of the specificity of TRIM5α-mediated capture of retroviral capsids. Our results from the phylogenetic trees constructed for VCE and CR suggested that the macaques' ability to inhibit SIV replication became gradually stronger if they carried corresponding alleles in four clades (clades4-7). More interesting, in clade3, both novel allele pairs (4E100a, 10E147a) and allele pairs (7R17b and 13R11b), which had the strong ability to inhibit SIV replication, originated from the same ancestral allele, suggesting that the novel alleles might play a key role to determine an animal's ability to inhibit SIV/HIV replication. However, further studies are needed to increase our understanding of the genetic background of TRIM5α in these two macaque species.

The genetic composition of populations of cynomolgus macaques (Macaca fascicularis) used in biomedical research

BACKGROUND

The genetic composition of cynomolgus macaques used in biomedical research is not as well-characterized as that of rhesus macaques.

METHODS

Populations of cynomolgus macaques from Sumatra, Corregidor, Mauritius, Singapore, Cambodia, and Zamboanga were analyzed using 24 STRs.

RESULTS

The Sumatran and Cambodian populations exhibited the highest allelic diversity, while the Mauritian population exhibited the lowest. Sumatran cynomolgus macaques were the most genetically similar to all others, consistent with an Indonesian origin of the species. The high diversity among Cambodian animals may result from interbreeding with rhesus macaques. The Philippine and Mauritian samples were the most divergent from other populations, the former due to separation from the Sunda Shelf by deepwater and the latter due to anthropogenic translocation and extreme founder effects.

CONCLUSIONS

Investigators should verify their research subjects' origin, ancestry, and pedigree to minimize risks to biomedical experimentation from genetic variance stemming from close kinship and mixed ancestry as these can obscure treatment effects.

Selection of new appropriate reference genes for RT-qPCR analysis via transcriptome sequencing of cynomolgus monkeys (Macaca fascicularis)

In the investigation of the expression levels of target genes, reverse transcription quantitative real-time polymerase chain reaction (RT-qPCR) is the most accurate and widely used method. However, a normalization step is a prerequisite to obtain accurate quantification results from RT-qPCR data. Therefore, many studies regarding the selection of reference genes have been carried out. Recently, these studies have involved large-scale gene analysis methods such as microarray and next generation sequencing. In our previous studies, we analyzed large amounts of transcriptome data from the cynomolgus monkey. Using a modification of this large-scale transcriptome sequencing dataset, we selected and compared 12 novel candidate reference genes (ARFGAP2, ARL1, BMI1, CASC3, DDX3X, MRFAP1, ORMDL1, RSL24D1, SAR1A, USP22, ZC3H11A, and ZRANB2) and 4 traditionally used reference genes (ACTB, GAPDH, RPS19, and YWHAZ) in 13 different whole-body tissues by the 3 well-known programs geNorm, NormFinder, and BestKeeper. Combined analysis by these 3 programs showed that ADP-ribosylation factor GTPase activating protein 2 (ARFGAP2), morf4 family associated protein 1 (MRFAP1), and ADP-ribosylation factor-like 1 (ARL1) are the most appropriate reference genes for accurate normalization. Interestingly, 4 traditionally used reference genes were the least stably expressed in this study. For this reason, selection of appropriate reference genes is vitally important, and large-scale analysis is a good method for finding new candidate reference genes. Our results could provide reliable reference gene lists for future studies on the expression of various target genes in the cynomolgus monkey.

Selection of appropriate reference genes for RT-qPCR analysis in a streptozotocin-induced Alzheimer's disease model of cynomolgus monkeys (Macaca fascicularis)

Reverse transcription quantitative real-time polymerase chain reaction (RT-qPCR) has been widely used to quantify relative gene expression because of the specificity, sensitivity, and accuracy of this technique. In order to obtain reliable gene expression data from RT-qPCR experiments, it is important to utilize optimal reference genes for the normalization of target gene expression under varied experimental conditions. Previously, we developed and validated a novel icv-STZ cynomolgus monkey model for Alzheimer's disease (AD) research. However, in order to enhance the reliability of this disease model, appropriate reference genes must be selected to allow meaningful analysis of the gene expression levels in the icv-STZ cynomolgus monkey brain. In this study, we assessed the expression stability of 9 candidate reference genes in 2 matched-pair brain samples (5 regions) of control cynomolgus monkeys and those who had received intracerebroventricular injection of streptozotocin (icv-STZ). Three well-known analytical programs geNorm, NormFinder, and BestKeeper were used to choose the suitable reference genes from the total sample group, control group, and icv-STZ group. Combination analysis of the 3 different programs clearly indicated that the ideal reference genes are RPS19 and YWHAZ in the total sample group, GAPDH and RPS19 in the control group, and ACTB and GAPDH in the icv-STZ group. Additionally, we validated the normalization accuracy of the most appropriate reference genes (RPS19 and YWHAZ) by comparison with the least stable gene (TBP) using quantification of the APP and MAPT genes in the total sample group. To the best of our knowledge, this research is the first study to identify and validate the appropriate reference genes in cynomolgus monkey brains. These findings provide useful information for future studies involving the expression of target genes in the cynomolgus monkey.

Population and landscape genetics of an introduced species (M. fascicularis) on the island of Mauritius

The cynomolgus macaque, Macaca fascicularis, was introduced onto the island of Mauritius in the early 17(th) century. The species experienced explosive population growth, and currently exists at high population densities. Anecdotes collected from nonhuman primate trappers on the island of Mauritius allege that animals from the northern portion of the island are larger in body size than and superior in condition to their conspecifics in the south. Although previous genetic studies have reported Mauritian cynomolgus macaques to be panmictic, the individuals included in these studies were either from the southern/central or an unknown portion of the island. In this study, we sampled individuals broadly throughout the entire island of Mauritius and used spatial principle component analysis to measure the fine-scale correlation between geographic and genetic distance in this population. A stronger correlation between geographic and genetic distance was found among animals in the north than in those in the southern and central portions of the island. We found no difference in body weight between the two groups, despite anecdotal evidence to the contrary. We hypothesize that the increased genetic structure among populations in the north is related to a reduction in dispersal distance brought about by human habitation and tourist infrastructure, but too recent to have produced true genetic differentiation.

SIV genome-wide pyrosequencing provides a comprehensive and unbiased view of variation within and outside CD8 T lymphocyte epitopes

Deep sequencing technology is revolutionizing our understanding of HIV/SIV evolution. It is known that acute SIV sequence variation within CD8 T lymphocyte (CD8-TL) epitopes is similar among MHC-identical animals, but we do not know whether this persists into the chronic phase. We now determine whether chronic viral variation in MHC-identical animals infected with clonal SIV is similar throughout the entire coding sequence when using a sensitive deep sequencing approach. We pyrosequenced the entire coding sequence of the SIV genome isolated from a unique cohort of four SIVmac239-infected, MHC-identical Mauritian cynomolgus macaques (MCM) 48 weeks after infection; one MCM in the cohort became an elite controller. Among the three non-controllers, we found that genome-wide sequences were similar between animals and we detected increased sequence complexity within 64% of CD8-TL epitopes when compared to Sanger sequencing methods. When we compared sequences between the MHC-matched controller and the three non-controllers, we found the viral population in the controller was less diverse and accumulated different variants than the viral populations in the non-controllers. Importantly, we found that initial PCR amplification of viral cDNA did not significantly affect the sequences detected, suggesting that data obtained by pyrosequencing PCR-amplified viral cDNA accurately represents the diversity of sequences replicating within an animal. This demonstrates that chronic sequence diversity across the entire SIV coding sequence is similar among MHC-identical animals with comparable viral loads when infected with the same clonal virus stock. Additionally, our approach to genome-wide SIV sequencing accurately reflects the diversity of sequences present in the replicating viral population. In sum, our study suggests that genome-wide pyrosequencing of immunodeficiency viruses captures a thorough and unbiased picture of sequence diversity, and may be a useful approach to employ when evaluating which sequences to include as part of a vaccine immunogen.