Characteristics (including hematology values) and treatments (with or without LPS) of Mus musculus CD-1 and Peromyscus leucopus LL stock

Total white blood cells, neutrophils, and lymphocytes of Mus musculus (M) and Peromyscus leucopus (P) with or without (control; C) treatment with 10 µg lipopolysaccharide (LPS; L) per g body mass 4 h previous. The box plots of left and center panels show values of individual animals and compiled median, quartiles, and extreme values. The linear regressions of the right panel are color-coded according to the species and treatment designations. The outlier value for a M. musculus control (MM17) was excluded from the linear regression for that group.

Principle component analysis of genome-wide RNA-seq data of P. leucopus or M. musculus with or without (blue dot) treatment with LPS 4 h previous. The individual animals listed in Table 1 are indicated on the graphs. The insets indicate the size and color of the symbol for the experimental condition (LPS-treated or control).

Gene Ontology (GO) term clusters associated with up-regulated genes (upper panels) and down-regulated genes (lower panels) of P. leucopus (left panels) and M. musculus (right panels) treated with LPS in comparison with untreated controls of each species. The scale for the x-axes for the panels was determined by the highest -log10 p values in each of the 4 sets. The horizontal bar color, which ranges from white to dark brown through shades of yellow through orange in between, is a schematic representation of the -log10 p values.

Scatter plot with linear regression of pairs of log2-transformed mean fold-changes between LPS-treated and control P. leucopus by male and female sex. The 5012 reference transcripts from the genome reference set are defined in the text and listed in Table S2 and Dryad Table D5. The coefficient of determination (R2), the 95% upper and lower prediction limits for the regression line, and distributions of the values on the x-and y-axes are shown. Selected genes for which their x-y coordinates fall outside the limits of prediction are labeled. Cxcl2, Ibsp, Saa3, Saa5, Sbno2, Serpine1, Slpi, and Steap1 were noted as up-regulated DEGs for the groups with both sexes (Table S1).

Scatter plots with linear regression of pairs of log-transformed (ln) normalized RNA-seq reads for selected coding sequences for control P. leucopus and M. musculus (left panel) and LPS-treated P. leucopus and M. musculus (right panel). The R2 values and selected genes (each with a different symbol) are indicated in each graph.

Targeted RNA-seq of blood of Peromyscus leucopus and Mus musculus 4 hours after intraperitoneal injection of lipopolysaccharide (LPS) or saline control

RT-qPCR of blood of LPS-treated P.leucopus and M.musculus

Taegeted RNA-seq of P. leucopus blood in 12 h and 1 ug/g LPS experiment

Box plots of log-transformed ratios of four pairs of gene transcripts from targeted RNA-seq analysis of blood of P. leucopus (P) or M. musculus (M) with (L) or without (C) treatment with LPS. The values are from Table S4. Upper left, Nos2/Arg1. Upper right, IL12/IL10. Lower left, Akt1/Akt2. Lower right, Cd14/Cd16.

Transcripts of interferon-gamma and interleukin-1 beta by targeted RNA-seq of the blood of P. leucopus (P) or M. musculus (M) with (L) or without (C) treatment with LPS. The top panels are box plots of the individual values. The lower left panel is a scatter plot of interleukin-1 ý on interferon-ψ values. The lower right panel is a Discriminant Analysis of these pairs of values where Factor 1 corresponds to interferon-gamma, and Factor 2 corresponds to interleukin-1 beta. Values for analysis are in Table S4.

Hematology, cytokines, and targeted RNA-seq of LPS-treated and control rats

Normalized transcripts of Nos2, Ifng, and Cd69 in targeted RNA-seq analysis of blood of P. leucopus (P) or M. musculus (M) with (L) or without (C) treatment with LPS. Upper left: scatter plot of individual values for Nos2 on Ifng with linear regression curve and coefficient of determination (R2). Upper right: natural logarithm (ln) of ratios of Nos2/Argi1 on Ifng/IL1b with regression curve and R2. Lower left: Box plots of individual values of normalized transcripts of Cd69. Lower right: Scatter plot of Ifng on Cd69 with separate regression curves and R2 values for M. musculus and P. leucopus. Values for analysis are in Table S4. Box plots for Nos2 and Arg1 are provided in Figure S5, and box plots for Ifng and Il1b are provided in Figure 7.

Co-variation between transcripts for selected PRRs and ISGs in the blood of P. leucopus (P) or M. musculus (M) with (L) or without (C) LPS treatment. Top panel: matrix of coefficients of determination (R2) for combined P. leucopus and M. musculus data. PRRs are indicated by yellow fill and ISGs by blue fill on horizontal and vertical axes. Shades of green of the matrix cells correspond to R2 values, where cells with values less than 0.30 have white fill and those of 0.90-1.00 have deepest green fill. Bottom panels: scatter plots of log-transformed normalized Mx2 transcripts on Rigi (left), Ifih1 (center), and Gbp4 (right). The linear regression curves are for each species. For the right-lower graph the result from the General Linear Model (GLM) estimate is also given. Values for analysis are in Table S4; box plots for Gbp4, Irf7, Isg15, Mx2, and Oas1 are provided in Figure S6.

Scatter plots of endogenous retrovirus (ERV) Env and Gag-pol protein gene transcription (left) and association of ERV Env with Rigi transcription (right) in the blood of P. leucopus (Pero; P), M. musculus (Mus; M), or R. norvegicus (Rattus) with (L) or without (control; C) treatment with LPS. In right panel the linear regression curve and coefficients of determination (R2) for P. leucopus and M. musculus are shown. Values for analysis are in Table S4; box plots for ERV Env and ERV Gag-pol are provided in Figure S7.

Targeted RNA-seq of P. leucopus with and without Borrelia hermsii infection

Summary of distinguishing features of transcriptional responses in the blood between P. leucopus and M. musculus 4 h after treatment with LPS. There is semi-quantitative representation of relative transcription of selected coding sequences or ratios of transcription for selected pairs of genes in the blood.