Conserved and divergent gene regulatory programs of the mammalian neocortex.

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

Evolutionary and biomedical insights from a marmoset diploid genome assembly.

The accurate and complete assembly of both haplotype sequences of a diploid organism is essential to understanding the role of variation in genome functions, phenotypes and diseases1. Here, using a trio-binning approach, we present a high-quality, diploid reference genome, with both haplotypes assembled independently at the chromosome level, for the common marmoset (Callithrix jacchus), an primate model system that is widely used in biomedical research2,3. The full spectrum of heterozygosity between the two haplotypes involves 1.36% of the genome-much higher than the 0.13% indicated by the standard estimation based on single-nucleotide heterozygosity alone. The de novo mutation rate is 0.43 × 10-8 per site per generation, and the paternal inherited genome acquired twice as many mutations as the maternal. Our diploid assembly enabled us to discover a recent expansion of the sex-differentiation region and unique evolutionary changes in the marmoset Y chromosome. In addition, we identified many genes with signatures of positive selection that might have contributed to the evolution of Callithrix biological features. Brain-related genes were highly conserved between marmosets and humans, although several genes experienced lineage-specific copy number variations or diversifying selection, with implications for the use of marmosets as a model system.

Comparative cellular analysis of motor cortex in human, marmoset and mouse.

The primary motor cortex (M1) is essential for voluntary fine-motor control and is functionally conserved across mammals1. Here, using high-throughput transcriptomic and epigenomic profiling of more than 450,000 single nuclei in humans, marmoset monkeys and mice, we demonstrate a broadly conserved cellular makeup of this region, with similarities that mirror evolutionary distance and are consistent between the transcriptome and epigenome. The core conserved molecular identities of neuronal and non-neuronal cell types allow us to generate a cross-species consensus classification of cell types, and to infer conserved properties of cell types across species. Despite the overall conservation, however, many species-dependent specializations are apparent, including differences in cell-type proportions, gene expression, DNA methylation and chromatin state. Few cell-type marker genes are conserved across species, revealing a short list of candidate genes and regulatory mechanisms that are responsible for conserved features of homologous cell types, such as the GABAergic chandelier cells. This consensus transcriptomic classification allows us to use patch-seq (a combination of whole-cell patch-clamp recordings, RNA sequencing and morphological characterization) to identify corticospinal Betz cells from layer 5 in non-human primates and humans, and to characterize their highly specialized physiology and anatomy. These findings highlight the robust molecular underpinnings of cell-type diversity in M1 across mammals, and point to the genes and regulatory pathways responsible for the functional identity of cell types and their species-specific adaptations.

Genetic characterization of a captive marmoset (Callithrix jacchus) colony using genotype-by-sequencing.

The marmoset is a fundamental nonhuman primate model for the study of aging, neurobiology, and many other topics. Genetic management of captive marmoset colonies is complicated by frequent chimerism in the blood and other tissues, a lack of tools to enable cost-effective, genome-wide interrogation of variation, and historic mergers and migrations of animals between colonies. We implemented genotype-by-sequencing (GBS) of hair follicle derived DNA (a minimally chimeric DNA source) of 82 marmosets housed at the Southwest National Primate Research Center (SNPRC). Our primary goals were the genetic characterization of our marmoset population for pedigree verification and colony management and to inform the scientific community of the functional genetic makeup of this valuable resource. We used the GBS data to reconstruct the genetic legacy of recent mergers between colonies, to identify genetically related animals whose relationships were previously unknown due to incomplete pedigree information, and to show that animals in the SNPRC colony appear to exhibit low levels of inbreeding. Of the >99,000 single-nucleotide variants (SNVs) that we characterized, >9800 are located within gene regions known to harbor pathogenic variants of clinical significance in humans. Overall, we show the combination of low-resolution (sparse) genotyping using hair follicle DNA is a powerful strategy for the genetic management of captive marmoset colonies and for identifying potential SNVs for the development of biomedical research models.

Cellular-resolution gene expression profiling in the neonatal marmoset brain reveals dynamic species- and region-specific differences.

Precise spatiotemporal control of gene expression in the developing brain is critical for neural circuit formation, and comprehensive expression mapping in the developing primate brain is crucial to understand brain function in health and disease. Here, we developed an unbiased, automated, large-scale, cellular-resolution in situ hybridization (ISH)-based gene expression profiling system (GePS) and companion analysis to reveal gene expression patterns in the neonatal New World marmoset cortex, thalamus, and striatum that are distinct from those in mice. Gene-ontology analysis of marmoset-specific genes revealed associations with catalytic activity in the visual cortex and neuropsychiatric disorders in the thalamus. Cortically expressed genes with clear area boundaries were used in a three-dimensional cortical surface mapping algorithm to delineate higher-order cortical areas not evident in two-dimensional ISH data. GePS provides a powerful platform to elucidate the molecular mechanisms underlying primate neurobiology and developmental psychiatric and neurological disorders.

The Verriest Lecture: Pathways to color in the eye and brain.

In common with the majority of New World monkeys, marmosets show polymorphic color vision by allelic variation of X-chromosome genes encoding opsin pigments in the medium/long wavelength range. Male marmosets are thus obligate dichromats ("red-green color blind"), whereas females carrying distinct alleles on X chromosomes show one of three trichromatic phenotypes. Marmosets thus represent a "natural knock-out" system enabling comparison of red-green color vision in dichromatic and trichromatic visual systems. Further, study of short-wave (blue) cone pathways in marmosets has provided insights into primitive visual pathways for depth perception and attention. These investigations represent a parallel line to clinical research on color vision defects that was pioneered in studies by Guy Verreist, whom we honor in this eponymous lecture.

Integrated analysis of long-noncoding RNA and circular RNA expression in Peste-des-Petits-Ruminants Virus (PPRV) infected marmoset B lymphocyte (B95a) cells.

Peste-des-Petits-Ruminants (PPR) or goat plague is an important viral disease of sheep and goats caused by the small ruminant morbilli virus or PPR virus (PPRV). Long non coding RNAs (lncRNA) and circular RNAs (circRNA) play a pivotal role in several biological processes including regulation of virus-host interactions. The present study explored the expression of lncRNA, circRNA and their functions in PPRV infected B-lymphocyte (B95a) cells. The results revealed a total of 4531 lncRNA and 2348 circRNA expression in both mock and PPRV infected samples. Analysis of differentially expressed (DE) RNA identified 123 DE-lncRNA and 39 DE-circRNA as significantly dysregulated. Functional analysis of cis-target genes of DE-lncRNA indicated activation of TCF dependent WNT signaling and PKN1 stimulated transcription process. Interactions (sponging) of microRNA (miRNA) revealed 344 DE-lncRNA-miRNA and 93 DE-circRNA-miRNA pairs. The competing endogenous RNA (ceRNA) network of lncRNA/circRNA-miRNA-mRNA in PPRV infected B95a cells was represented by 69 ceRNA pairs. We validated the DE-circRNA by targeted amplification and sequencing of back spliced junctions (BSJs). The present study revealed a profile of lncRNA, circRNA and their potential ceRNA network in PPRV infection. The results provide insight for better understanding of PPRV-host interactions.

Mitogenomic phylogeny of Callithrix with special focus on human transferred taxa.

BACKGROUND: Callithrix marmosets are a relatively young primate radiation, whose phylogeny is not yet fully resolved. These primates are naturally para- and allopatric, but three species with highly invasive potential have been introduced into the southeastern Brazilian Atlantic Forest by the pet trade. There, these species hybridize with each other and endangered, native congeners. We aimed here to reconstruct a robust Callithrix phylogeny and divergence time estimates, and identify the biogeographic origins of autochthonous and allochthonous Callithrix mitogenome lineages. We sequenced 49 mitogenomes from four species (C. aurita, C. geoffroyi, C. jacchus, C. penicillata) and anthropogenic hybrids (C. aurita x Callithrix sp., C. penicillata x C. jacchus, Callithrix sp. x Callithrix sp., C. penicillata x C. geoffroyi) via Sanger and whole genome sequencing. We combined these data with previously published Callithrix mitogenomes to analyze five Callithrix species in total. RESULTS: We report the complete sequence and organization of the C. aurita mitogenome. Phylogenetic analyses showed that C. aurita was the first to diverge within Callithrix 3.54 million years ago (Ma), while C. jacchus and C. penicillata lineages diverged most recently 0.5 Ma as sister clades. MtDNA clades of C. aurita, C. geoffroyi, and C. penicillata show intraspecific geographic structure, but C. penicillata clades appear polyphyletic. Hybrids, which were identified by phenotype, possessed mainly C. penicillata or C. jacchus mtDNA haplotypes. The biogeographic origins of mtDNA haplotypes from hybrid and allochthonous Callithrix were broadly distributed across natural Callithrix ranges. Our phylogenetic results also evidence introgression of C. jacchus mtDNA into C. aurita. CONCLUSION: Our robust Callithrix mitogenome phylogeny shows C. aurita lineages as basal and C. jacchus lineages among the most recent within Callithrix. We provide the first evidence that parental mtDNA lineages of anthropogenic hybrid and allochthonous marmosets are broadly distributed inside and outside of the Atlantic Forest. We also show evidence of cryptic hybridization between allochthonous Callithrix and autochthonous C. aurita. Our results encouragingly show that further development of genomic resources will allow to more clearly elucidate Callithrix evolutionary relationships and understand the dynamics of Callithrix anthropogenic introductions into the Brazilian Atlantic Forest.

An ancient retroviral RNA element hidden in mammalian genomes and its involvement in co-opted retroviral gene regulation.

BACKGROUND: Retroviruses utilize multiple unique RNA elements to control RNA processing and translation. However, it is unclear what functional RNA elements are present in endogenous retroviruses (ERVs). Gene co-option from ERVs sometimes entails the conservation of viral cis-elements required for gene expression, which might reveal the RNA regulation in ERVs. RESULTS: Here, we characterized an RNA element found in ERVs consisting of three specific sequence motifs, called SPRE. The SPRE-like elements were found in different ERV families but not in any exogenous viral sequences examined. We observed more than a thousand of copies of the SPRE-like elements in several mammalian genomes; in human and marmoset genomes, they overlapped with lineage-specific ERVs. SPRE was originally found in human syncytin-1 and syncytin-2. Indeed, several mammalian syncytin genes: mac-syncytin-3 of macaque, syncytin-Ten1 of tenrec, and syncytin-Car1 of Carnivora, contained the SPRE-like elements. A reporter assay revealed that the enhancement of gene expression by SPRE depended on the reporter genes. Mutation of SPRE impaired the wild-type syncytin-2 expression while the same mutation did not affect codon-optimized syncytin-2, suggesting that SPRE activity depends on the coding sequence. CONCLUSIONS: These results indicate multiple independent invasions of various mammalian genomes by retroviruses harboring SPRE-like elements. Functional SPRE-like elements are found in several syncytin genes derived from these retroviruses. This element may facilitate the expression of viral genes, which were suppressed due to inefficient codon frequency or repressive elements within the coding sequences. These findings provide new insights into the long-term evolution of RNA elements and molecular mechanisms of gene expression in retroviruses.

Single cell transcriptome analysis of human, marmoset and mouse embryos reveals common and divergent features of preimplantation development.

The mouse embryo is the canonical model for mammalian preimplantation development. Recent advances in single cell profiling allow detailed analysis of embryogenesis in other eutherian species, including human, to distinguish conserved from divergent regulatory programs and signalling pathways in the rodent paradigm. Here, we identify and compare transcriptional features of human, marmoset and mouse embryos by single cell RNA-seq. Zygotic genome activation correlates with the presence of polycomb repressive complexes in all three species, while ribosome biogenesis emerges as a predominant attribute in primate embryos, supporting prolonged translation of maternally deposited RNAs. We find that transposable element expression signatures are species, stage and lineage specific. The pluripotency network in the primate epiblast lacks certain regulators that are operative in mouse, but encompasses WNT components and genes associated with trophoblast specification. Sequential activation of GATA6, SOX17 and GATA4 markers of primitive endoderm identity is conserved in primates. Unexpectedly, OTX2 is also associated with primitive endoderm specification in human and non-human primate blastocysts. Our cross-species analysis demarcates both conserved and primate-specific features of preimplantation development, and underscores the molecular adaptability of early mammalian embryogenesis.

Gene expression in response to optical defocus of opposite signs reveals bidirectional mechanism of visually guided eye growth.

Myopia (nearsightedness) is the most common eye disorder, which is rapidly becoming one of the leading causes of vision loss in several parts of the world because of a recent sharp increase in prevalence. Nearwork, which produces hyperopic optical defocus on the retina, has been implicated as one of the environmental risk factors causing myopia in humans. Experimental studies have shown that hyperopic defocus imposed by negative power lenses placed in front of the eye accelerates eye growth and causes myopia, whereas myopic defocus imposed by positive lenses slows eye growth and produces a compensatory hyperopic shift in refractive state. The balance between these two optical signals is thought to regulate refractive eye development; however, the ability of the retina to recognize the sign of optical defocus and the composition of molecular signaling pathways guiding emmetropization are the subjects of intense investigation and debate. We found that the retina can readily distinguish between imposed myopic and hyperopic defocus, and identified key signaling pathways underlying retinal response to the defocus of different signs. Comparison of retinal transcriptomes in common marmosets exposed to either myopic or hyperopic defocus for 10 days or 5 weeks revealed that the primate retina responds to defocus of different signs by activation or suppression of largely distinct pathways. We also found that 29 genes differentially expressed in the marmoset retina in response to imposed defocus are localized within human myopia quantitative trait loci (QTLs), suggesting functional overlap between genes differentially expressed in the marmoset retina upon exposure to optical defocus and genes causing myopia in humans. These findings identify retinal pathways involved in the development of myopia, as well as potential new strategies for its treatment.

An improved de novo genome assembly of the common marmoset genome yields improved contiguity and increased mapping rates of sequence data.

BACKGROUND: The common marmoset (Callithrix jacchus) is one of the most studied primate model organisms. However, the marmoset genomes available in the public databases are highly fragmented and filled with sequence gaps, hindering research advances related to marmoset genomics and transcriptomics. RESULTS: Here we utilize single-molecule, long-read sequence data to improve and update the existing genome assembly and report a near-complete genome of the common marmoset. The assembly is of 2.79 Gb size, with a contig N50 length of 6.37 Mb and a chromosomal scaffold N50 length of 143.91 Mb, representing the most contiguous and high-quality marmoset genome up to date. Approximately 90% of the assembled genome was represented in contigs longer than 1 Mb, with approximately 104-fold improvement in contiguity over the previously published marmoset genome. More than 98% of the gaps from the previously published genomes were filled successfully, which improved the mapping rates of genomic and transcriptomic data on to the assembled genome. CONCLUSIONS: Altogether the updated, high-quality common marmoset genome assembly provide improvements at various levels over the previous versions of the marmoset genome assemblies. This will allow researchers working on primate genomics to apply the genome more efficiently for their genomic and transcriptomic sequence data.

NHP-immunome: A translational research-oriented database of non-human primate immune system proteins.

We are currently living the advent of a new age for medicine in which basic research is being quickly translated into marketable drugs, and the widespread access to genomics data is allowing the design and implementation of personalized solutions to medical conditions. Non-human primates (NHP) have gained an essential role in drug discovery and safety testing due to their close phylogenetic relationship to humans. In this study, a collection of well characterized genes of the human immune system was used to define the orthology-based immunome in four NHP species, with carefully curated annotations available based on multi-tissue RNA-seq datasets. A broad variation in the frequency of expressed protein isoforms was observed between species. Finally, this analysis also revealed the lack of expression of at least four different chemokines in new-world primates. In addition, transcripts corresponding to four genes including interleukin 12 subunit alpha were expressed in humans but no other primate species analyzed. Access to the non-human primate immunome is available in http://www.fidic.org.co:90/proyecto/.

Maternal weight affects placental DNA methylation of genes involved in metabolic pathways in the common marmoset monkey (Callithrix jacchus).

Accumulating evidence suggests that dysregulation of placental DNA methylation (DNAm) is a mechanism linking maternal weight during pregnancy to metabolic programming outcomes. The common marmoset, Callithrix jaccus, is a platyrrhine primate species that has provided much insight into studies of the primate placenta, maternal condition, and metabolic programming, yet the relationships between maternal weight and placental DNAm are unknown. Here, we report genome-wide DNAm from term marmoset placentas using reduced representation bisulfite sequencing. We identified 74 genes whose DNAm pattern is associated with maternal weight during gestation. These genes are predominantly involved in energy metabolism and homeostasis, including the regulation of glycolytic and lipid metabolic processes pathways.

Anomalous Pigmentation in Invasive and Native Marmosets, Callithrix jacchus, Callithrix penicillata (Primates, Callitrichidae), and Their Hybrids in Brazil.

Leucism is the lack or reduction in pigmentation in the most or parts of the body, but not in the eyes and body extremities. It is extremely rare in primates and has never been reported for Callithrix, a genus endemic to Brazil. We searched for individuals of Callithrix jacchus and C. penicillata with pigmentation anomalies in a systematic survey of three protected areas in the Atlantic Forest, within museum collections in Brazil, and opportunistically during field studies. Since 2008, we have recorded 8 individuals with leucism in small urban and periurban forest patches. Four were from native populations of C. penicillata in Cerrado savannahs and of C. jacchus in the Caatinga xeric scrubland, and 4 were from populations of hybrids between C. jacchus and C. penicillata in invaded areas in the coastal Atlantic Forest. We found no pigmentation abnormalities in museum specimens. We hypothesize that the observed leucism may be linked to inbreeding within the native range, but to hybridization within the invaded range, and discuss the likely ecological consequences to leucistic individuals.

Characterization of major histocompatibility complex-related molecule 1 sequence variants in non-human primates.

The major histocompatibility complex (MHC) class I-related molecule, MR1, presents vitamin B metabolites from bacteria and yeast to mucosal-associated invariant T (MAIT) cells. Despite the evolutionary conservation of MR1, we do not know whether different allele variants of MR1 exist within the nonhuman primate (NHP) populations that are commonly used for biomedical research. In this study, we identified 21 distinct MR1 nucleotide sequences representing 32 different alleles across five different NHP populations. The majority of the alleles conferring amino acid changes (allele variants) were found in or near the alpha-1 domain of the mature MR1 protein. We expressed four of the most commonly observed MR1 allele variants in 293T cells, and we found that each variant could present bacterial metabolites on the cell surface. We successfully induced cytokine production in macaque MAIT cells stimulated with 293T cells expressing the four most common MR1 allele variants, demonstrating the usefulness of these cell lines to study MAIT cell activity. Our data suggests that MR1 is not monomorphic, but that there are multiple MR1 alleles in NHPs. The materials we describe here will be valuable for characterizing differences in MR1 antigen presentation and MAIT cell function in NHPs.

Discovery of a new repeat family in the Callithrix jacchus genome.

We identified a novel repeat family, termed Platy-1, in the Callithrix jacchus (common marmoset) genome that arose around the time of the divergence of platyrrhines and catarrhines and established itself as a repeat family in New World monkeys (NWMs). A full-length Platy-1 element is ∼100 bp in length, making it the shortest known short interspersed element (SINE) in primates, and harbors features characteristic of non-LTR retrotransposons. We identified 2268 full-length Platy-1 elements across 62 subfamilies in the common marmoset genome. Our subfamily reconstruction and phylogenetic analyses support Platy-1 propagation throughout the evolution of NWMs in the lineage leading to C. jacchus Platy-1 appears to have reached its amplification peak in the common ancestor of current day marmosets and has since moderately declined. However, identification of more than 200 Platy-1 elements identical to their respective consensus sequence, and the presence of polymorphic elements within common marmoset populations, suggests ongoing retrotransposition activity. Platy-1, a SINE, appears to have originated from an Alu element, and hence is likely derived from 7SL RNA. Our analyses illustrate the birth of a new repeat family and its propagation dynamics in the lineage leading to the common marmoset over the last 40 million years.

Developmental and family history-based analysis of congenital fused labia phenotype in the captive common marmoset (Callithrix jacchus).

BACKGROUND: Congenital fused labia (CFL) is defined as a failure or significant delay in the opening of the juvenile sealed labia majora. This phenotype is known to be variably common in adult captive female marmosets but has never been investigated in detail before. MATERIALS AND METHODS: Here, we define, describe and quantify the variations in the degree of closure of the vulva in 122 captive marmosets (Callithrix jacchus) from 1.2 to 42 months old and include colony analysis. RESULTS: There was a negative correlation between the degree of labial fusion and animal age after prepubertal period (P < 0.05). CFL females had higher number CFL relatives (4.3 ± 0.6 vs 2.4 ± 0.5 for non-CFL, P < 0.05) and more external ancestors compared to non-CFL (P < 0.05). CONCLUSIONS: Our results therefore suggest that CFL phenotype is most likely associated with epigenetic effects induced by the captive environment and colony management strategy of extensive crossing of family lines to promote heterozygosity.

How many pygmy marmoset (Cebuella Gray, 1870) species are there? A taxonomic re-appraisal based on new molecular evidence.

The pygmy marmoset, Cebuella pygmaea, the smallest of the New World monkeys, has one of the largest geographical distributions of the Amazonian primates. Two forms have been recognized: Cebuella pygmaea pygmaea (Spix, 1823), and C. p. niveiventris Lönnberg, 1940. In this study, we investigated if the separation of pygmy marmosets into these two clades can be corroborated by molecular data. We also examine and compare coloration of the pelage in light of the new molecular results. We analyzed the mtDNA cytochrome b gene and, for the first time for any Neotropical primate, we used a reduced representation genome sequencing approach (ddRADseq) to obtain data for recently collected, geographically representative samples from the Rio Japurá, a northern tributary of the Rio Solimões and from the Javarí, Jutaí, Juruá, Madeira and Purus river basins, all tributaries south of the Solimões. We estimated phylogenies and diversification times under both maximum likelihood and Bayesian inference criteria. Our analysis showed two highly supported clades, with intraclade divergences much smaller than interclade divergences, indicating two species of Cebuella: one from the Rio Japurá and one to the south of Solimões. The interpretation of our results in light of the current taxonomy is not trivial however. Lönnberg stated that the type of Spix's pygmy marmoset (type locality 'near Tabatinga') was obtained from the south of the Solimões, and his description of the distinct niveiventris from Lago Ipixuna, south of the Solimões and several hundred kilometres east of Tabatinga, was based on a comparison with specimens that he determined as typical pygmaea that were from the upper Rio Juruá (south of the Solimões). As such it remains uncertain whether the name pygmaea should be applicable to the pygmy marmosets north of the Rio Solimões (Tabatinga type locality) or south (near Tabatinga but across the Solimões). Finally, our analysis of pelage coloration revealed three phenotypic forms: (1) south of the Rio Solimoes, (2) Eirunepé-Acre, upper Juruá basin; and (3) Japurá. More samples from both sides of Solimões in the region of Tabatinga will be necessary to ascertain the exact type locality for Spix's pygmaea and to resolve the current uncertainties surrounding pygmy marmoset taxonomy.