Tabula Microcebus: A transcriptomic cell atlas of mouse lemur, an emerging primate model organism

Mouse lemurs are the smallest, fastest reproducing, and among the most abundant primates, and an emerging model organism for primate biology, behavior, health and conservation. Although much has been learned about their physiology and their Madagascar ecology and phylogeny, little is known about their cellular and molecular biology. Here we used droplet- and plate-based single cell RNA-sequencing to profile 226,000 cells from 27 mouse lemur organs and tissues opportunistically procured from four donors clinically and histologically characterized. Using computational cell clustering, integration, and expert cell annotation, we defined and biologically organized over 750 mouse lemur molecular cell types and their full gene expression profiles. These include cognates of most classical human cell types, including stem and progenitor cells, and the developmental programs for spermatogenesis, hematopoiesis, and other adult tissues. We also described dozens of previously unidentified or sparsely characterized cell types and subtypes. We globally compared cell type expression profiles to define the molecular relationships of cell types across the body, and explored primate cell type evolution by comparing mouse lemur cell profiles to those of the homologous cells in human and mouse. This revealed cell type specific patterns of primate cell specialization even within a single tissue compartment, as well as many cell types for which lemur provides a better human model than mouse. The atlas provides a cellular and molecular foundation for studying this primate model organism, and establishes a general approach for other emerging model organisms.

Microenvironment Impacts the Molecular Architecture and Interactivity of Resident Cells in Marmoset Brain

The microenvironments of the brain consist of specialized cell types that together influence physiological functions in health and pathological outcomes in disease. Despite apparent differences in the density of neurons and oligodendrocytes in various milieus, such as gray matter (GM) and white matter (WM), the extent of structural and functional heterogeneity of other resident cells remains unclear. We profiled RNA in ~500,000 nuclei from 19 tissue types across the central nervous system of the healthy adult common marmoset (Callithrix jacchus) and mapped 87 identified subclusters (including neurons, glia, and vasculature) spatially onto a 3D MRI atlas. We performed cross-species comparison, explored regulatory pathways, surveyed cellular determinants of neurological disorders, and modeled regional intercellular communication. We found spatially segregated microglia, oligodendrocyte lineage cells, and astrocytes in WM and GM. WM-glia are diverse, are enriched with genes involved in stimulus response and biomolecule modification, and interact with other resident cells more extensively than their GM counterparts. GM-glia preserve the expression of developmental morphogens into adulthood and share 6 differentially enriched transcription factors that restrict the transcriptome complexity. Our work in marmoset, an experimentally tractable animal model with >5 times more WM volume and complexity than mouse, identifies novel WM-glia subtypes and their contributions to different neurological disorders. A companion Callithrix jacchus Primate Cell Atlas (CjPCA) is available through an online portal https://cjpca.ninds.nih.gov to facilitate data exploration.

Structure of the Receptor Binding Domain of EnvP(b)1, an Endogenous Retroviral Envelope Protein Expressed in Human Tissues

ABSTRACT EnvP(b)1 is an endogenous retroviral envelope gene found in human and other primate genomes. We report EnvP(b)1 sequences in primate genomes consistent with an integration event between 40 and 71 million years ago. Using a highly specific polyclonal antiserum raised against the putative receptor binding domain (RBD) of human EnvP(b)1, we detected expression in human placenta, ovaries, and thymus. We found that EnvP(b)1 is proteolytically processed, and using cell-cell fusion assays in multiple primate cell lines, we demonstrated that extant EnvP(b)1 proteins from a variety of primate genomes are fusogenic. This work supports the idea that EnvP(b)1 is under purifying selection and its fusogenic activity has been maintained for over 40 million years. We determined the structure of the RBD of human EnvP(b)1, which defines structural similarities with extant leukemia viruses, despite little sequence conservation. This structure highlights a common scaffold from which novel receptor binding specificities likely evolved. The evolutionary plasticity of this domain may underlie the diversity of related Envs in circulating viruses. IMPORTANCE Organisms can access genetic and functional novelty by capturing viral elements within their genomes, where they can evolve to drive new cellular or organismal processes. We demonstrate that a retroviral envelope gene, EnvP(b)1, has been maintained and its fusion activity preserved for 40 to 71 million years. It is expressed as a protein in multiple healthy human tissues. We determined the structure of its inferred receptor binding domain and compared it with the same domain in modern viruses. We found a common conserved architecture that underlies the varied receptor binding activity of divergent Env genes. The modularity and versatility of this domain may underpin the evolutionary success of this clade of fusogens.

The Virology of Taterapox Virus In Vitro

Taterapox virus (TATV) is phylogenetically the closest related virus to variola—the etiological agent of smallpox. Despite the similarity, few studies have evaluated the virus. In vivo, TATV can infect several animals but produces an inapparent infection in wild-type mice; however, TATV does cause morbidity and mortality in some immunocompromised strains. We employed in vitro techniques to compare TATV to ectromelia (ECTV) and vaccinia (VACV) viruses. Both ECTV and TATV replicate efficiently in primate cell lines but TATV replicates poorly in murine cells lines. Furthermore, TATV induces cytopathic effects, but to a lesser extent than ECTV, and changes cytoskeletal networks differently than both ECTV and VACV. Bioinformatic studies revealed differences in several immunomodulator open reading frames that could contribute to the reduced virulence of TATV, which were supported by in vitro cytokine assays.

Isolation and Characterization of a Novel Gammaherpesvirus from a Microbat Cell Line

ABSTRACT Bats are of significant interest as reservoirs for zoonotic viral pathogens; however, tools to dissect bat-virus interactions are limited in availability. This study serendipitously identified, in an established bat cell line, a fully replication-competent gammaherpesvirus; determined the complete genome sequence of the virus; and generated a viral transcript map. This virus can replicate in select human and nonhuman primate cell lines. However, analyses of viral sequences support a bat origin for this virus; we therefore refer to the virus as bat gammaherpesvirus 8 (BGHV8). The viral genome contains unique open reading frames that likely encode modulators of bat innate and adaptive immune signaling pathways and expresses viral miRNAs. The virus and its gene products should provide a unique tool to dissect both bat and gammaherpesvirus biology. While employing deep sequencing and de novo assembly to characterize the mRNA transcript profile of a cell line derived from the microbat Myotis velifer incautus, we serendipitously identified mRNAs encoding proteins with a high level of identity to herpesviruses. A majority were closely related to proteins of equine herpesvirus 2 (EHV-2), a horse gammaherpesvirus. We demonstrated by electron microscopy the presence of herpesvirus-like particles in the microbat cells. Passage of supernatants from microbat cells to Vero cells resulted in syncytium formation, and expression of viral genes and amplification of viral DNA were demonstrated by quantitative PCR. Susceptibility of human cell lines to productive infection was also demonstrated. Next-generation sequencing and de novo assembly of the viral genome from supernatants from Vero cells yielded a single contig of approximately 130 kb with at least 77 open reading frames (ORFs), predicted microRNAs (miRNAs), and a gammaherpesvirus genomic organization. Phylogenic analysis of the envelope glycoprotein (gB) and DNA polymerase (POLD1) revealed similarity to multiple gammaherpesviruses, including those from as-yet-uncultured viruses of the Rhadinovirus genus that were obtained by deep sequencing of bat tissues. Moreover, the assembled genome revealed ORFs that share little or no homology to known ORFs in EHV-2 but are similar to accessory proteins of other gammaherpesviruses. Some also have striking homology to predicted Myotis bat proteins. Cumulatively, this study provides the first isolation and characterization of a replication-competent bat gammaherpesvirus. IMPORTANCE Bats are of significant interest as reservoirs for zoonotic viral pathogens; however, tools to dissect bat-virus interactions are limited in availability. This study serendipitously identified, in an established bat cell line, a fully replication-competent gammaherpesvirus; determined the complete genome sequence of the virus; and generated a viral transcript map. This virus can replicate in select human and nonhuman primate cell lines. However, analyses of viral sequences support a bat origin for this virus; we therefore refer to the virus as bat gammaherpesvirus 8 (BGHV8). The viral genome contains unique open reading frames that likely encode modulators of bat innate and adaptive immune signaling pathways and expresses viral miRNAs. The virus and its gene products should provide a unique tool to dissect both bat and gammaherpesvirus biology.

Detection of EBV, HBV, HCV, HIV-1, HTLV-I and -II, and SMRV in Human and Other Primate Cell Lines

The high prevalence of contaminated cell cultures suggests that viral contaminations might be distributed among cultures. We investigated more than 460 primate cell lines for Epstein-Barr (EBV), hepatitis B (HBV), hepatitis C (HCV), human immunodeficiency virus type 1 (HIV-1), human T-cell leukemia/lymphoma virus I and II (HTLV-I/-II), and squirrel monkey retrovirus (SMRV) infections for risk assessment. None of the cell lines were infected with HCV, HIV-1, or HTLV-I/-II. However, one cell line displayed reverse transcriptase activity. Thirty-nine cell lines harbored EBV DNA sequences. Studies on the lytic phase of EBV revealed that five cell lines produce EBV particles and six further cell lines produced EBV upon stimulation. One cell line contained an integrated HBV genome fragment but showed no virus production. Six cell lines were SMRV-infected. Newly established cell lines should be tested for EBV infections to detect B-lymphoblastoid cell lines (B-LCL). B-LCLs established with EBV from cell line B95-8 should be tested for SMRV infections.

Role of the Alpha/Beta Interferon Response in the Acquisition of Susceptibility to Poliovirus by Kidney Cells in Culture

ABSTRACT Replication of poliovirus (PV) is restricted to a few sites, including the brain and spinal cord. However, this neurotropism is not conserved in cultured cells. Monkey kidney cells become susceptible to PV infection after cultivation in vitro, and cell lines of monolayer cultures from almost any tissue of primates are susceptible to PV infection. These observations suggest that cellular changes during cultivation are required for acquisition of susceptibility. The molecular basis for the cellular changes during this process is not known. We investigated the relationship between PV susceptibility and interferon (IFN) response in primary cultured kidney and liver cells derived from transgenic mice expressing human PV receptor and in several primate cell lines. Both kidneys and liver in vivo showed rapid IFN response within 6 h postinfection. However, monkey and mouse kidney cells in culture and primate cell lines, which were susceptible to PV, did not show such rapid response or showed no response at all. On the other hand, primary cultured liver cells, which were partially resistant to infection, showed rapid IFN induction. The loss of IFN inducibility in kidney cells was associated with a decrease in expression of IFN-stimulated genes involved in IFN response. Mouse kidney cells pretreated with a small dose of IFN, in turn, restored IFN inducibility and resistance to PV. These results strongly suggest that the cells in culture acquire PV susceptibility during the process of cultivation by losing rapid IFN response that has been normally maintained in extraneural tissues in vivo.

Contributions of Matrix and Large Protein Genes of the Measles Virus Edmonston Strain to Growth in Cultured Cells as Revealed by Recombinant Viruses

ABSTRACT The Edmonston strain of measles virus (MV) was obtained by sequential passages of the original isolate in various cultured cells. Although attenuated in vivo, it grows efficiently in most primate cell lines. Previous studies have revealed that MV tropism cannot be solely explained by the use of CD150 and/or CD46 as a cellular receptor. In order to evaluate the contributions of individual genes of the Edmonston strain to growth in cultured cells, we generated a series of recombinant viruses in which part of the genome of the clinical isolate IC-B (which uses CD150 as a receptor) was replaced with the corresponding sequences of the Edmonston strain. The recombinant virus possessing the Edmonston hemagglutinin (H) gene (encoding the receptor-binding protein) grew as efficiently in Vero cells as the Edmonston strain. Those viruses having either the matrix (M) or large (L) protein gene from the Edmonston strain could also replicate well in Vero cells, although they entered them at low efficiencies. P64S and E89K substitutions were responsible for the ability of the M protein to make virus grow efficiently in Vero cells, while the first half of the Edmonston L gene was important for better replication. Despite efficient growth in Vero cells, the recombinant viruses with these mutations had growth disadvantage in CD150-positive lymphoid B95a cells. Thus, not only the H gene but also the M and L genes contribute to efficient replication of the Edmonston strain in some cultured cells.