Skip to main content
Advertisement

Main menu

  • Home
  • Articles
    • Newest Articles
    • Current Issue
    • Methods & Resources
    • Archive
    • Subjects
  • Collections
  • Submit
    • Submit a Manuscript
    • Author Guidelines
    • License, Copyright, Fee
    • FAQ
    • Why Submit
  • About
    • About Us
    • Editors & Staff
    • Board Members
    • Licensing and Reuse
    • Reviewer Guidelines
    • Privacy Policy
    • Advertise
    • Contact Us
    • LSA LLC
  • Alerts
  • Other Publications
    • EMBO Press
    • The EMBO Journal
    • EMBO reports
    • EMBO Molecular Medicine
    • Molecular Systems Biology
    • Rockefeller University Press
    • Journal of Cell Biology
    • Journal of Experimental Medicine
    • Journal of General Physiology
    • Cold Spring Harbor Laboratory Press
    • Genes & Development
    • Genome Research

User menu

  • My alerts

Search

  • Advanced search
Life Science Alliance
  • Other Publications
    • EMBO Press
    • The EMBO Journal
    • EMBO reports
    • EMBO Molecular Medicine
    • Molecular Systems Biology
    • Rockefeller University Press
    • Journal of Cell Biology
    • Journal of Experimental Medicine
    • Journal of General Physiology
    • Cold Spring Harbor Laboratory Press
    • Genes & Development
    • Genome Research
  • My alerts
Life Science Alliance

Advanced Search

  • Home
  • Articles
    • Newest Articles
    • Current Issue
    • Methods & Resources
    • Archive
    • Subjects
  • Collections
  • Submit
    • Submit a Manuscript
    • Author Guidelines
    • License, Copyright, Fee
    • FAQ
    • Why Submit
  • About
    • About Us
    • Editors & Staff
    • Board Members
    • Licensing and Reuse
    • Reviewer Guidelines
    • Privacy Policy
    • Advertise
    • Contact Us
    • LSA LLC
  • Alerts
  • Follow lsa Template on Twitter
Resource
Source Data
Transparent Process
Open Access

Conservation of cell-intrinsic immune responses in diverse nonhuman primate species

Jenna M Gaska, View ORCID ProfileLance Parsons, Metodi Balev, Ann Cirincione, Wei Wang, Robert E Schwartz, View ORCID ProfileAlexander Ploss  Correspondence email
Jenna M Gaska
1Lewis Thomas Laboratory, Department of Molecular Biology, Princeton University, Princeton, NJ, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Lance Parsons
2Carl Icahn Laboratory, Lewis-Sigler Institute for Integrative Genomics, Princeton University, Princeton, NJ, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Lance Parsons
Metodi Balev
1Lewis Thomas Laboratory, Department of Molecular Biology, Princeton University, Princeton, NJ, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Ann Cirincione
1Lewis Thomas Laboratory, Department of Molecular Biology, Princeton University, Princeton, NJ, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Wei Wang
2Carl Icahn Laboratory, Lewis-Sigler Institute for Integrative Genomics, Princeton University, Princeton, NJ, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Robert E Schwartz
3Weill Cornell Medical College, Belfer Research Building, New York, NY, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Alexander Ploss
1Lewis Thomas Laboratory, Department of Molecular Biology, Princeton University, Princeton, NJ, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Alexander Ploss
  • For correspondence: aploss@princeton.edu
Published 24 October 2019. DOI: 10.26508/lsa.201900495
  • Article
  • Figures & Data
  • Info
  • Metrics
  • Reviewer Comments
  • PDF
Loading

Abstract

Differences in immune responses across species can contribute to the varying permissivity of species to the same viral pathogen. Understanding how our closest evolutionary relatives, nonhuman primates (NHPs), confront pathogens and how these responses have evolved over time could shed light on host range barriers, especially for zoonotic infections. Here, we analyzed cell-intrinsic immunity of primary cells from the broadest panel of NHP species interrogated to date, including humans, great apes, and Old and New World monkeys. Our analysis of their transcriptomes after poly(I:C) transfection revealed conservation in the functional consequences of their response. In mapping reads to either the human or the species-specific genomes, we observed that with the current state of NHP annotations, the percent of reads assigned to a genetic feature was largely similar regardless of the method. Together, these data provide a baseline for the cell-intrinsic responses elicited by a potent immune stimulus across multiple NHP donors, including endangered species, and serve as a resource for refining and furthering the existing annotations of NHP genomes.

Introduction

Many of the pathogens with a significant impact on human morbidity and mortality—such as HIV, hepatitis C virus, hepatitis B virus, and yellow fever virus—have a narrow host range limited to humans and select nonhuman primate (NHP) species (1). Broadly speaking, the host tropism of such viral pathogens can be determined by (i) the absence or incompatibility of (a) host factor(s) needed for part(s) of the viral life cycle; (ii) the presence of dominant restriction factors; and/or (iii) differences in immune responses. Understanding the molecular basis of host tropism can inform our understanding of intra- and interspecies transmission and aid in the generation of improved animal models and clinical therapies. Here, we analyzed how cell-intrinsic immune responses compare across humans and NHP species. As our closest evolutionary relatives, the tactics by which NHPs have evolved to confront pathogens and the corresponding ways those pathogens modulate host immune defenses—as highlighted by hotspots for positive selection in genes involved with immune/defense responses in humans and select NHPs ((2, 3, 4, 5), reviewed in reference 6)—could shed light on host range barriers pertinent to studying zoonotic infections and using NHPs as biomedical research models.

In the present study, we purposefully chose a cell type that could be easily collected from a wide array of multiple NHP species, was conducive to studies of cell-intrinsic immune responses and could be used for future generation of induced pluripotent stem cells to study host–pathogen interactions in desired cell lineages. The latter criterion would be especially advantageous for studying pathogens with both highly limited host and tissue tropism, such as hepatitis C virus, for which acquiring hepatocytes from NHPs can be extremely challenging for both ethical and financial reasons. Dermal fibroblasts (DFs) met all these requirements and could be obtained through biopsies or existing repositories. As an added advantage, we used only primary DFs which do not have the cell cycle dysregulation or disruption of immune signaling as observed for immortalized cells such as human hepatoma cell lines (7, 8).

We selected three donors from eight different NHP species covering ∼35 million years of evolution (Fig 1A and Table 1): great apes (Pan troglodytes, Pan paniscus, Gorilla gorilla, and Pongo pygmaeus), because of their close genetic relationship to humans but whose endangered status precludes their use in most scientific studies; three Old World monkey species (Papio anubis, Macaca nemestrina, and Macaca mulatta); and one New World species (Saimiri sciureus), which are all commonly used in biomedical research. As a more distant point in evolutionary time (ca. 65 million years since divergence), we also included mouse. For our initial studies comparing cell-intrinsic immunity in these DFs from diverse species, we aimed to provoke an immune response independent of a given species’ susceptibility and/or permissivity to a virus. Thus, we used polyinosinic:polycytidylic acid (poly(I:C)), a synthetic dsRNA analog that can be transfected into cells and potently stimulates an IFN-mediated response by resembling a pathogen-associated molecular pattern characteristic of RNA virus infections (Fig 1B) (9, 10, 11). This allowed us to compare the repertoire of cell-intrinsic responses available to DFs from each of these species unhindered by pathogen-mediated modulation that can occur as early as viral binding and entry.

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1. Comparing the cell-intrinsic immune response to poly(I:C) across multiple species.

(A) Schematic representation of the phylogenetic relationship of species included in this study. (B) Overview of the experimental workflow. (C, D, E) Using RT-qPCR, IFNβ (C), MX1 (D), and OASL (E) mRNA expression were assessed in primate and mouse primary DFs relative to the housekeeping gene HPRT1 at 12 and/or 24 h after poly(I:C) (∼53 ng/cm2) or mock transfection. Values shown are the fold changes relative to mock-transfected cells. Each data point represents an individual well of cells and is colored by donor (see Fig S2 for the color code used). Circles represent data from 24-well experiments, whereas squares are from the transfections performed in six-well plates that underwent RNA-Seq. Data information: In (C, D, E), bars depict the mean with SD. Note that the lower error bar for the orangutan samples at both time points and the pig-tailed macaque at 24 h are missing because they approach a value of 0. Where P ≤ 0.05, the P value is indicated. In (C), an unpaired t test was used to compare timepoints, with Welch’s correction performed where the SD between time points differed by more than a factor of 2. In (D, E), an ordinary one-way ANOVA was performed followed by Dunnett’s multiple comparisons test with a single pooled variance and the human samples serving as the control against which all comparisons were made. *P ≤ 0.05; **P ≤ 0.01; ***P ≤ 0.001; ****P ≤ 0.0001.

Source data are available for this figure.

Source Data for Figure 1[LSA-2019-00495_SdataF1.pdf]

View this table:
  • View inline
  • View popup
Table 1.

Known information for each DF donor.

We subsequently analyzed the transcriptomic profile of these cells, generating data in response to a mimic of viral infection in the broadest panel of NHP species compiled to date. Although the evolution of select factors in response to host–pathogen conflict has been explored (12, 13), there has not been a comprehensive analysis of cell-intrinsic responses in nonimmune cells across a diverse range of NHPs. Much of comparative primate immunology focuses on the major cellular constituents of innate and adaptive immunity (14, 15, 16, 17, 18, 19) and highly variable genetic loci such as those encoding MHC class I molecules (20, 21, 22), immunoglobulins (23), and MHC class I–related (MR1) molecules (24). Furthermore, such studies often use a limited number of NHP species and focus on differences at the genetic level versus expression. To our knowledge, the closest study to ours examining the functional consequences of differences and similarities in innate immunity in NHPs stimulated primary monocytes from chimpanzee, rhesus, and human with LPS (25). Two more recent studies also provided insights into cell-intrinsic immunity but on a much larger evolutionary scale, including human or human and rhesus macaque as the only primates (26, 27).

The immense transcriptomic dataset we have acquired and analyzed in the course of this work fulfills a need for not only better understanding cell-intrinsic immunity in our closest evolutionary relatives but also providing additional sequencing information from NHP species, including those that are endangered. In this study, we mapped the data from each species both to the human genome and to the species-specific genomes, creating the added opportunity of comparing these two approaches in parallel. Overall, we observed a high level of conservation in the functional responses of the NHP DFs to poly(I:C) and found that with the current state of genome annotations, the percent of reads assigned to a genetic feature were largely similar between the two mapping methods. We anticipate that making these data available will enhance the continued efforts to more fully annotate NHP genomes and to aid in the identification of novel transcripts unobserved in humans. In addition, we believe these unique data will greatly facilitate evolutionary analyses, especially in the area of comparative immunology.

Results and Discussion

Expression of select interferon-stimulated genes (ISGs) in response to poly(I:C) varies across diverse species

Triggering innate immune sensing pathways characteristically provokes the induction of hundreds of ISGs, which collectively create an antiviral state. However, it remains unclear how similar such responses are across NHP species. To initially test whether poly(I:C) would elicit a cell-intrinsic response in the DFs we had from 10 different species (each represented by three donors), we first assessed transcription via RT-qPCR of IFNβ at 12 and 24 h post-transfection and two well-characterized ISGs, myxoma resistance protein 1 (MX1) and 2′-5′-oligoadenylate synthetase-like (OASL, Oasl1 in mice), 24 h post-transfection using species-specific primers (Figs 1C–E, S1, and S2). IFNβ mRNA had on average less than a 10-fold change in expression relative to mock-transfected cells for most species (Fig 1C). However, bonobo, gorilla, orangutan, and squirrel monkey did exceed this for at least one of the time points tested. Bonobo and gorilla had significantly higher IFNβ expression at 12 h compared with 24 h, whereas the opposite trend was observed for rhesus macaque. The two ISGs examined demonstrated far greater fold changes relative to mock-transfected cells, underscoring the prominence of the downstream response to IFN at 24 h post-transfection. All species, with the exception of mouse, had at least an average 100-fold increase relative to mock in MX1 mRNA, although there was marked variation amongst rhesus donors. Average levels of OASL mRNA were also elevated across the primate donors, albeit with greater donor variation, but all three mouse donors still displayed minimal change. The donor–donor variation we observed was not surprising as the NHP fibroblasts were acquired from outbred donors, with at least one male and one female represented for each species. As these data represent only a small sample of the hundreds of ISGs (28, 29, 30) whose expression could be changing after poly(I:C) transfection, we performed RNA-Seq on total mRNA isolated from these DFs 24 h post-transfection for a more comprehensive view. Unlike recent studies that looked earlier at 4 h post-immune stimulation (26, 27), we wanted to analyze effects further downstream as the cell-intrinsic response was amplified over time.

Figure S1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S1. Extended kinetics of IFNβ expression in select donors.

For select donors, time points were collected in addition to those at 12 and 24 h post-transfection with poly(I:C) (∼53 ng/cm2) or mock transfection as shown in Fig 1. Each data point represents an individual well of cells and is colored by donor (Fig S2). Circles represent data from 24-well experiments whereas squares are from the transfections performed in six-well plates that underwent RNA-Seq. Bars depict the mean with SD.

Figure S2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S2. Color coding system used throughout the article for the individual donors.

Using species-specific genomes increases read alignment but to regions with no assigned features

In working with such a diverse array of species, one challenge was choosing the genomes upon which to map the NHP-derived RNA-Seq reads. Since the genomes for human and mouse are well-annotated and species-specific resources available for downstream analyses, we felt confident in aligning reads from these species to their respective genomes. However, for the NHP species, to make our analysis more thorough, we used two approaches in parallel: mapping all primate-derived reads to the human genome as previously performed (31, 32, 33, 34) (hereby referred to as “human method”) or to their respective genomes as they currently exist on Ensembl (Table 2) (“species method”). Human and mouse reads averaged above 90% alignment and ∼80% assignment. For the great ape species examined, the percent alignment was generally the same regardless of the genome used, exceeding 85% (Fig S3). For the Old World monkeys, there was a small but noticeable decline when mapping to the human genome, with percent alignment close to 85% for the macaques and dipping under 80% for samples from olive baboon. Aligning reads from each of these species to their respective genomes increased these values, especially for rhesus and pig-tailed macaque, which approached 90% and 95%, respectively. Moving further out evolutionarily, squirrel monkey had the lowest percent alignment to the human genome, but still most reads did map, centering around 70%. Not surprisingly, more than 90% of reads were aligned when the squirrel monkey genome was used.

View this table:
  • View inline
  • View popup
Table 2.

Reference genomes used.

Figure S3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S3. The percentage of RNA-Seq reads aligned and assigned for each species after mapping to either the human genome or the species-specific genome.

Total RNA was extracted from DFs 24 h after poly(I:C) or mock transfection followed by cDNA library preparation and sequencing. Samples originated from nine species, as shown grouped here, with three donors per species and each donor transfected in triplicate with either poly(I:C) or mock transfected. The donor ID and transfection condition (“mock” or “treated” with poly(I:C)) are indicated on the X axis of each grouped species plot. The NHP species reads were mapped to the human genome (“human method”) or their species-specific genome (“species method”) (see Table 2) and then the percent alignment and assignment to each genome determined by RNA STAR and featureCounts, respectively. As indicated in the figure legend, bars in shades of red indicate percent alignment and blue percent assignment. The darker shade of each is for the mapping performed with the human genome and the lighter shade the mapping performed with the species-specific genome for each NHP species.

Source data are available for this figure.

Source Data for Figure S3[LSA-2019-00495_SdataFS3.csv]

Despite these variations in alignment depending on the reference genome used, the percent of reads assigned to a genome feature by the “human method” tended to be higher (up by ∼5–6%) or comparable with that of the “species method.” Only for some chimpanzee, pig-tailed macaque and squirrel monkey samples did the species mapping improve the percent assignment (once more up by ∼5–6%). Orangutan and olive baboon demonstrated the greatest differences, with the “species method” as much as ∼10 or 15 percentage points lower, respectively, than the “human method.”

Taking a closer look at the relationship between the percent of reads assigned to a genome feature and the mapping method, we observed that the benefit of the “species method” over the “human method,” whereby the percentage of reads failing to be assigned a feature because of poor mapping quality declined, was less advantageous than expected as an increased percentage of the reads that mapped were to regions with no assigned feature (Fig S4). Notably, this was consistent across NHP species, irrespective of evolutionary relationship to humans. Thus, the percentage of reads assigned a feature was generally comparable between the two mapping methods or better by the “human” method. The NHP genome annotations as they currently exist resulted in on average anywhere from ∼12–22% of reads mapping to sequence with no assigned feature (for comparison, <5% of human and mouse reads fell in this category after alignment to their respective genomes). In contrast, less than 5.5% of reads fell in this category when using the “human method.” Identifying what these unassigned features might be beyond the extensive homologous sequence searching across existing genomes already performed by Ensembl is outside the scope of this article. These data highlight the weaknesses of the current NHP genome annotations and the core difficulty in performing comparative analyses, as the more divergent, species-specific sequences will still be missed regardless of the genome used. We hope that sequencing datasets like ours from multiple donors of multiple NHP species will aid in improving such annotations.

Figure S4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S4. Comparison of feature assignment between the human and species mapping methods.

Total RNA was extracted from DFs 24 h after poly(I:C) or mock transfection followed by cDNA library preparation and sequencing. Each grouped bar graph represents an NHP species, with results from reads aligned to the human genome shown in the left column and those from reads aligned to the species-specific genome (see Table 2) in the right column. The donor ID and transfection condition (“mock” or “treated” with poly(I:C)) are indicated on the X axis of each plot. featureCounts was used to determine the percent of reads that were assigned to an existing genome feature (red), the percent that mapped but with no known feature to assign the reads to (orange), and the percent that were unassigned because of low mapping quality (blue).

Source data are available for this figure.

Source Data for Figure S4[LSA-2019-00495_SdataFS4.csv]

Species cluster by evolutionary relationship when looking at overall transcriptomic response

Across all NHP species examined, most genes had a one-to-one human ortholog, with the second largest category of genes those without a known ortholog (Fig S5). For both the “human” and the “species” methods, we first limited the resultant mapped read counts to only those genes that had one-to-one orthologs across all species (“common orthologs”), resulting in a common denominator of 11,677 genes (Supplemental Data 1). Using such an approach, we were able to directly compare all eight NHP species with one another, finding that by principal component analysis (PCA), the samples clustered according to the evolutionary relationship of the species, with great apes, Old World monkeys, and the single New World monkey species forming their own distinct groupings regardless of whether they were mapped by the “human” or the “species” method (Fig 2A and B). For both methods, the orangutan samples formed a group slightly removed from the other great ape species. However, only by the “species” method did we observe resolution of the Old World monkey species into three distinct groups unlike by the “human” method. This could be, at least in part, attributed to the larger differences in percent assignment of reads, especially between olive baboon and the two macaques, when using the “species” method (Fig S3).

Figure S5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S5. Types of human orthologs found across NHP species.

The ENSEMBL gene IDs for each NHP species under examination were gathered and the human ortholog, if any, retrieved from bioMart along with the ortholog “type” (none, one-to-one, one-to-many, or many-to-many). The number of NHP ENSEMBL gene IDs falling into each of these ortholog “type” categories are shown here in a stacked bar graph, with each column representing a different species and the colors within each column representing a different ortholog type. chimp, chimpanzee; orang, orangutan; rhesus, rhesus macaque; ptmac, pig-tailed macaque; bab, olive baboon; and sqmonk, squirrel monkey.

Supplemental Data 1.

List of human Ensembl IDs for genes which have a one-to-one ortholog across all NHP species (“common orthologs”).[LSA-2019-00495_Supplemental_Data_1.xls]

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2. Overall transcriptomic responses to poly(I:C) cluster by species’ phylogenetic relationships and innate immune genes by treatment condition.

(A, B, C, D) RNA sequencing reads were aligned to the human genome (A, C) or the respective NHP species from which the samples were derived (B, D). The resultant counts from either mapping method were then limited to genes with a one-to-one human ortholog (based off Ensembl ID) across all species, resulting in ∼11,000 genes (referenced as “common orthologs” in the main text). For counts mapped to the NHP genomes, the human ortholog Ensembl ID was used for all subsequent analysis for ease of comparison. These counts were normalized using DESeq2 default options and finally transformed using the regularized log2 function in DESeq2. (A, B, C, D) The first two components of the PCA are depicted here for either all ∼11,000 genes (“common orthologs”) (A, B) or for a subset of these genes which overlapped with those found in the database InnateDB (“innate genes”) (C, D). Each point represents an individual donor and is colored by species as shown in the figure legend, with shades of blue indicating great ape species, shades of green Old World Monkeys, and purple the single New World Monkey used in this study. Circles indicate mock-transfected samples and triangles poly(I:C)-transfected samples. (E, F) The counts from mapping to either the human or the species-specific genome were limited on a species-by-species basis to the Ensembl IDs that had a one-to-one human ortholog (referenced as “species orthologs” in the main text). The counts of these genes were then processed in DESeq2 to determine DGE, specifically comparing the DGE profile of each NHP species (poly(I:C)-transfected versus mock-transfected) to that of human. Thus, genes from NHP species that significantly differed from human in their response to poly(I:C) could be determined and are shown here as pseudo-violin plots from either the human genome (E) or the species-specific genome (F) mapping method. Each point represents an individual gene and is colored by species as shown in the figure legend. Data information: in (E, F), the differentially expressed genes shown have Padj ≤ 0.05.

Innate immune genes strongly up-regulated by poly(I:C) transfection across species

As we had transfected the cells with poly(I:C), a synthetic dsRNA mimic known to induce cell-intrinsic immune responses, we wanted to focus further on known innate immunity genes and how their expression compared across these primate species. We limited our “common orthologs” output to a manually curated set of genes involved with innate immunity, InnateDB (35). This database includes 988 distinct gene symbols for human, 656 of which fell into our “common orthologs” set as one-to-one orthologs across all species (Supplemental Data 2). Examining these genes across our samples by PCA, we now observed clear separation of samples based on whether they had been transfected with poly(I:C) (Fig 2C and D, mouse shown in Fig S6), indicating the strong contribution of treatment in differentiating samples based off these genes. Within the treated and mock samples, the great ape samples did cluster separate from the monkey species, although mapping by the “species method” did result in greater separation between Old and New World monkey species. The squirrel monkey samples, our only New World monkey species under consideration, were closer than the Old World monkeys to the great ape cluster, especially after the “species method,” suggesting higher similarity of these genes’ expression with the great apes even though they are more distantly related.

Supplemental Data 2.

The 656 genes from the InnateDB database that were among the “common orthologs” listed in Supplemental Data 1.[LSA-2019-00495_Supplemental_Data_2.xls]

Figure S6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S6. PCA plots of mouse samples.

Reads from the murine samples were aligned to the mouse genome and the species annotations used for all analysis. These counts were normalized using DESeq2 default options and finally transformed using the regularized log2 function in DESeq2. The first two components of the PCA are depicted here for either all genes (left) (not limited by any type of orthologs with human) or for a subset of innate immune response genes from the database InnateDB (“innate genes”) (right). InnateDB has both human and murine genes listed, so whereas the primate samples used the human genes, here we were able to use the murine genes. Each point represents an individual donor, colored as shown in the legend. Circles indicate mock-transfected samples and triangles poly(I:C)-transfected samples.

We examined the sample-to-sample distances of these same data to see if we still observed a closer relationship between the great apes and squirrel monkey when not limited to just the first two principal components. Hierarchical clustering of the sample-to-sample distances was highly similar regardless of mapping method, clearly distinguishing between mock and poly(I:C)-transfected samples as well as the distinct nature of poly(I:C)-transfected samples from gorilla donor PR00107 (Fig S7). In addition, despite the PCA plot for the “human method” showing less separation between the New and Old World monkeys, we were able to confirm the greater similarity of the squirrel monkey samples to the great apes in regards to this gene set (Fig S7). Although squirrel monkey was the most distant of our NHP species from humans, the clustering of this species with the great apes, instead of out beyond the Old World monkeys, suggests strong conservation of the response in these innate immune genes after poly(I:C) transfection. However, the exact basis for these observations cannot be determined, especially without examining more New World monkey species to see if this is broadly applicable.

Figure S7.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S7. Cluster dendrogram of samples based off read counts for InnateDB genes with a one-to-one human ortholog across all NHP species examined.

Hierarchical clustering was determined for the sample–sample distances for the normalized, regularized log2-transformed counts that were used to generate Fig 2 PCA plots of InnateDB genes with a one-to-one human ortholog across all NHP species. (A, B) The human method results are shown in (A) and the species method results in (B). Each tip label indicates the donor ID, treatment condition and replicate number. The tip labels are colored according to the larger classifications of the species (i.e., great ape or Old World Monkey) as used throughout the article. The tips themselves are colored based on treatment condition (orange for mock-transfection and red for poly(I:C)-transfection).

Determining the NHP response to poly(I:C) distinctive from that of human

Although informative for broader analyses, using the genes with a one-to-one human ortholog across all eight species was far more limiting than if we did so on a species-by-species basis (referred to as “species orthologs”). Such an approach substantially increased the possible number of Ensembl IDs from 11,766 to 15,793-22,040 (Table 3; note that after determining the one-to-one human ortholog for a given NHP gene, the human Ensembl ID was used). In looking at the differential gene expression (DGE) profiles (poly(I:C)-transfected versus mock-transfected) for each species using “species orthologs,” we observed similar distributions, with more genes increasing versus decreasing in expression after poly(I:C) transfection compared with mock-transfected cells (Fig S8, and Supplemental Data 3 and 4). Although the range of expression for up-regulated genes was larger compared with down-regulated genes, most genes had a |log2(fold change)| <3 relative to mock-transfected cells. Because we were most interested in how the NHP species compared with human, we used the latter as our baseline to assess DGE for each species (Padj ≤ 0.05), finding on average ∼2,800 genes differentially expressed compared with human among the great ape species and ∼4,800 among the monkeys (Fig 2E and F, and Supplemental Data 5).

View this table:
  • View inline
  • View popup
Table 3.

Number of Ensembl ID annotations for each NHP species that is listed as having a one-to-one human ortholog.

Figure S8.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S8. Summary of DGE in diverse species after poly(I:C) transfection.

poly-A–containing RNA transcripts isolated from DFs mock-transfected or transfected 24 h with poly(I:C) were analyzed by RNA-Seq. (A, B) Shown are pseudo-violin plots of the distribution of differentially expressed genes (depicted as the log2(foldchange) in expression between poly(I:C)-transfected samples and mock samples) in human and NHP samples after mapping RNA-Seq reads to either the human (A) or species-specific (B) genomes. In both cases, genes were limited on a species-specific basis to those that have a one-to-one human ortholog. Hence, the top panel of (A) and (B) show the expression of genes from the NHP species and the lower panel that of the “related” human DGEs, which were limited to the same genes as those for the corresponding NHP species. Each pseudo-violin plot is colored according to species as shown in the figure legend. (C) Pseudo-violin plot of the DGE for the RNA-Seq reads from mice, which were all aligned to the mouse genome. Data information: in (A, B, C), the differentially expressed genes shown have Padj ≤ 0.05.

Supplemental Data 3.

These files are the DGE (poly(I:C)–transfected versus mock-transfected) outputs for each of the species after mapping to the human genome and limiting the genes for each NHP species to those that have a one-to-one human ortholog. The Ensembl IDs and gene information used are that of human. For each NHP species, there is an accompanying human DGE profile corresponding to the same limited list of one-to-one orthologs. Please see the README file for this dataset for further details.[LSA-2019-00495_Supplemental_Data_3.zip]

Supplemental Data 4.

These files are the DGE (poly(I:C)–transfected versus mock-transfected) outputs for each of the species after mapping to the species-specific genome and limiting the genes for each NHP species to those that have a one-to-one human ortholog. The Ensembl IDs and gene information used are that of human. For each NHP species, there is an accompanying human DGE profile corresponding to the same limited list of one-to-one orthologs. Please see the README file for this dataset for further details.[LSA-2019-00495_Supplemental_Data_4.zip]

Supplemental Data 5.

These files represent the DGE of each NHP species relative to the human DGE profile in response to poly(I:C) treatment. As mentioned in the main text, the genes for each NHP species were limited to those that have a one-to-one human ortholog. The Ensembl IDs and gene information used are that of human. Please see the README file for this dataset for further details.[LSA-2019-00495_Supplemental_Data_5.zip]

To find genes in each species with markedly different expression in response to poly(I:C) from that of human, we limited these DGEs to those with a |log2(fold change)| ≥ 3 and then compared their occurrence across the species. By either mapping method, more genes overall were differentially expressed compared with human in accordance with the evolutionary distance of the species from human (Fig 3, horizontal bar graphs). Furthermore, ∼80–89% of these genes overlapped between the two mapping methods (Table 4). A fraction of these genes were unique to each NHP species (Fig 3, first eight columns), with proportionally the greatest in two of the species more distantly related to human: squirrel monkey and olive baboon. By either mapping method, most genes unique to a given species had higher expression relative to human except for rhesus macaque, where the opposite trend was observed, and squirrel monkey and pig-tailed macaque, which each had an almost equal number of genes with higher and lower expression (Fig 3, pseudo-violin plots). In all cases, most genes found as uniquely different from human in a given species or species grouping were one-to-one orthologs across all species (Fig S9). Thus, these genes’ status as “unique” was not simply a result of their exclusion by our “species orthologs.”

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3. Overview of genes in NHP species and subgroups that differ markedly in their response to poly(I:C) compared with that of human.

(A, B) For the genes shown in Fig 2E and F, the RNA-Seq reads were mapped to either the human genome (A) or the species-specific genome (B). The resultant gene counts from each mapping method were then limited on a species-by-species basis to those for which there is a one-to-one human ortholog. DGE upon poly(I:C) transfection for each of the NHP species was then assessed relative to that of human. The vertical bar graph in the upper panel displays the number of genes which were differentially expressed (Padj ≤ 0.05, |log2(fold change)| ≥ 3) by a given species or group of species, as indicated by the matrix below the graph, compared with human. Each bar, thus, represents a unique group of genes with the matrix beneath indicating in which species these genes are expressed, similar to what is represented by the “intersection” of circles in a Venn diagram (as indicated by the inset Venn diagram in (A)). For ease of presentation, only the intersections with 14 or more genes are shown, except for the intersection of seven genes that differ from human in their expression upon treatment in all the species. The points in the matrix correspond to the different species (color-coded) from which DFs were sourced, with shades of blue indicating great ape species, shades of green Old World Monkeys and purple the single New World Monkey used in this study. The inset horizontal bar graph shows the total number of differentially expressed genes for each species based on the cutoff used (Padj ≤ 0.05, |log2(fold change)| ≥ 3). The pseudo-violin plot below the matrix summarizes the differential log2(fold change) expression relative to human for the genes included in the corresponding vertical bar. For genes expressed by multiple species (i.e., all except the first eight columns), the average log2(fold change) across species for each gene was used. For the first eight columns of genes unique to each species, the dots are colored in accordance with the horizontal bar graph. An orange asterisk with corresponding orange points in the pseudo-violin plot highlights the gene groups of particular interest that are in common across larger NHP groupings: all Old World monkeys, all monkeys, or all species.

View this table:
  • View inline
  • View popup
Table 4.

Overlap of genes between mapping methods for each NHP species significantly different in their transcriptomic response to poly(I:C) compared with that of human.

Figure S9.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S9. Most genes that are “uniquely” (i.e., only observed in one species or in a certain subset of species) significantly differentially expressed (Padj ≥ 0.05) compared with human in response to poly(I:C) have a one-to-one human ortholog in the other NHP species.

These data correspond with that shown in Fig 3 and delineates whether the genes that were found to be “unique” in their differential expression after poly(I:C) treatment compared with human was due to those genes having a one-to-one human ortholog only in that particular species or species group. Because the determination of one-to-one orthologs was done on a species-by-species basis before the analysis summarized in Fig 3, some genes could be found as “unique” to a species just because there was no one-to-one human ortholog in any of the other species. The upper left panel shows the “unique” genes amongst all the monkey species, the upper right panel that for the Old World monkeys and the bottom two rows those that were “unique” to individual species. The columns are colored based on the species as done throughout the rest of the article and as shown in the figure legend. The first column for each grouped plot (or columns in the case of the upper two panels) is that of the species given in the plot title. The y axis shows the number of genes that have a one-to-one human ortholog for each species. The maximum value is that of the first column (or group of columns for the upper two panels) and corresponds to the number of genes contained within that intersection as shown in Fig 3.

In addition, of the possible combinations of species with common genes differentially expressed compared with human, we focused on NHP groupings of evolutionary significance (Fig 3, columns with orange asterisks): Old World monkeys, all monkey species, and all species (there were no genes common to all the great apes). For these genes, the average differential expression relative to human was determined for the species in each grouping and tended by either mapping method to be lower than that of human (Fig 3, pseudo-violin plots). Our finding that of the 3,000–6,000+ genes differentially expressed (Padj ≤ 0.05) upon poly(I:C) treatment in the NHPs <500 were differentially expressed compared with human at our cutoff highlights the similarity of these species with human.

For all these genes differentially expressed compared with human in individual NHP species or groups of species, we then examined the DGE profiles we had generated of poly(I:C) versus mock-transfected cells for each species (Supplemental Data 3 and 4). Thus, we aimed to find the basis for each gene’s significantly different expression compared with human. For example, genes with a negative log2(fold change) in comparing humans versus NHPs could be because the gene’s expression was unchanged by poly(I:C) transfection in NHPs but increased in humans, was decreased upon poly(I:C)-transfection in NHPs but less so or not at all in human, or did increase in NHPs after poly(I:C) transfection but at a magnitude lower than that for human. We prepared heat maps to answer this question for the genes “unique” to a given species (Fig S10) and then the three groups of interest (Old World monkeys, all monkeys, and all species) after mapping by the “human method” (Fig 4A) or the “species method” (Fig 4B). For these three groupings, regardless of mapping method, most genes that strongly increased upon poly(I:C) transfection in humans showed weak changes that were often nonsignificant in the NHP species for that group (Supplemental Data 6). With a few exceptions, when a gene’s expression was significantly different between treatment conditions in the NHP species, the expression of the NHP genes tended to follow the same trend as humans but at a weaker magnitude. In the smaller group of genes different across all species compared with human, all the genes except one demonstrated decreased or negligible change in expression after poly(I:C) transfection in the NHP species but increased expression in humans. The exception, an InnateDB gene phosphoinositide-3-kinase adaptor protein 1 (PIK3AP1; plays a role in immune cell development and controlling cytokine production), stood out as being up-regulated, albeit at a lower magnitude than in humans, in squirrel monkey, gorilla, and bonobo but negligibly changed in the remaining species. Between mapping methods for these smaller subgroups, most genes overlapped.

Figure S10.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S10. Differential expression profiles on a species-by-species basis for genes significantly different from human and unique to individual NHP species after poly(I:C) transfection.

(A, B) For the genes found to be significantly differentially expressed compared with human (Padj ≤ 0.05, absolute value [log2(foldchange)] ≥ 3) among individual NHP species, the DGE between poly(I:C)-transfected and mock-transfected samples was determined for each NHP species and is shown with the accompanying expression profile for the same genes in human. (A, B) Because read counts were mapped to either the human (A) or the species-specific genome (B) and the genes limited on a species-by-species basis to the ENSEMBL IDs that had a one-to-one human ortholog (referenced as “species orthologs” in the main text), each NHP species has a corresponding human DGE profile limited to those same genes. Thus, the human profiles were generated from the same set of biological samples just limited to different lists of genes for each species. For presentation purposes of the multi-species groups, the human values shown are from the olive baboon-ortholog-limited expression profile. Because these genes were found across multiple species, the expression pattern only marginally differs between the various accompanying human profiles. Genes that are decreased in expression upon poly(I:C) treatment relative to mock are shown in orange and those increased in purple, with the intensity of the color corresponding to the magnitude of the change (expressed as log2(foldchange)). The names of genes which are from the InnateDB list are shown in bold italics.

Figure 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4. Differential expression profiles on a species-by-species basis for genes significantly different from human in subgroups of NHP species after poly(I:C) transfection.

(A, B) For the genes found to be significantly differentially expressed compared with humans (Padj ≤ 0.05, |log2(fold change)| ≥ 3) among all Old World monkeys, all monkeys, or all species (as shown in Fig 3), the DGE between poly(I:C)-transfected and mock-transfected samples was determined for each NHP species and is shown with the accompanying expression profile for the same genes in human. Because read counts were mapped to either the human (A) or the species-specific genome (B) and the genes limited on a species-by-species basis to the Ensembl IDs that had a one-to-one human ortholog (referenced as “species orthologs” in the main text), each NHP species has a corresponding human DGE profile limited to those same genes. Thus, the human profiles were generated from the same set of biological samples just limited to different lists of genes for each species. For presentation purposes of the multispecies groups, the human values shown are from the olive baboon ortholog–limited expression profile. Because these genes were found across multiple species, the expression pattern only marginally differs between the various accompanying human profiles. Genes that are decreased in expression upon poly(I:C) treatment relative to mock are shown in orange and those increased in purple, with the intensity of the color corresponding to the magnitude of the change (expressed as log2(fold change)). The names of genes which are from the InnateDB list are shown in bold italics.

Supplemental Data 6.

These files give the log2(fold change) values (poly(I:C)–transfected versus mock-transfected) for the genes either uniquely differentially expressed compared with humans upon poly(I:C) treatment in a single NHP species or in a group of NHP species (i.e., Old World Monkeys). Please see the README file for this dataset for further details.[LSA-2019-00495_Supplemental_Data_6.zip]

For the genes “unique” to each of the great ape species, most were significantly different upon poly(I:C) transfection, whereas for these same genes in humans, the change was negligible (Fig S10 and Supplemental Data 6). The number of nondifferentially expressed genes was more evenly split between human and each of the macaque species (Fig S10 and Supplemental Data 6). In accordance with the summary data in Fig 3, most genes in the rhesus samples had lower log2(fold change) values compared with humans, with up-regulation occurring in humans but nonsignificant or low-level increases in rhesus. The higher number of genes for squirrel monkey and olive baboon is likely reflective of their more distant relationship to human (Fig S10 and Supplemental Data 6). For both of these NHP species, there was a subgroup of genes whose expression was clearly up- or down-regulated after poly(I:C) treatment that were nonsignificantly or negligibly different in the human samples. Likewise, the opposite trend was seen at either extreme of the human genes, with squirrel monkey and, to a lesser extent olive baboon, showing smaller or nonsignificant changes in genes that were strongly up- or down-regulated in human.

Enrichment for cell-intrinsic immune response pathways conserved across all species

Biological outcomes are not necessarily dictated by individual genes but rather the concerted efforts of multiple genes acting in tandem. Although innate immune genes were present amongst genes differentially expressed in NHP but not in human samples, it was unclear if this handful of players was resulting in large-scale differences in the cell-intrinsic immune responses of these various species to poly(I:C). Furthermore, as only one time point post-transfection was tested in this study, it is impossible to determine if genes with altered expression from humans are differentially stimulated by poly(I:C) or the result of varying transcriptional kinetics. For a broader overview of these DGE profiles, we performed gene set analysis, whereby data are assessed for an enrichment of genes related in some way, such as their participation in the same pathway or biological function (36, 37, 38, 39). Because each species was limited to its own number of one-to-one orthologs, we prepared a human DGE profile for each species limited in the same way to ease comparisons. We used the R package Generally Applicable Gene-set Enrichment (GAGE) (40) to determine the Kyoto Encyclopedia of Genes and Genomes (KEGG) (41, 42, 43) pathways with the most significantly altered gene expression after poly(I:C) transfection. To reflect the nature of pathways in biological systems, bidirectional changes in gene expression were considered when determining pathway enrichment. For each NHP species besides rhesus macaque and regardless of the mapping method, genes were significantly enriched (qval ≤ 0.09) for multiple pathways known to be important for cell-intrinsic immune responses, either directly or indirectly, including Toll-like receptor (TLR) signaling, nucleotide-binding oligomerization domain (NOD)-like receptor signaling, cytosolic DNA sensing, JAK-STAT signaling, antigen processing and presentation, and chemokine signaling (Fig 5; the latter two were also significant for mice, Table 5). Although all the NHPs had pathway profiles that corresponded well with the accompanying human profile, rhesus macaque failed to meet the significance cutoff for any pathway. This was not due to limiting the list of genes to one-to-one orthologs as evidenced by the accompanying human profile, which was limited to the same gene list. The top-ranking pathways were the same as the other species but likely did not pass the significance threshold because of the lower number of differentially expressed genes relative to the other species. In some instances, such as NOD-like receptor signaling for olive baboon and retinoic acid-inducible gene I–like receptor signaling for bonobo, the significance value fell just outside of the cutoff, indicating that although perhaps not to the same extent as the other NHPs, these pathways are still important. Apoptosis and osteoclast differentiation genes stood out as clearly significant in their enrichment only for squirrel monkey and some of the corresponding human profiles. Although apoptosis is less unexpected, osteoclast differentiation could be appearing as a result of the common origin of osteoclasts and immune cells from hematopoietic cells and the role of TLR stimulation in disrupting osteoclast differentiation (44).

Figure 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 5. Enrichment of genes involved with cell-intrinsic immune pathways dominate the response to poly(I:C) across all species.

(A, B) Read counts mapped to the human (A) or species-specific genomes (B) were limited on a species-by-species basis to the Ensembl IDs that had a one-to-one human ortholog (referenced as “species orthologs” in the main text). The resultant genes were then analyzed using GAGE, with the only genes excluded being those that lacked an Entrez ID (the IDs used by GAGE) and/or had a missing Padj value. A heat map is shown of pathways (the row labels) with a q value ≤ 0.09 for at least one species (the column labels), with the cells colored by −log10(q value) so that cells meeting the cutoff are shades of purple and those that do not are shades of orange, with white set at the cutoff of −log10(0.09). Because genes were limited on a species-by-species basis to those that had a one-to-one human ortholog, each NHP species has a human profile that was limited to that same gene set (the eight right most columns of each heat map). The reads were derived from the same experimental samples—only the downstream analysis differs in terms of the genes that were included by their ortholog status.

Source data are available for this figure.

Source Data for Figure 5[LSA-2019-00495_SdataF5.xls]

View this table:
  • View inline
  • View popup
Table 5.

Top 12 pathways for which genes are enriched in mice after poly(I:C) transfection.

Together, these data mimicking a viral infection serve as a first step in comparing the cell-intrinsic immune responses of the most diverse panel of NHP species to date. In using a synthetic compound as a proxy for RNA virus infection, we were able to analyze transcriptomic responses without having to account for species-specific permissiveness and/or susceptibility to a pathogen. While there were of course differences in individual gene expression profiles after poly(I:C) transfection across species, these ultimately resulted in similar outcomes, activating pathways important to cell-intrinsic immunity. We compared our findings to a recent study of rhesus, mouse, and human primary DFs transfected with poly(I:C) for 4 h (26). Although our time point was 20 h later, some of these early genes (filtered for primates to those with a one-to-one human ortholog) remained up-regulated in these same species at 24 h (Padj ≤ 0.05, log2(foldchange) ≥ 3) − ∼22% for rhesus, ∼38% for mouse, and ∼60% for human (Supplemental Data 7). This comparison highlights the rapid nature of cell-intrinsic immune responses and the sustained transcription of these responses promoted by continued ISG production. However, it also suggests species–species differences in the kinetics of such responses that will be an important area for further research. Similarly, among the 62 genes that were commonly increased in response to type I IFN treatment in species far more divergent from the ones we studied (bats, chicken, horse, sheep, rat, cow, human, and pig) (27), from 50% (rhesus) to 77% (squirrel monkey) of the genes that had a one-to-one human ortholog in each of our NHP species (>90% of the 62 genes) were also up-regulated 24 h post-poly(I:C) treatment using the same cutoffs as those of Shaw et al (Padj ≤ 0.05, log2(foldchange) ≥ 2; Supplemental Data 8). In comparison, our human samples showed ∼80%.

Supplemental Data 7.

These files contain the DGE profiles from our experiments for the genes found to be differentially expressed in mouse, rhesus and human by Hagai et al. in their 2018 work (26).[LSA-2019-00495_Supplemental_Data_7.zip]

Supplemental Data 8.

These files contain the DGE profiles from our experiments for the genes identified as “Core ISGs” in the work of Shaw et al in 2017 (27).[LSA-2019-00495_Supplemental_Data_8.zip]

Although we observed a high level of similarity across these NHP species, this is not to say that in the context of a viral infection, all species would mount the same response. In transfecting the cells with poly(I:C), we provoked a response that is not tempered by the antagonism of viral proteins and thus closer to what the potential cell-intrinsic immune response of these cells would be if left uninhibited by the actions of viral proteins. As a result, the data presented here do not presume to predict the susceptibility of these species to particular viral pathogens or whether other cell types from these species would react the same way. We tested the assumption that the species’ response to poly(I:C) would be highly similar as this had not previously been formally demonstrated for all of these species. Thus, although striking differences between species would have been of great interest, a quantitative demonstration of their similarity is still very useful. Knowing the unfettered capabilities of these cells in terms of cell-intrinsic immunity provides a solid foundation to explore how this more maximal response is potentially dampened by the interplay of host and viral proteins.

In addition, as alluded to above, because we only examined one time point, the kinetics of the transcriptomic response provides another area for future analysis that could demonstrate interesting differences across species. These data will serve as an important baseline for such work, showing what transcriptional responses are possible and how they are impacted upon exposure to a pathogen. Furthermore, because the current study only examined poly-A–containing RNA transcripts, we cannot exclude the role of other RNA species, such as long non-coding RNAs that are becoming of increasing interest in NHP genomics ((45), reviewed in reference 46).

Finally, our parallel workflow allowed us to compare the benefits of mapping to the human versus the species-specific genome for these NHP species. The broader outcomes, such as the pathway enrichment analysis, were strongly similar between the two methods, and as noted earlier, the additional alignment gained by using species-specific genomes was often to regions with no features yet assigned. As a further comparison of the two mapping methods, we examined the top 500 most significantly differentially expressed genes ultimately identified after either the “species method” or the “human method.” Depending on the species, 88–93% of the top 500 most significant genes overlapped between the two mapping methods (Fig S11 and Supplemental Data 9). For each set of 500 genes from a given mapping method, we compared the log2(fold change) and significance for those same genes from the alternative mapping method. The log2(fold change) values were highly correlated and we observed for each species less than 15 genes that were (i) significantly differentially expressed by one approach (Padj ≥ 0.05) but not the other or (ii) had a log2(fold change) value that was 1.5 times higher or lower in one mapping method versus the other and also had a |log2(fold change)| ≥ 2 in at least one of the mapping methods. For genes that met these criteria, the differing expression between mapping methods was due to a variety of reasons. Reads which mapped to a region that had overlapping annotations in one of the reference genomes but not the other caused those reads to be excluded for not mapping to a unique feature. In other instances, multiple transcripts of the same gene were annotated covering a large region in one reference but not the other, so although alignment was observed by either mapping reference, an annotation only existed for one of the genomes.

Figure S11.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S11. Comparing the top 500 most significantly differentially expressed genes between the two mapping methods.

In the left column are the top 500 most significantly (based off Padj value) differentially expressed genes found after using the species-specific genome for mapping, with the genes arranged by increasing log2(foldchange) for the species mapping method (shown in purple). The corresponding log2(foldchange) for the same genes after using the human genome as the mapping reference are shown in green. The right column shows the inverse, with the coloring maintained, where the top 500 most significantly differentially expressed genes after mapping to the human genome are shown by increasing log2(foldchange). Each row is a different species and the Spearman coefficient of the log2(foldchange) values (after omitting any missing values of which there were no more than five) in each plot is shown. Data information: genes are labeled with their gene symbol if any of the following criteria were met: (i) the ratio of the log2(foldchange) from each method for a given gene was less than 0.67 or greater than 1.5 AND the absolute value of the log2(foldchange) by at least one of the mapping methods was ≥2, (ii) the log2(foldchange) was determined as “NA” by the mapping method not used to create the list of 500 most significant genes, and (iii) the Padj value was missing or nonsignificant. For genes in the final category, the point and label are shown enlarged. Where no gene symbol was listed, the human ENSEMBL ID was used.

Supplemental Data 9.

These are the log2(fold change) and Padj values for the genes in NHP species that were found to be in the top 500 most significantly differentially expressed genes after using the human genome mapping approach (file name includes “HumanTop500Sig”) versus those after using the species-specific genome mapping approach (file name includes “SpeciesTop500Sig”). In each file, the Padj and log2(fold change) values are shown for the same genes but by the other mapping approach. These are the data that were used to generate Fig S11.[LSA-2019-00495_Supplemental_Data_9.zip]

These findings are in line with the general observation that the human genome is more well annotated than those of NHPs and underscores the importance of improving the feature assignment of these sequences. RNA-Seq datasets such as ours can help identify novel exons or in some cases aid in rectifying feature assignments at regions where exons were incorrectly annotated. For example, as alluded to above, a given gene in the human genome may have more annotated transcript variants that can sometimes cover a larger region and lead to read assignment, whereas the orthologous NHP gene, with sometimes only one listed variant, has no reads assigned to it as the reads map but to an area outside of the limited region where the single variant was annotated. In making these data available, we anticipate that they will facilitate continued efforts to annotate the genomes of NHP species and to identify additional transcript variants of already annotated genes as well as completely novel transcripts unobserved in humans. As with the “1,000 Genomes Project” for human samples (47), having more sequencing data available from multiple individuals for a given species aids in identifying single-nucleotide polymorphisms and other sources of genomic variation to improve annotations. We explicitly strove to include at least one male and one female with varying ages across the three donors. This is far more reflective of the natural variance that exists in any outbred species population and makes our sequencing reads all the more useful, especially for genomes, such as that of orangutan, sourced from a single individual. We hope our present analysis and this data collection more broadly will serve as a springboard for additional evolutionary analysis and comparisons.

Materials and Methods

DFs

Mouse DFs from C57BL/6 mice were purchased from Cell Biologics. Primary human DFs were obtained from American Type Culture Collection (PCS-201-012). All the NHP DFs were purchased from the Coriell Institute for Medical Research with the exception of two squirrel monkey donors (SQMA and SQMB), which were derived from skin biopsies generously provided by Robert Lanford (Southwest Biomedical Research Center). For any of the DF donor species that are on the United States Fish and Wildlife Services endangered or threatened species list, the Coriell Institute has the appropriate and required documentation of breeding records (Supplemental File 1) indicating captivity before or birth after November 18, 1976. The squirrel monkey skin biopsies were obtained after terminal necropsy using Princeton University Institutional Animal Care and Use Committee–approved protocols (#1930) and shipped overnight on wet ice to Princeton University. Upon arrival, the skin was prepped and DFs isolated according to previously published protocols (48). In brief, the skin biopsies were scraped to remove connective tissue, cut into smaller pieces, and digested overnight at 4°C in HBSS without Ca2+ and Mg2+ (Thermo Fisher Scientific), containing 1 ml dispase (5,000 caseinolytic units/ml; Corning) for every 9 ml of HBSS containing final concentrations of 100 mg/ml streptomycin, 100 U/ml penicillin, and 250 ng/ml amphotericin B (HyClone). After digestion, the epidermis was removed and discarded, whereas the remaining dermis was cut into smaller pieces less than a few square millimeter in area. These pieces were moistened with DMEM and pressed into a six-well plate scored with a razor blade. The dermis was maintained in DMEM containing 10% FBS and 1% vol/vol penicillin/streptomycin solution at 37°C, 5% CO2. Media was changed every 4–5 d and fibroblast growth was typically observed within 1 wk of culture. Once sufficient outgrowth had occurred, the dermis was removed from the plate and the fibroblasts expanded into larger cultures. Complete donor information as provided by these sources can be found in Table 4.

Supplemental File 1.

The breeding records and documentation for DFs acquired from species on the United States Fish and Wildlife Services endangered or threatened species list.[LSA-2019-00495_Supplemental_file_1.pdf]

Cell culture

All cells, unless otherwise stated, were grown under standard conditions in DMEM (Thermo Fisher Scientific) containing 10% FBS and 1% vol/vol penicillin/streptomycin. Upon reaching confluency, the cells were trypsinized with 0.05% trypsin/EDTA and split 1:3.

polyI:C transfections

DFs were transfected with ∼0.05 μg of high molecular weight poly(I:C) (Invivogen) per square centimeter (0.5 μg/well in a six-well format in triplicate for samples undergoing RNA-Seq and 0.1 μg/well in duplicate in a 24-well format for relative RT-qPCR validation experiments) using X-tremeGENE HP DNA Transfection Reagent (Roche) 1 μl per μg of poly(I:C) in Opti-MEM (Thermo Fisher Scientific). Mock transfections were performed in parallel under the same conditions minus poly(I:C). Collected cell lysates (350 μl volume) were immediately frozen at −80°C until RNA extraction performed.

RNA extraction, cDNA library preparation, and RNA sequencing

Total RNA was isolated from DFs using the RNeasy Mini Kit (QIAGEN) according to the manufacturer’s instructions. For all samples undergoing RNA-Seq, the quality and concentration of RNA was assessed on an Agilent 2100 Bioanalyzer (Agilent Technologies). All samples subsequently sequenced had an RNA Integrity Number of ≥8. The poly-A–containing RNA transcripts in the total RNA samples were converted to cDNA and amplified after the Smart-seq2 method (49). Sequencing libraries were made from the amplified cDNA samples using the Nextera kit (Illumina), assigning a unique barcode to each of the libraries to be sequenced together. The cDNA samples and libraries were examined on the Bioanalyzer (Agilent) DNA HS chips for size distribution and quantified by Qubit fluorometer (Invitrogen). The RNA-Seq libraries were pooled together in equal amounts and sequenced on the Illumina HiSeq 2500 Rapid Flowcell as single-end 75-nt reads following the standard protocol, giving a range of 15–20 million reads per sample. Raw sequencing reads were filtered by Illumina HiSeq Control Software and only pass-filter reads were used for further analysis.

RT-qPCR of select ISGs

RT-qPCR of total RNA isolated from the 24-well format poly(I:C) transfections was performed using the Luna Universal One-Step RT-qPCR kit (New England BioLabs, Inc.) according to the manufacturer’s directions. Primer sequences can be found in Table 6. In brief, a master mix was prepared that comprised 10 μl of 2× Luna Universal One-Step Reaction Mix (2×), 1 μl of 20× Luna WarmStart RT Enzyme Mix, 0.8 μl of a 10 μM stock of each primer, and 5.4 μl of nuclease-free water per reaction. Each well received 18 μl of the appropriate master mix and 2 μl of the RNA being assayed. The following PCR program was then run on an Applied Biosystems Step One Plus qPCR machine (Life Technologies): denatured at 55°C for 10 min, 95°C for 1 min, followed by 40 cycles of 95°C for 10 s and 60°C for 60 s. Last, a melt curve was performed at 95°C for 10 s, 65°C for 10 s, 95°C for 10 s, and 50°C for 5 s.

View this table:
  • View inline
  • View popup
Table 6.

Primers for RT-qPCR of IFNβ, MX1, OASL, and HPRT1 for the species used in this study.

RNA-Seq analysis

Using the Galaxy system (50) provided by Princeton University, short reads were aligned to the human or the originating species’ genomes (see Table 2) using RNA STAR (51, 52) (Galaxy version 2.6.0b-1) with default parameters. Counts were generated using featureCounts with default settings (53) (Galaxy version 1.6.3+galaxy2), downloaded, and read into R (54) version 3.5.2 (December 20, 2018) using scripts run in RStudio (55) version 1.1.463. DGE was determined using DESeq2 (56) (version 1.22.2) using the standard filters and nbinomWaldTest (57). The design used to model the samples was ∼species + species:donor.n + treatment + species:treatment with various contrasts set as shown in the R code and as described in the README files for the Datasets EV3–5. Transcript counts were normalized using DESeq2 default options and transformed using the regularized log2 function in DESeq2 before PCA plotting. Results were extracted from the DESeq2 analysis and annotated using Bioconductor’s AnnotationDbi (version 1.44.0) org.Hs.eg.db and org.Mm.eg.db. For further details concerning the R packages used and the specific conditions used for analysis, R code can be accessed at https://github.com/aploss/polyIC-dermal-fibroblasts-RNA-Seq.

Statistical analysis

Statistical analysis for RT-qPCR data was performed with GraphPad Prism software as indicated in the figure legends.

Data access

The RNA-Seq data generated in this study are deposited in the National Center for Biotechnology Information Gene Expression Omnibus database (accession number GSE105160).

Acknowledgements

This work was supported by grants from the National Institutes of Health (R01 AI107301 to A Ploss, R21AI117213 to A Ploss and RE Schwartz), a Research Scholar Award from the American Cancer Society (RSG-15-048-01-MPC to A Ploss), Princeton University, a research grant through the Princeton University Center for Health and Wellbeing Health Grand Challenge Program to A Ploss and JM Gaska, and an Investigator in Pathogenesis Award by the Burroughs Wellcome Fund (101539) to A Ploss. JM Gaska was supported by the National Institute of General Medical Sciences of the National Institutes of Health under grant number T32GM007388. We also thank Jennifer Miller and Jessica Wiggins of the Lewis-Sigler Institute for Integrative Genomics for their help in assessing the RNA quality, preparing the cDNA libraries, and performing RNA-Seq. We thank members of the Ploss lab for critical discussions and comments on this project. We thank Matthew Cahn (Princeton University) for his assistance in facilitating the submission of our raw and processed data to National Center for Biotechnology Information Gene Expression Omnibus (accession number GSE105160).

Author Contributions

  • JM Gaska: conceptualization, formal analysis, funding acquisition, investigation, visualization, methodology, and writing—original draft and editing.

  • L Parsons: formal analysis, methodology, and writing—review and editing.

  • M Balev: formal analysis.

  • A Cirincione: formal analysis.

  • W Wang: formal analysis, investigation, and methodology.

  • RE Schwartz: conceptualization.

  • A Ploss: conceptualization, supervision, funding acquisition, project administration, and writing—original draft, review, and editing.

Conflict of Interest Statement

The authors declare that they have no conflict of interest.

  • Received July 20, 2019.
  • Revision received October 14, 2019.
  • Accepted October 15, 2019.
  • © 2019 Gaska et al.
Creative Commons logoCreative Commons logohttps://creativecommons.org/licenses/by/4.0/

This article is available under a Creative Commons License (Attribution 4.0 International, as described at https://creativecommons.org/licenses/by/4.0/).

References

  1. 1.↵
    1. Douam F,
    2. Gaska JM,
    3. Winer BY,
    4. Ding Q,
    5. von Schaewen M,
    6. Ploss A
    (2015) Genetic dissection of the host tropism of human-tropic pathogens. Annu Rev Genet 49: 21–45. doi:10.1146/annurev-genet-112414-054823
    OpenUrlCrossRef
  2. 2.↵
    1. Deschamps M,
    2. Laval G,
    3. Fagny M,
    4. Itan Y,
    5. Abel L,
    6. Casanova JL,
    7. Patin E,
    8. Quintana-Murci L
    (2016) Genomic signatures of selective pressures and introgression from archaic hominins at human innate immunity genes. Am J Hum Genet 98: 5–21. doi:10.1016/j.ajhg.2015.11.014
    OpenUrlCrossRefPubMed
  3. 3.↵
    1. Nielsen R,
    2. Bustamante C,
    3. Clark AG,
    4. Glanowski S,
    5. Sackton TB,
    6. Hubisz MJ,
    7. Fledel-Alon A,
    8. Tanenbaum DM,
    9. Civello D,
    10. White TJ, et al.
    (2005) A scan for positively selected genes in the genomes of humans and chimpanzees. PLoS Biol 3: e170. doi:10.1371/journal.pbio.0030170
    OpenUrlCrossRefPubMed
  4. 4.↵
    1. Voight BF,
    2. Kudaravalli S,
    3. Wen X,
    4. Pritchard JK
    (2006) A map of recent positive selection in the human genome. PLoS Biol 4: e72. doi:10.1371/journal.pbio.0040154
    OpenUrlCrossRefPubMed
  5. 5.↵
    1. Enard D,
    2. Depaulis F,
    3. Roest Crollius H
    (2010) Human and non-human primate genomes share hotspots of positive selection. PLoS Genet 6: e1000840. doi:10.1371/journal.pgen.1000840
    OpenUrlCrossRefPubMed
  6. 6.↵
    1. Sironi M,
    2. Cagliani R,
    3. Forni D,
    4. Clerici M
    (2015) Evolutionary insights into host-pathogen interactions from mammalian sequence data. Nat Rev Genet 16: 224–236. doi:10.1038/nrg3905
    OpenUrlCrossRefPubMed
  7. 7.↵
    1. Sheahan T,
    2. Jones CT,
    3. Ploss A
    (2010) Advances and challenges in studying hepatitis C virus in its native environment. Expert Rev Gastroenterol Hepatol 4: 541–550. doi:10.1586/egh.10.53
    OpenUrlCrossRefPubMed
  8. 8.↵
    1. Wilkening S,
    2. Stahl F,
    3. Bader A
    (2003) Comparison of primary human hepatocytes and hepatoma cell line Hepg2 with regard to their biotransformation properties. Drug Metab Dispos 31: 1035–1042. doi:10.1124/dmd.31.8.1035
    OpenUrlAbstract/FREE Full Text
  9. 9.↵
    1. Sumpter R, Jr.,
    2. Loo YM,
    3. Foy E,
    4. Li K,
    5. Yoneyama M,
    6. Fujita T,
    7. Lemon SM,
    8. Gale M, Jr.
    (2005) Regulating intracellular antiviral defense and permissiveness to hepatitis C virus RNA replication through a cellular RNA helicase, RIG-I. J Virol 79: 2689–2699. doi:10.1128/jvi.79.5.2689-2699.2005
    OpenUrlAbstract/FREE Full Text
  10. 10.↵
    1. De Clercq E
    (2006) Interferon and its inducers: A never-ending story: “Old” and “new” data in a new perspective. J Infect Dis 194: S19–S26. doi:10.1086/505351
    OpenUrlCrossRefPubMed
  11. 11.↵
    1. Field AK,
    2. Tytell AA,
    3. Lampson GP,
    4. Hilleman MR
    (1967) Inducers of interferon and host resistance. II. Multistranded synthetic polynucleotide complexes. Proc Natl Acad Sci U S A 58: 1004–1010. doi:10.1073/pnas.58.3.1004
    OpenUrlFREE Full Text
  12. 12.↵
    1. Areal H,
    2. Abrantes J,
    3. Esteves PJ
    (2011) Signatures of positive selection in toll-like receptor (TLR) genes in mammals. BMC Evol Biol 11: 368. doi:10.1186/1471-2148-11-368
    OpenUrlCrossRefPubMed
  13. 13.↵
    1. Wlasiuk G,
    2. Nachman MW
    (2010) Adaptation and constraint at toll-like receptors in primates. Mol Biol Evol 27: 2172–2186. doi:10.1093/molbev/msq104
    OpenUrlCrossRefPubMed
  14. 14.↵
    1. Carter DL,
    2. Shieh TM,
    3. Blosser RL,
    4. Chadwick KR,
    5. Margolick JB,
    6. Hildreth JE,
    7. Clements JE,
    8. Zink MC
    (1999) CD56 identifies monocytes and not natural killer cells in rhesus macaques. Cytometry 37: 41–50. doi:10.1002/(sici)1097-0320(19990901)37:1<41::aid-cyto5>3.0.co;2-4
    OpenUrlCrossRefPubMed
  15. 15.↵
    1. Brown KN,
    2. Barratt-Boyes SM
    (2009) Surface phenotype and rapid quantification of blood dendritic cell subsets in the rhesus macaque. J Med Primatol 38: 272–278. doi:10.1111/j.1600-0684.2009.00353.x
    OpenUrlCrossRefPubMed
  16. 16.↵
    1. Steiper ME,
    2. Young NM
    (2006) Primate molecular divergence dates. Mol Phylogenet Evol 41: 384–394. doi:10.1016/j.ympev.2006.05.021
    OpenUrlCrossRefPubMed
  17. 17.↵
    1. Fabre PH,
    2. Rodrigues A,
    3. Douzery EJ
    (2009) Patterns of macroevolution among primates inferred from a supermatrix of mitochondrial and nuclear DNA. Mol Phylogenet Evol 53: 808–825. doi:10.1016/j.ympev.2009.08.004
    OpenUrlCrossRefPubMed
  18. 18.↵
    1. Webster RL,
    2. Johnson RP
    (2005) Delineation of multiple subpopulations of natural killer cells in rhesus macaques. Immunology 115: 206–214. doi:10.1111/j.1365-2567.2005.02147.x
    OpenUrlCrossRefPubMed
  19. 19.↵
    1. Rogers KA,
    2. Jayashankar L,
    3. Scinicariello F,
    4. Attanasio R
    (2008) Nonhuman primate IgA: Genetic heterogeneity and interactions with CD89. J Immunol 180: 4816–4824. doi:10.4049/jimmunol.180.7.4816
    OpenUrlAbstract/FREE Full Text
  20. 20.↵
    1. Kita YF,
    2. Hosomichi K,
    3. Kohara S,
    4. Itoh Y,
    5. Ogasawara K,
    6. Tsuchiya H,
    7. Torii R,
    8. Inoko H,
    9. Blancher A,
    10. Kulski JK, et al.
    (2009) MHC class I A loci polymorphism and diversity in three Southeast Asian populations of cynomolgus macaque. Immunogenetics 61: 635–648. doi:10.1007/s00251-009-0390-y
    OpenUrlCrossRef
  21. 21.↵
    1. Otting N,
    2. de Vos-Rouweler AJ,
    3. Heijmans CM,
    4. de Groot NG,
    5. Doxiadis GG,
    6. Bontrop RE
    (2007) MHC class I A region diversity and polymorphism in macaque species. Immunogenetics 59: 367–375. doi:10.1007/s00251-007-0201-2
    OpenUrlCrossRefPubMed
  22. 22.↵
    1. Zhang GQ,
    2. Ni C,
    3. Ling F,
    4. Qiu W,
    5. Wang HB,
    6. Xiao Y,
    7. Guo XJ,
    8. Huang JY,
    9. Du HL,
    10. Wang JF, et al.
    (2012) Characterization of the major histocompatibility complex class I A alleles in cynomolgus macaques of Vietnamese origin. Tissue Antigens 80: 494–501. doi:10.1111/tan.12024
    OpenUrlCrossRefPubMed
  23. 23.↵
    1. Ramesh A,
    2. Darko S,
    3. Hua A,
    4. Overman G,
    5. Ransier A,
    6. Francica JR,
    7. Trama A,
    8. Tomaras GD,
    9. Haynes BF,
    10. Douek DC, et al.
    (2017) Structure and diversity of the rhesus macaque immunoglobulin loci through multiple de novo genome assemblies. Front Immunol 8: 1407. doi:10.3389/fimmu.2017.01407
    OpenUrlCrossRef
  24. 24.↵
    1. Ellis-Connell AL,
    2. Kannal NM,
    3. Balgeman AJ,
    4. O’Connor SL
    (2019) Characterization of major histocompatibility complex-related molecule 1 sequence variants in non-human primates. Immunogenetics 71: 109–121. doi:10.1007/s00251-018-1091-1
    OpenUrlCrossRef
  25. 25.↵
    1. Barreiro LB,
    2. Marioni JC,
    3. Blekhman R,
    4. Stephens M,
    5. Gilad Y
    (2010) Functional comparison of innate immune signaling pathways in primates. PLoS Genet 6: e1001249. doi:10.1371/journal.pgen.1001249
    OpenUrlCrossRefPubMed
  26. 26.↵
    1. Hagai T,
    2. Chen X,
    3. Miragaia RJ,
    4. Rostom R,
    5. Gomes T,
    6. Kunowska N,
    7. Henriksson J,
    8. Park JE,
    9. Proserpio V,
    10. Donati G, et al.
    (2018) Gene expression variability across cells and species shapes innate immunity. Nature 563: 197–202. doi:10.1038/s41586-018-0657-2
    OpenUrlCrossRef
  27. 27.↵
    1. Shaw AE,
    2. Hughes J,
    3. Gu Q,
    4. Behdenna A,
    5. Singer JB,
    6. Dennis T,
    7. Orton RJ,
    8. Varela M,
    9. Gifford RJ,
    10. Wilson SJ, et al.
    (2017) Fundamental properties of the mammalian innate immune system revealed by multispecies comparison of type I interferon responses. PLoS Biol 15: e2004086. doi:10.1371/journal.pbio.2004086
    OpenUrlCrossRefPubMed
  28. 28.↵
    1. de Veer MJ,
    2. Holko M,
    3. Frevel M,
    4. Walker E,
    5. Der S,
    6. Paranjape JM,
    7. Silverman RH,
    8. Williams BR
    (2001) Functional classification of interferon-stimulated genes identified using microarrays. J Leukoc Biol 69: 912–920. doi:10.1189/jlb.69.6.912
    OpenUrlCrossRefPubMed
  29. 29.↵
    1. Schoggins JW,
    2. Wilson SJ,
    3. Panis M,
    4. Murphy MY,
    5. Jones CT,
    6. Bieniasz P,
    7. Rice CM
    (2011) A diverse range of gene products are effectors of the type I interferon antiviral response. Nature 472: 481–485. doi:10.1038/nature09907
    OpenUrlCrossRefPubMed
  30. 30.↵
    1. Indraccolo S,
    2. Pfeffer U,
    3. Minuzzo S,
    4. Esposito G,
    5. Roni V,
    6. Mandruzzato S,
    7. Ferrari N,
    8. Anfosso L,
    9. Dell’Eva R,
    10. Noonan DM, et al.
    (2007) Identification of genes selectively regulated by IFNs in endothelial cells. J Immunol 178: 1122–1135. doi:10.4049/jimmunol.178.2.1122
    OpenUrlAbstract/FREE Full Text
  31. 31.↵
    1. Benjamin AM,
    2. Nichols M,
    3. Burke TW,
    4. Ginsburg GS,
    5. Lucas JE
    (2014) Comparing reference-based RNA-Seq mapping methods for non-human primate data. BMC Genomics 15: 570. doi:10.1186/1471-2164-15-570
    OpenUrlCrossRef
  32. 32.↵
    1. Peng X,
    2. Thierry-Mieg J,
    3. Thierry-Mieg D,
    4. Nishida A,
    5. Pipes L,
    6. Bozinoski M,
    7. Thomas MJ,
    8. Kelly S,
    9. Weiss JM,
    10. Raveendran M, et al.
    (2015) Tissue-specific transcriptome sequencing analysis expands the non-human primate reference transcriptome resource (NHPRTR). Nucleic Acids Res 43: D737–D742. doi:10.1093/nar/gku1110
    OpenUrlCrossRefPubMed
  33. 33.↵
    1. Pipes L,
    2. Li S,
    3. Bozinoski M,
    4. Palermo R,
    5. Peng X,
    6. Blood P,
    7. Kelly S,
    8. Weiss JM,
    9. Thierry-Mieg J,
    10. Thierry-Mieg D, et al.
    (2013) The non-human primate reference transcriptome resource (NHPRTR) for comparative functional genomics. Nucleic Acids Res 41: D906–D914. doi:10.1093/nar/gks1268
    OpenUrlCrossRefPubMed
  34. 34.↵
    1. Hornett EA,
    2. Wheat CW
    (2012) Quantitative RNA-seq analysis in non-model species: Assessing transcriptome assemblies as a scaffold and the utility of evolutionary divergent genomic reference species. BMC Genomics 13: 361. doi:10.1186/1471-2164-13-361
    OpenUrlCrossRefPubMed
  35. 35.↵
    1. Breuer K,
    2. Foroushani AK,
    3. Laird MR,
    4. Chen C,
    5. Sribnaia A,
    6. Lo R,
    7. Winsor GL,
    8. Hancock RE,
    9. Brinkman FS,
    10. Lynn DJ
    (2013) InnateDB: Systems biology of innate immunity and beyond--recent updates and continuing curation. Nucleic Acids Res 41: D1228–D1233. doi:10.1093/nar/gks1147
    OpenUrlCrossRefPubMed
  36. 36.↵
    1. Hung JH,
    2. Yang TH,
    3. Hu Z,
    4. Weng Z,
    5. DeLisi C
    (2012) Gene set enrichment analysis: Performance evaluation and usage guidelines. Brief Bioinform 13: 281–291. doi:10.1093/bib/bbr049
    OpenUrlCrossRefPubMed
  37. 37.↵
    1. Subramanian A,
    2. Tamayo P,
    3. Mootha VK,
    4. Mukherjee S,
    5. Ebert BL,
    6. Gillette MA,
    7. Paulovich A,
    8. Pomeroy SL,
    9. Golub TR,
    10. Lander ES, et al.
    (2005) Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A 102: 15545–15550. doi:10.1073/pnas.0506580102
    OpenUrlAbstract/FREE Full Text
  38. 38.↵
    1. Khatri P,
    2. Sirota M,
    3. Butte AJ
    (2012) Ten years of pathway analysis: Current approaches and outstanding challenges. PLoS Comput Biol 8: e1002375. doi:10.1371/journal.pcbi.1002375
    OpenUrlCrossRefPubMed
  39. 39.↵
    1. de Leeuw CA,
    2. Neale BM,
    3. Heskes T,
    4. Posthuma D
    (2016) The statistical properties of gene-set analysis. Nat Rev Genet 17: 353–364. doi:10.1038/nrg.2016.29
    OpenUrlCrossRefPubMed
  40. 40.↵
    1. Luo W,
    2. Friedman MS,
    3. Shedden K,
    4. Hankenson KD,
    5. Woolf PJ
    (2009) GAGE: Generally applicable gene set enrichment for pathway analysis. BMC Bioinformatics 10: 161. doi:10.1186/1471-2105-10-161
    OpenUrlCrossRefPubMed
  41. 41.↵
    1. Kanehisa M,
    2. Furumichi M,
    3. Tanabe M,
    4. Sato Y,
    5. Morishima K
    (2017) KEGG: New perspectives on genomes, pathways, diseases and drugs. Nucleic Acids Res 45: D353–D361. doi:10.1093/nar/gkw1092
    OpenUrlCrossRefPubMed
  42. 42.↵
    1. Kanehisa M,
    2. Goto S
    (2000) KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res 28: 27–30. doi:10.1093/nar/28.1.27
    OpenUrlCrossRefPubMed
  43. 43.↵
    1. Kanehisa M,
    2. Sato Y,
    3. Kawashima M,
    4. Furumichi M,
    5. Tanabe M
    (2016) KEGG as a reference resource for gene and protein annotation. Nucleic Acids Res 44: D457–D462. doi:10.1093/nar/gkv1070
    OpenUrlCrossRefPubMed
  44. 44.↵
    1. Takami M,
    2. Kim N,
    3. Rho J,
    4. Choi Y
    (2002) Stimulation by toll-like receptors inhibits osteoclast differentiation. J Immunol 169: 1516–1523. doi:10.4049/jimmunol.169.3.1516
    OpenUrlAbstract/FREE Full Text
  45. 45.↵
    1. Field AR,
    2. Jacobs FMJ,
    3. Fiddes IT,
    4. Phillips APR,
    5. Reyes-Ortiz AM,
    6. LaMontagne E,
    7. Whitehead L,
    8. Meng V,
    9. Rosenkrantz JL,
    10. Olsen M, et al.
    (2019) Structurally conserved primate LncRNAs are transiently expressed during human cortical differentiation and influence cell-type-specific genes. Stem Cell Rep 12: 245–257. doi:10.1016/j.stemcr.2018.12.006
    OpenUrlCrossRef
  46. 46.↵
    1. Awan HM,
    2. Shah A,
    3. Rashid F,
    4. Shan G
    (2017) Primate-specific long non-coding RNAs and microRNAs. Genomics Proteomics Bioinformatics 15: 187–195. doi:10.1016/j.gpb.2017.04.002
    OpenUrlCrossRef
  47. 47.↵
    1. 1000 Genomes Project Consortium,
    2. Auton A,
    3. Brooks LD,
    4. Durbin RM,
    5. Garrison EP,
    6. Kang HM,
    7. Korbel JO,
    8. Marchini JL,
    9. McCarthy S,
    10. McVean GA, et al.
    (2015) A global reference for human genetic variation. Nature 526: 68–74. doi:10.1038/nature15393
    OpenUrlCrossRefPubMed
  48. 48.↵
    1. Aasen T,
    2. Izpisua Belmonte JC
    (2010) Isolation and cultivation of human keratinocytes from skin or plucked hair for the generation of induced pluripotent stem cells. Nat Protoc 5: 371–382. doi:10.1038/nprot.2009.241
    OpenUrlCrossRefPubMed
  49. 49.↵
    1. Picelli S,
    2. Faridani OR,
    3. Bjorklund AK,
    4. Winberg G,
    5. Sagasser S,
    6. Sandberg R
    (2014) Full-length RNA-seq from single cells using smart-seq2. Nat Protoc 9: 171–181. doi:10.1038/nprot.2014.006
    OpenUrlCrossRefPubMed
  50. 50.↵
    1. Afgan E,
    2. Baker D,
    3. van den Beek M,
    4. Blankenberg D,
    5. Bouvier D,
    6. Cech M,
    7. Chilton J,
    8. Clements D,
    9. Coraor N,
    10. Eberhard C, et al.
    (2016) The galaxy platform for accessible, reproducible and collaborative biomedical analyses: 2016 update. Nucleic Acids Res 44: W3–W10. doi:10.1093/nar/gkw343
    OpenUrlCrossRefPubMed
  51. 51.↵
    1. Dobin A,
    2. Davis CA,
    3. Schlesinger F,
    4. Drenkow J,
    5. Zaleski C,
    6. Jha S,
    7. Batut P,
    8. Chaisson M,
    9. Gingeras TR
    (2013) STAR: Ultrafast universal RNA-seq aligner. Bioinformatics 29: 15–21. doi:10.1093/bioinformatics/bts635
    OpenUrlCrossRefPubMed
  52. 52.↵
    1. Dobin A,
    2. Gingeras TR
    (2015) Mapping RNA-seq reads with STAR. Curr Protoc Bioinformatics 51: 11.14.1–19. doi:10.1002/0471250953.bi1114s51
    OpenUrlCrossRefPubMed
  53. 53.↵
    1. Liao J,
    2. Li X,
    3. Wong TY,
    4. Wang JJ,
    5. Khor CC,
    6. Tai ES,
    7. Aung T,
    8. Teo YY,
    9. Cheng CY
    (2014) Impact of measurement error on testing genetic association with quantitative traits. PLoS One 9: e87044. doi:10.1371/journal.pone.0087044
    OpenUrlCrossRefPubMed
  54. 54.↵
    1. R Core Team
    (2018) R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing. https://www.R-project.org/
  55. 55.↵
    1. R Studio Team
    (2016) RStudio: Integrated Development for R. Boston, MA: RStudio, Inc. http://www.rstudio.com/
  56. 56.↵
    1. Love MI,
    2. Huber W,
    3. Anders S
    (2014) Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol 15: 550. doi:10.1186/s13059-014-0550-8
    OpenUrlCrossRefPubMed
  57. 57.↵
    1. Anders S,
    2. Pyl PT,
    3. Huber W
    (2015) HTSeq: A python framework to work with high-throughput sequencing data. Bioinformatics 31: 166–169. doi:10.1093/bioinformatics/btu638
    OpenUrlCrossRefPubMed
PreviousNext
Back to top
Download PDF
Article Alerts
Sign In to Email Alerts with your Email Address
Email Article

Thank you for your interest in spreading the word on Life Science Alliance.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Conservation of cell-intrinsic immune responses in diverse nonhuman primate species
(Your Name) has sent you a message from Life Science Alliance
(Your Name) thought you would like to see the Life Science Alliance web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Primate cell-intrinsic immunity
Jenna M Gaska, Lance Parsons, Metodi Balev, Ann Cirincione, Wei Wang, Robert E Schwartz, Alexander Ploss
Life Science Alliance Oct 2019, 2 (5) e201900495; DOI: 10.26508/lsa.201900495

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Share
Primate cell-intrinsic immunity
Jenna M Gaska, Lance Parsons, Metodi Balev, Ann Cirincione, Wei Wang, Robert E Schwartz, Alexander Ploss
Life Science Alliance Oct 2019, 2 (5) e201900495; DOI: 10.26508/lsa.201900495
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
Issue Cover

In this Issue

Volume 2, No. 5
October 2019
  • Table of Contents
  • Cover (PDF)
  • About the Cover
  • Masthead (PDF)
Advertisement

Jump to section

  • Article
    • Abstract
    • Introduction
    • Results and Discussion
    • Materials and Methods
    • Acknowledgements
    • References
  • Figures & Data
  • Info
  • Metrics
  • Reviewer Comments
  • PDF

Subjects

  • Immunology
  • Microbiology, Virology & Host Pathogen Interaction
  • Molecular Biology

Related Articles

  • No related articles found.

Cited By...

  • No citing articles found.
  • Google Scholar

More in this TOC Section

  • Proteome profiling in heart valve disease
  • IARA: a curated genetic atlas of human spliceosome
  • Concordance of MERFISH with single-cell RNA-Seq
Show more Resource

Similar Articles

EMBO Press LogoRockefeller University Press LogoCold Spring Harbor Logo

Content

  • Home
  • Newest Articles
  • Current Issue
  • Archive
  • Subject Collections

For Authors

  • Submit a Manuscript
  • Author Guidelines
  • License, copyright, Fee

Other Services

  • Alerts
  • Twitter
  • RSS Feeds

More Information

  • Editors & Staff
  • Reviewer Guidelines
  • Feedback
  • Licensing and Reuse
  • Privacy Policy

ISSN: 2575-1077
© 2023 Life Science Alliance LLC

Life Science Alliance is registered as a trademark in the U.S. Patent and Trade Mark Office and in the European Union Intellectual Property Office.