Skip to main content

ORIGINAL RESEARCH article

Front. Immunol., 12 December 2022
Sec. Viral Immunology
This article is part of the Research Topic Host-Virus Interaction at the Omics and Ecology Levels View all 7 articles

Chikungunya patient transcriptional signatures faithfully recapitulated in a C57BL/6J mouse model

  • 1Department of Infection and Inflammation, Queensland Institute of Medical Research, Berghofer Medical Research Institute, Brisbane, QLD, Australia
  • 2Pathology Advanced Translational Research Unit, Department of Pathology and Laboratory Medicine, Emory University School of Medicine, Atlanta, GA, United States
  • 3Department of Clinical and Toxicological Analyses, School of Pharmaceutical Sciences, University of São Paulo, São Paulo, Brazil
  • 4Global Virus Network (GVN) Center of Excellence, Australian Infectious Disease Research Centre, Brisbane, QLD, Australia

Introduction: An adult wild-type C57BL/6J mouse model of chikungunya virus (CHIKV) infection and disease has been extensively used to study the alphaviral arthritic immunopathology and to evaluate new interventions. How well mouse models recapitulate the gene expression profiles seen in humans remains controversial.

Methods: Herein we perform a comparative transcriptomics analysis using RNA-Seq datasets from the C57BL/6J CHIKV mouse model with datasets obtained from adults and children acutely infected with CHIKV.

Results: Despite sampling quite different tissues, peripheral blood from humans and feet from mice, gene expression profiles were quite similar, with an overlap of up to ≈50% for up-regulated single copy orthologue differentially expressed genes. Furthermore, high levels of significant concordance between mouse and human were seen for immune pathways and signatures, which were dominated by interferons, T cells and monocyte/macrophages. Importantly, predicted responses to a series of anti-inflammatory drug and biologic treatments also showed cogent similarities between species.

Discussion: Comparative transcriptomics and subsequent pathway analysis provides a detailed picture of how a given model recapitulates human gene expression. Using this method, we show that the C57BL/6J CHIKV mouse model provides a reliable and representative system in which to study CHIKV immunopathology and evaluate new treatments.

Introduction

Chikungunya virus (CHIKV) is a mosquito transmitted alphavirus responsible for sporadic outbreaks of rheumatic disease, the most recent of which (2004-2019) resulted in >10 million cases in >100 countries on four continents (1). The disease is primarily associated with symmetrical polyarthralgia/polyarthritis, with other acute symptoms including fever, rash, myalgia and fatigue. A series of severe atypical manifestations are also recognized, with hospitalization rates ranging from 0.6% to 13% of cases and mortality rate estimates ranging from 0.024% to 0.7%. Chronic disease is also recognized and is primarily characterized by polyarthralgia/polyarthritis, but can also include depression, fatigue and alopecia, with ≈0.3-20% of patients reporting symptoms at 1 year post onset (1). Treatment for arthritogenic alphaviruses usually involves paracetamol (acetaminophen) and/or non-steroidal anti-inflammatory drugs, which can provide symptom relief (2). The rapid accumulation of cases, the occasionally high attack rates (up to 30-75% in certain populations) and the ensuing economic burden, have prompted widespread efforts to develop vaccines (37), antiviral therapies (8, 9) and more effective anti-inflammatory treatments (2).

Development of new interventions usually involves use of mouse models as early evaluation tools. Although a number of animal models have been described (1012), the non-lethal, adult C57BL/6J mouse model of viremia and arthritic foot swelling (1315) has been widely adopted internationally (1623). This model has been used to investigate the virology, immunobiology and immunopathology of CHIKV infections (12, 14, 2430). Given the considerable debate regarding how faithfully mouse models recapitulate the transcriptional responses seen in humans (10, 3138), we recently developed a series of bioinformatic methods for investigating how well the transcriptional responses in a given mouse model mimic those seen in humans (39). For human CHIKV patients, RNA-Seq data is available for peripheral blood from acutely infected adults and children (40, 41). For the adult C57BL/6J mouse model, RNA-Seq data is available for both lymph nodes and arthritic feet during acute infection (14, 15). RNA-Seq provides an excellent tool for detailed characterization of any given mouse model to ascertain how faithfully it recapitulates human disease responses and to determine its reliability for evaluation and development of new interventions. Herein we use bioinformatic methods to evaluate the adult wild-type C57BL/6J mouse model of CHIKV infection and arthritic disease. Despite comparing transcriptional signatures from the peripheral blood of human patients with arthritic feet from infected mice, nearly half the genes significantly up-regulated by infection were shared between species. Pathway analyses also illustrated highly significant concordance for inflammatory responses, and predicted multiple drugs and biologics that have seen evaluation for CHIKV rheumatic disease. By these criteria, this mouse model of CHIKV shows a high level of consensus with human disease.

Results

RNA-Seq data sets for CHIKV infection in mice and humans

The mouse data sets were derived from pooled hind feet and inguinal lymph nodes of CHIKV and mock infected adult female C57BL/6J mice at 2 days post infected (dpi) (peak viremia) and 7 dpi (peak arthritic foot swelling) (15) (Table 1). The raw data (fastq files) were reanalyzed herein using STAR, RSEM and EdgeR, with a q<0.05 filter applied to provide Differentially Expressed Genes (DEGs) (Tables S1A, D, G, J, M, P). For each of these four DEG lists, a mouse-human orthologue DEG list (orthoDEGs) (Tables S1B, E, H, K, N, Q) and a single copy orthologue DEG list (scoDEGs) was generated (Tables S1C, F, I, L, O, R).

TABLE 1
www.frontiersin.org

Table 1 Origins of human and mouse gene expression datasets.

The data sets from infected human adults were obtained from individual whole blood of adult CHIKV patients where blood was collected at or within 2 days of disease onset (n=30), with healthy adults providing control samples (n=20) (41) (Table 1). DEG, orthoDEG and scoDEG lists were generated (Tables S1A–C).

The data sets from human children were obtained from individual whole blood of acute pediatric CHIKV cases, with the controls being convalescent whole blood (40) (Table 1). Multi-Dimensional Scaling (MDS) and Molecular Degree of Perturbation analyses identified 12 and 20 samples, respectively, as outliers (Figure S1), which were removed leaving 140 accessions (n=71 after collapse of technical replicates) for n=35 infected and n=36 convalescent samples (Table 1). DEG, orthoDEG and scoDEG lists were generated (Tables S1D–F).

Single copy orthologues, principle component analysis and hierarchical clustering

When expression data from all single copy orthologues from all samples were displayed on a MDS plot, a clear separation emerged between Adults, Children and Mice samples, with differences between these 3 groups much larger than differences between infected and uninfected samples (Figure 1A). The source of the samples (Adult, Children or Mice), rather than the presence of a CHIKV infection, thus had the dominant role in determining the gene expression profiles.

FIGURE 1
www.frontiersin.org

Figure 1 Principal component and hierarchical cluster analyses, viral loads, DEGs and scoDEGs. (A) Principal component analysis (PCA) plot (PC1 - Dim1, PC2 - Dim2) for all accessions using all single copy orthologues. (B) Hierarchical cluster analysis using log2 TMM-normalized read counts of the top 500 single copy orthologues ranked by PC1 and PC2 loadings. (C) The number of viral reads in each accession expressed as a percentage of the total number of reads aligning to protein coding genes in the human or mouse genome in the same sample. Cross-bars represent the mean for each group. (D) Histogram showing the number of DEGs and scoDEGs for each group after application of the indicated significance cutoffs (q, FDR).

To more fully understand the main sources of variation between the groups, a hierarchical clustering analysis was undertaken (Figure 1B). The largest contribution to variation in mRNA expression patterns between the groups was species (Figure 2B, blue & red), consistent with previous analyses of COVID-19 mouse models (39). Thereafter, the largest contributors to variation were mouse tissue type (feet vs. LNs) (Figure 2B, Tissue, purple & green) and Adults vs. Children (Figure 2B, Dataset, brown & pale blue). Viral load (see below) and days post infection provided comparatively minor contributions to the Euclidian distance, consistent with Figure 1A and previous analyses (39).

FIGURE 2
www.frontiersin.org

Figure 2 ScoDEG overlaps. (A) Upset plot of overlaps of upregulated scoDEGs (from Figure 1D). The overlap between the union scoDEGs (scoDEGs from all mouse groups) and human adults and children. (B) The percent overlap between scoDEGs for all groups. Corresponding Euler diagrams in Figure S4. * Highest overlaps, percent of up-regulated scoDEGs identified in Mice Ft 2 and 7 dpi that were also identified in human groups. # For the same groups, but the percentage of up-regulated scoDEGs identified in humans that were also identified in Mice Ft 2 and 7 dpi. (C) Top two Euler diagrams; scoDEG overlaps using the FDR (q) cutoffs described in Figure 1D, but where all scoDEGs from all mouse groups are concatenated (union scoDEGs). Red boxes; overlap of scoDEGs identified in humans that were also identified in any mouse group. Bottom two Euler diagrams; same as above except using a cutoff of q<0.01 for scoDEG list generation for all groups.

Blood transcription modules and peripheral blood lymphopenia

The sizable separation between Adult and Children (Figure 1A) is perhaps surprising, although CHIKV arthropathy is generally viewed as less severe in children when compared with adults (1, 42). Early innate responses are also reported to be more vigorous in children (43) and dominated by monocyte driven responses (40). Consistent with the latter findings, Gene Set Enrichment Analysis (GSEA) using Blood Transcription Modules (BTMs) (44) and gene lists ranked by fold change, revealed more significant GSEAs for Children than for Adults, often for BTMs associated with monocytes (Figure S2). Significant GSEAs for BTMs associated with antiviral responses were also evident in both gene sets, and importantly, both gene sets showed signatures consistent with lymphopenia (Figure S2). Transient peripheral blood lymphopenia is a known feature of acute CHIKV infections in primates (4547).

Comparable viral loads for mice feet and human groups

The reference genomes for mouse and human were augmented to include the CHIKV viral genome (Reunion Island isolate, LR2006_OPY1, GenBank KT449801), with viral read counts providing quantitation of the viral RNA load within each sample. Viral read counts are presented as a percentage of read counts aligning to mouse or human protein coding genes, with the results illustrating that both mouse and human groups covered an overlapping ≈4 log range (Figure 1C). The mean percent viral read counts for human samples was comparable to those seen in mice feet, with LNs showing 2-3 logs lower viral reads (Figure 1C).

As the latent period for CHIKV in humans is 2-6 days, both human and murine studies represent samples taken within a week of infection and encompass both the period of acute viral replication and the, often abrupt and fulminant, onset of rheumatic disease (13, 40, 41, 48). Onset of disease (usually fever and arthropathy) in humans occurs around the time of peak viremia (49, 50). Blood for both human data sets was obtained within 2 days of symptom onset (40, 41). The lower viral loads in Children (Figure 1C) are likely associated with the faster clearance of viremia as a result of the aforementioned stronger early innate immune responses (43).

DEGs and scoDEGs

The number of biological replicates for the human studies were much higher than for the mice study [which pooled RNA from 4 mice for each of the 3 replicates (15)] (Figures 1A, C). As a result, DEG numbers were substantially higher for the human studies. The q value (FDR) cutoff was thus adjusted so (i) that human and mouse groups had broadly similar number of scoDEGs, allowing human-mouse comparisons without introducing bias from different gene set sizes, and (ii) that DEG numbers were below 8000, the limit for Ingenuity Pathway Analysis (IPA) (see below) (Figure 1D).

When the human DEG lists were curated to contain only scoDEGs, nearly a third of the DEGs were lost (Figure 1D). These lost human non-scoDEGs were analyzed further to ensure that important information was not lost from this process. Nearly 40% of non-scoDEGs had an ENS prefix (Ensembl) rather than a HUGO ID (Figure S3A), and primarily represent long non-coding RNA and pseudogenes. These genes also had, on average, lower read counts (Figure S3B). Finally, IPA analysis of the non-scoDEGs provided very similar UpStream Regulator (USR) z-scores as those identified for scoDEGs (Figure S3C). Thus removal of non-scoDEGs from the human DEG lists would appear not to excessively remove important information.

Overlap between mouse and human scoDEGs

For the up-regulated scoDEGs described in Figure 1D, the scoDEG overlaps between the groups is illustrated in an Upset plot (Figure 2A). For many scoDEGs there was no overlap between mouse and human groups (Figure 2A, black bars and lines/circles). Using a series of Euler diagrams (Figure S4), the best scoDEG overlaps between species was seen for mice feet were 48% (2 dpi) and 42% (7 dpi) of up-regulated scoDEGs were also identified in human blood (Figure 2B, asterisks). However, for the same groups, only 28% and 33% of scoDEGs identified in humans were up-regulated scoDEGs in mice feet for 2 and 7 dpi, respectively (Figure 2B, hash symbols). Overlap percentages were lower for down-regulated scoDEGs, and lower again for LNs (Figures 2B, S4). That down-regulated scoDEGs show lower levels of overlap may reflect infection of different cell types (e.g. fibroblasts in feet (1) vs. monocytes in blood (51)), but likely also reflects CHIKV-induced lymphopenia in human peripheral blood (Figure S2) versus arthritic infiltrates in mice feet (13).

The aforementioned analysis (Figure 2B, hash symbols) argues that the majority of the scoDEGs identified in peripheral blood of humans were not scoDEGs in mice. However, when the union of mouse scoDEGs was compared with humans, percentage overlaps increased, with 49% and 40% of scoDEGs identified in humans seen in at least one mice group, for up and down-regulated scoDEGs, respectively (Figure 2C, red boxes). These percentages did not increase when the FDR (q) values were changed to q<0.01 (for DEG cutoff) for the human groups (Figure 2C, bottom Euler diagrams).

Given the different source of infected material (feet for mice and peripheral blood for humans) (Table 1), overlaps of nearly 50% for up-regulated scoDEGs (Figure 2B, asterisks; Figure 2C, red box) might be considered a relatively high level of consensus between species (39).

Reciprocal GSEAs show species concordance

Reciprocal GSEAs were performed using the gene lists ranked by fold change and the DEG sets for each group. For mouse/human single-copy orthologues with different HUGO IDs between species (≈8%), the mouse gene symbols were changed to their orthologous human equivalent in the DEG sets and the orthologue lists. This allowed GSEAs to be undertaken for mouse vs. human orthologue sets. A separate GSEA was performed for each up- and down-regulated DEG set.

The human up-regulated DEGs were significantly enriched with positive Normalized Enrichment Scores (NES) in 5/8 of the ranked mice gene lists (Figure 3A, blue circles; Table S2). The up-regulated DEGs from mice were significantly enriched with positive NES in 6/8 of the ranked human gene lists (Figure 3A, green circles; Table S2). The numbers were 5/8 (human-mouse) and 1/8 (mouse-human) for down-regulated DEGs and negative NES scores (Figure 3A, blue and green crosses; Table S2). Five significant GSEAs found signatures for up-regulated DEGs in the down-regulated genes or down-regulated DEGs in the up-regulated genes (Figure 3, opposite direction GSEAs); a result that likely reflects, in part, the aforementioned leukopenia. Thus although some incongruities were evident, these GSEA results argue that genes significantly up-regulated in mice, often showed significant enrichment in human ranked gene lists and vice versa.

FIGURE 3
www.frontiersin.org

Figure 3 GSEAs between groups and using ImmuneSigDB gene sets. (A) Reciprocal GSEAs using genes lists ranked by fold change and up- and down-regulated DEG lists from each group. Most GSEAs were significant (left of q=0.05 line), with the number of significant GSEAs over the number of GSEAs provided top right (q<0.05 & same direction). Five GSEAs gave opposite direction; e.g. up-regulated orthoDEG sets significantly enriched in the down-regulated ranked gene list (Opposite direction GSEAs). (B) GSEAs using up-regulated gene sets from ImmuneSigDB and gene lists ranked by fold change for all the groups. GSEAs significant in at least one group are segregated by cell type (mentioned in the ImmuneSigDB annotation) and ranked by Adults. Blue represents up-regulated ImmuneSigDB gene sets that are enriched in the down-regulated genes in the ranked gene lists. (C) For the GSEAs shown in b, NES for Adults vs. Ft 7 dpi and Children vs. Ft 7 dpi were plotted for the indicated cell types and Pearson correlation coefficients and significance values obtained.

GSEAs using ImmuneSigDB gene sets

ImmuneSigDB is a compendium of 5219 immunology-specific gene sets that can be used to interrogate mouse and human ranked gene lists using GSEAs to identify immune signatures (37). A ranked gene list was used as input for each GSEA. As before, mouse/human single-copy orthologues with HUGO IDs that differed between species were changed in the mouse lists to their orthologous human equivalents. GSEAs using up-regulated ImmuneSigDB gene sets that were significantly enriched in at least one ranked gene list were used to compare groups (n=1616) (Table S3).

The NES are plotted for GSEAs that reached significance (q<0.05) in at least one group, with the NES ranked by Adults (Figure 3B). The results were also grouped by cell type according to the specific cell type mentioned in the ImmuneSigDB gene set annotations (Figure 3B) (37, 39). GSEAs with negative NES were prominent in the human groups (Figure 3B, Negative NES, Blue), an observation consistent with the peripheral blood lymphopenia (45, 52) (Figure S2). The inferred lymphopenia signature was less prominent in Children (Figure 3B, Negative NES, Blue), perhaps associated with the lower viral loads (Figure 1C). As might be expected, lymphopenia was not a feature of mouse groups, as arthritic feet and LN were sampled, rather than peripheral blood (1, 5355). This GSEA incongruity (Figure 3B, Negative NES, Blue) is thus more likely due to the source of the samples, rather than representing a difference between species.

For the GSEAs with positive NES scores there was considerably more concordance between mouse and human groups, with the exception of Mice LN 7 dpi (Figure 3B, Red). Mice LN 7 dpi showed the lowest viral loads of all the groups (Figure 1C), with the viremia usually no longer detectable by 6 dpi in this model, although high tissue viral titers remain in feet (13, 15). The larger number of ImmuneSigDB GSEAs with positive NES scores for Mice Ft 7 dpi (Figure 3B, Red), is likely due to the overt foot swelling and pronounced, CD4 T cell-dependent, monocyte/macrophage dominated, inflammatory infiltrate seen in feet at this time (13, 5557).

For the main cell types in Figure 3B, Pearson correlations were determined for the NES scores for Adults vs. Mice Ft 7 dpi and Children vs. Mice Ft 7 dpi (Figure 3C). Correlations were highly significant for the main cell types involved in CHIKV arthritis; T cells (5658), monocyte/macrophages (45, 51, 55) and NK cells (59), and to a lesser extent neutrophils (54, 60), regulatory T cells (61, 62) and dendritic cells (13) (Figure 3C). Coefficients of determination (r2) were generally not so high (Figure 3C), likely reflecting the influence of the aforementioned lymphopenia.

Despite the confounding effect of lymphopenia vs. arthritic infiltrates, overall these analyses argue that during acute CHIKV infections, the dominant immune responses up-regulated in peripheral blood in humans are surprisingly similar to the dominant immune responses up-regulated in arthritic feet in mice.

High species concordance for inflammatory pathways

IPA accepts both human and mouse gene IDs, and the Up-Stream Regulator (USR) feature identifies the most likely up-stream causes (by p values and z-scores) of the transcriptional changes seen in any given DEG list. Using the DEG lists (Figure 1D) significant (p<0.05) IPA USRs were identified for each group (Table S4). USR are grouped by Molecule type, and for Cytokine and Transcription regulator USRs, the top 50 USRs by absolute z-score are shown ranked by Adults (Figure 4). A high level of concordance is apparent across species, with the exception of LN 7 dpi (Figure 4, Cytokine, Transcription regulator). The general congruence between Adults and Children (Figure 4) is consistent with a previous study showing a common pattern of immune signatures associated with CHIKV infections in humans across age groups (43). The contention is further supported by plots showing relatively minor differences in USR z-scores for Adults versus Children (Figure S5).

FIGURE 4
www.frontiersin.org

Figure 4 IPA USRs heatmaps. The DEG lists for all groups were analyses by IPA and the top 50 USRs (significant in at least 1 group) plotted by molecule type and ranked by USR z-scores for Adults. Cytokine and Transcription regulator – USRs for Adults ranked by absolute z-score. Biologic and Chemical drug – ranked from highest negative z-score (most potent predicted inhibitor). Red arrows indicates compounds referred to in the text. Abbreviations: AC acti. pp - adenylate cyclase-activating polypeptide; intest. pep. intestinal peptide; APP der. comp. II - 6-aminopyrazolopyrimidine derivative compound II; h.r.p.- hormone related protein; HAART - highly active antiretroviral therapy; tyr.kin.inh. tyrosine kinase inhibitor.

Most of the top cytokines and transcriptional regulators (Figure 4) are well described in the CHIKV literature, and as expected are dominated by interferon (IFN) and Th1 responses (1, 15, 40, 41, 54, 60, 6366). Perhaps surprising is the absence of CCL2/CCR2 as an USR, even though CCL2 is a highly significant DEG, and is well described for CHIKV infections (45, 55, 67). Conceivably, this is an under-annotation issue for IPA.

Species concordance for biologic and chemical drug USRs

A key question is whether new treatments evaluated in mice might provide reliable insights into therapeutic activity in humans. Biologic and Chemical drug USRs (identified by IPA as above) were ranked by negative NES (most likely to inhibit) in Adults (Figure 4, Biologic drug, Chemical drug). Many of the drugs identified herein have been the subject of independent evaluations in mice and/or humans (Figure 4, red arrows). Only very early subcutaneous edema is reduced in feet of IFNγ-/- mice, with marginally increased viremia (15), with the anti-IFNγ biologic, fontolizumab, to the best of our knowledge not evaluated for any alphavirus. Targeting IL1B with Anakinra shows some benefit in mice (68), but Anakinra has not been evaluated in CHIKV patients. Experience with anti-TNF biologics, such as infliximamb (anti-TNF monoclonal antibody) and etanercept (TNF receptor fusion protein), have been reported in CHIKV patients (69, 70) and may reduce duration of symptoms (71). Mouse experiments are limited to a weanling mouse model of Ross River virus (RRV) where etanercept treatment was started 1 day after infection (i.e. before antibody production), and caused exacerbated disease and lethality (72). The anti-CD80/CD86 monoclonal, abatacept, has shown promise for CHIKV arthritis (that peaks 6/7 dpi) in mice (73), but IPA analysis of human expression data from peripheral blood did not predict efficacy in CHIKV patients (Figure 4, abatacept z-score = 0). Bone loss (not generally a significant feature of CHIKV disease in humans) is reduced by anti-IL-6 antibody in the weanling mouse model of RRV (74). However, limited studies on the anti-IL6 receptor antibody, tocilizumab, have not shown benefit in CHIKV patients (71), nor does the IPA analysis predict efficacy (Figure 4, Biologic drug).

The top scoring Chemical drugs are SB203580 (Adezmapimod, a p38 MAP kinase inhibitor) and LY294002 (PI3K inhibitor) (Figure 4, Chemical drug), with both shown to have some anti-alphaviral activity in vitro (75, 76). Importantly, aspirin, a well annotated representative of non-steroidal anti-inflammatory drugs (NSAIDs) is identified across species, with NSAIDs the mainstay of treatment for alphaviral arthropathies (48, 77). Fingolimod (a sphingosine-1-phosphate receptor modulator that sequesters lymphocytes in lymph nodes) has shown efficacy against CHIKV arthropathy in mice (78). MCC950 is absent, despite activity in mice (27), as no annotation exists in IPA for this NLRP3 inhibitor. Also absent are disease modifying anti-rheumatic drugs (DMARDs) (e.g. methotrexate and sulfasalazine), with these well annotated within IPA; however, despite a number of studies, no clear benefit has been established in patients with acute CHIKV (77). Failure to formally exclude CHIKV patients that may also have underlying autoimmune arthritides may have complicated a number of DMARD studies (1). Corticosteroid use (79, 80) was similarly not highlighted by our analyses.

Taken together these analyses point to a high level of conservation for predicted drug sensitivities between CHIKV patients and the adult C57BL/6J CHIKV mouse model, with many of the drugs that were identified (Figure 4, red arrows) also the subject of investigations in mice and CHIKV patients.

Human USR pathways correlate better with mouse feet than mouse LNs

Pearson correlations were undertaken for pairwise comparisons of z-scores across all USR molecule types for Adults and Children versus mouse groups. Coefficients of determination (r2) for human versus mice feet (Figure 5A; Figure S6) were nearly always higher than for human versus mice LNs (Figure 5B), with LN 7 dpi again showing very poor correlations (Figure 5B, LN 7 dpi). Thus for acute CHIKV infections, arthritic feet from mice, rather than LNs, show a higher level of congruence with human peripheral blood across multiple pathways. This may reflect the fact that leukocytes, including antigen-specific T cells and other lymphocytes, extravasate from the peripheral blood to the sites of infection (joints) (81) and become responsible for the CHIKV inflammatory arthropathy (1).

FIGURE 5
www.frontiersin.org

Figure 5 Pearson correlations for IPA USRs. (A) The DEGs from Adults, Children, mice Ft 2 dpi, and mice Ft 7 dpi were analyzed by IPA to provide z-scores for significant USRs (p<0.05). The USRs were then grouped by Molecule type and Pearson correlations undertaken. For example, top left, Cytokine, compares Adults Cytokine USR z-scores with mice feet 2 dpi and 7 dpi Cytokine USR z-scores (correlation plot shown in Figure S6A). (B) As for a, but correlations with mice LN 2 and 7 dpi.

Specificity of CHIKV signatures

A potential criticism of this type of interspecies comparison is that the signatures being compared are simply a generic reflection of infection and inflammation. Therefore, by reanalyzing publically available RNA-Seq data, we compared the degree of correlation of pathway signatures between CHIKV-infected humans and CHIKV-infected mice with those of CHIKV-infected humans and two non-CHIKV sources of human inflammatory disease. Pearson coefficients of determination comparing Cytokine USR z-scores from (i) PBMCs of adults with acute CHIKV infection versus feet of CHIKV-infected mice, were considerably higher than (ii) PBMCs of adults with acute CHIKV infection versus PBMCs of children with acute bronchiolitis (82) (Figure 6; Cytokine). A similar result was obtained for comparisons of Cytokine USR z-scores from PBMCs of adults with acute CHIKV infection versus PBMCs of adults acutely infected with SARS-CoV-2 (83) (Figure 6). Using Cytokine USR z-scores from peripheral blood of Children with acute CHIKV infections instead of Adults, did not significantly change these results (Figure S7A). Thus the best determinant of congruence of cytokine signatures was CHIKV, rather than infection, species or tissue.

FIGURE 6
www.frontiersin.org

Figure 6 Pearson correlations for IPA USRs with other infections and rheumatoid arthritis. The DEGs from whole blood of CHIKV-infected Adults were analyzed by IPA to provide z-scores for significant USRs and USRs for (Molecule type) Cytokine, Biologic drug and Chemical drug. These z-scores were compared by Pearson correlations with z-scores from peripheral blood mononuclear cells (PBMC) from children with bronchiolitis and Adults infected with SARS-CoV-2, as well as synovial biopsies from rheumatoid arthritis (RA) patients. Coefficients of determination (r2) are plotted.

Given our previous work on the similarities between rheumatoid arthritis (RA) in humans and CHIKV arthritis in mice (based on microarray studies) (56), we accessed publically available RNA-Seq data from synovial biopsies of RA patients (compared with healthy controls), generated DEGs, obtained IPA USR z-scores and used these for the same correlations. Pearson coefficients of determination were higher than those seen for bronchiolitis and SARS-CoV-2, although did not reach the r2 values seen for correlations with mice groups (Figure 6). These results support the view that arthritic signatures are observable in peripheral blood (84), and reaffirm the view (also seen in Figure 4, drug) that drugs used in the treatment of RA may find utility in the treatment of alphaviral arthritides (56) (Figure S7B).

C57BL/6J versus C57BL/6N mice

The C57BL/6J mouse, bred at the Jackson Laboratory, is arguably the most commonly used mouse strain in medical research. There are, however, multiple C57BL/6 mouse strains (85), including the original C57BL/6N mice from which many knock-out mice are derived. The latter have a number of genetic differences from C57BL/6J mice (86, 87) and show an ameliorated CHIKV arthropathy, in part due to the presence of an intact Nicotinamide Nucleotide Transhydrogenase (Nnt) gene (14). MDS plots illustrate all mice infected feet groups cluster together closely when Adults and Children are included, illustrating that differences in transcriptional profiles between Adults and Children is much greater than between mouse strains (Figure S8A). When Children are removed, feet from infected C57BL/6N mice clearly segregate from C57BL/6J mice, although the Dim2 axis (on which segregation is evident) accounts for only 7% of the total variance (Figure S8B). Only when C57BL/6N are analyzed by themselves can the influence of Nnt be seen on the MDS plot (Figure S8C); the influence was explained in greater detail previously (14). This expression profiling argues that C57BL/6N mice do not ostensibly provide a better model of CHIKV, when compared with C57BL/6J mice, with overall differences between mouse strains actually quite minor when compared with differences between humans.

Discussion

Herein we provide a series of bioinformatic approaches that illustrate a validation process for the adult C57BL/6J mouse model of CHIKV infection and disease using publically available RNA-Seq data sets. Overlap of up-regulated scoDEGs reached nearly 50%, despite the fact that human data sets were derived from peripheral blood (with associated lymphopenia) and murine data sets were obtained from arthritic feet (that contain inflammatory infiltrates). This scoDEG overlap is higher than the ≈15-35% seen for a similar bioinformatic study for SARS-CoV-2 infections in humans and mouse lung tissues (39). For the latter study, only fixed human lung tissues were available and days post disease onset for the COVID-19 patients was also not known, perhaps contributing to a lower level of overlap. Either way, both studies showed a high level of congruence for inflammatory and immune pathways between mice and humans. Importantly, the predicted behavior of a range of anti-inflammatory CHIKV treatments also showed a high level of concordance between humans and mice, arguing that the CHIKV mouse model represents a reliable and representative model in which to evaluate anti-inflammatory interventions.

A limitation of this study is that we were unable to compare mouse and human studies of chronic CHIKV disease. Chronic CHIKV disease (primarily arthropathy) is well described in humans (88, 89); however, RNA-Seq data is not available and establishing cohorts of chronic CHIKV patients free of other underlying (often autoimmune) rheumatic conditions is often difficult (1, 90). RNA-Seq data for 30 dpi is available for the C57BL/6J mouse model, but this might be viewed as post-acute rather than chronic (91), and largely only portrays a diminution of acute responses (15). We were also unable to address severe manifestations of CHIKV, which can include lethality, usually in very young and elderly patients (1, 92). Although lethal mouse models of CHIKV exist, they involve use of young mice or GMO mice defective in the type I IFN pathway (64, 93, 94). RNA-Seq data is currently not publically available for such mouse models, nor for CHIKV patients with lethal outcomes. Another limitation is that current bioinformatic approaches are unable reliably to predict treatments that might exacerbate inflammation (54, 55) or compromise anti-viral immunity (47, 72, 95, 96). However, with respect to the latter, in most settings, diagnosis of alphaviral arthritis involves paired serology (97), with treatment thus generally initiated after protective anti-viral antibodies have been generated. Finally, there are many biological processes that are cell-type specific and/or for which pathway annotations are poor or not readily available. In such cases the whole tissue RNA-Seq approach and the pathway analysis tools used herein may often be uninformative. An example would be the interplay between apoptosis and autophagy (98, 99), where species-specific response have been reported (100).

Understanding how a mouse model does or does not recapitulate certain responses seen during human disease is often important for assessing the validity of mouse experiments for understanding aspects of infection and disease, and for evaluation of particular interventions. Given there are published reports of mouse models recapitulating responses to infection and inflammation both well (33, 39, 101) and poorly (10, 34, 102), RNA-Seq and subsequent pathway analyses provide a useful strategy for determining exactly how well any given model recapitulates specific gene expression patterns seen in human patients. The adult wild-type C57BL/6J mouse model of CHIKV shows up-regulation of many genes that are up-regulated in humans, as well as showing a high level of concordance for immune and inflammation pathways, arguing that it overall provides a highly representative model of human disease.

Materials and methods

Gene expression analysis

Raw sequence reads (see Table 1) were accessed from the National Centre for Biotechnology Information Sequence Read Archive using Aspera (IBM). Quality control of sequence reads was performed using FastQC v0.11.9 (103). Reads were quality-trimmed with a minimum Phred score cutoff of Q20, and size-selected with a minimum length of 36 nucleotides, using Trimmomatic v0.36 (104). Trimmed reads were aligned to either the GRCm39 vM26 or GRCh38 v37 reference genome for mouse and human datasets, respectively, using STAR aligner v2.7.1a (105). The KT449801.1 CHIKV genome (58) was appended to each reference genome prior to read-alignment to allow quantitation of viral reads. Quantitation of viral reads was performed using primary proper paired reads, with Samtools v1.9 (106). Gene expression was calculated using RSEM v1.3.1 (107). Differentially expressed genes were identified using EdgeR v3.34.0 (108) in R v4.1.0 (109). Mouse-human orthologues and single-copy mouse-human orthologues were obtained from the Ensembl database using BiomaRt v2.48.2 (110) in R. A proportion of mouse-human single-copy orthologues (~8%) have HUGO IDs that differ between species. Therefore, all mouse samples had the HUGO IDs of these genes converted to the human equivalent prior to performing cross-species comparisons. Gene expression profiles were visualized across all samples using multi-dimensional scaling and hierarchical clustering analysis in R. Data were plotted using R packages, Eulerr v6.1.0 (111), ggplot2 v3.3.5 (112), ggpubR v0.4.0 (113), ComplexUpset v1.3.3 (114), PCAtools v2.4.0 (115), and RColorBrewer v1.1 (116).

Reciprocal gene set enrichment analysis

An up- and a down-regulated DEG set was created for each group. A minimum absolute log2 fold-change cutoff was determined separately for each group by ranking all DEGs by absolute log2 fold-change (i.e. ignoring directionality) and then calculating the 50th percentile. For a DEG to be included in the set, it had to have either a log2 fold-change greater than the cutoff (up-regulated DEGs) or lower than the negative cutoff (down-regulated DEGs). A log2 fold-change ranked gene list was produced for each group using EdgeR as described above. Pairwise Gene Set Enrichment analyses were performed using the up- and down-regulated DEG sets and ranked gene lists, with GSEA v4.1.0 (117) with 100 permutations and the ‘no_collapse’ setting. HUGO IDs of single-copy orthologues were standardized between species as described above.

ImmuneSigDB gene set enrichment analysis

The ImmuneSigDB v7.4 gene set collection composed of 5219 immune-related gene sets was obtained from the Molecular Signatures Database. For each dataset and time point a GSEA was performed to test for enrichment of up-regulated ImmuneSigDB gene sets within the log2 fold-change ranked gene lists produced using EdgeR as described above. Data were visualized using pheatmap v1.0.12 (118) in R.

Pathway analysis

Pathway analysis was performed using Ingenuity Pathway Analysis (IPA) v65367011 (Qiagen) with default settings, using log2 fold-changes of DEGs produced using EdgeR as described above. Data were visualized using pheatmap v1.0.12 in R.

Comparison with non-CHIKV disease

The NCBI Bioprojects database was searched for projects containing the terms RNA-Seq, Homo sapiens, inflammation, SARS-CoV-2, PBMC, rheumatoid arthritis, and synovial. From the results of the search, three projects were chosen according to the following criteria: they were each associated with at least one peer-reviewed publication, had adequate controls (i.e. uninfected healthy individuals), and were generally similar to the CHIKV projects with respect to tissue, sequencing platform and method (i.e. bulk RNA-Seq). For the bronchiolitis dataset, all child samples were retained, and all infant samples were removed. For the SARS-CoV-2 dataset, all samples from healthy patients, and all samples from infected patients at zero days post testing PCR-positive, were retained. For the rheumatoid arthritis dataset, 20 age matched healthy and ‘established’ rheumatoid arthritis samples were retained. Differential expression was measured as described above, with a q-value cutoff of 0.05. Pathway analysis was performed using IPA as described above.

Statistics

Statistics were performed using R version 4.1.0. For gene expression and pathway enrichment data a Pearson’s correlation test was used in accordance with the central limit theorem. Correlations with a p-value < 0.05 were considered significant.

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

Author contributions

CB, FC undertook the bioinformatic analyses. HN and AS supervised the research and obtained the funding. AS and CB wrote the manuscript with input from all the authors. All authors contributed to the article and approved the submitted version.

Funding

The work was supported by an Investigator grant from the National Health and Medical Research Council (NHMRC) of Australia (APP1173880) awarded to AS. HN and FC received support by Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP) (2018/14933-2), Brazil.

Acknowledgments

The authors would like to thank Dr Gunter Hartel (Head of Statistics, QIMR Berghofer Institute of Medical Research) for help with statistics and Dr Viviana Lutzky (QIMR Berghofer MRI) for manuscript review and proof reading.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2022.1092370/full#supplementary-material

Supplementary Table 1 | DEGs, orthoDEGs, and scoDEGs from mice and humans in response to CHIKV infection.

Supplementary Table 2 | Reciprocal Gene Set Enrichment Analysis.

Supplementary Table 3 | ImmuneSigDB Gene Set Enrichment Analysis.

Supplementary Table 4 | Upstream Regulator annotations enriched in DEG lists.

Supplementary Table 5 | Information relating to three non-CHIKV human inflammatory disease datasets

References

1. Suhrbier A. Rheumatic manifestations of chikungunya: emerging concepts and interventions. Nat Rev Rheumatol (2019) 15(10):597–611. doi: 10.1038/s41584-019-0276-9

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Rodrigo C, Herath T, Wickramarachchi U, Fernando D, Rajapakse S. Treatment of chikungunya-associated joint pain: a systematic review of controlled clinical trials. Trans R Soc Trop Med Hyg (2022) 116(10):889–99. doi: 10.1093/trstmh/trac045

PubMed Abstract | CrossRef Full Text | Google Scholar

3. de Lima Cavalcanti TYV, Pereira MR, de Paula SO, Franca RFO. A review on chikungunya virus epidemiology, pathogenesis and current vaccine development. Viruses (2022) 14(5):969. doi: 10.3390/v14050969

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Bettis AA, L'Azou Jackson M, Yoon IK, Breugelmans JG, Goios A, Gubler DJ, et al. The global epidemiology of chikungunya from 1999 to 2020: A systematic literature review to inform the development and introduction of vaccines. PloS Negl Trop Dis (2022) 16(1):e0010069. doi: 10.1371/journal.pntd.0010069

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Wressnigg N, Hochreiter R, Zoihsl O, Fritzer A, Bezay N, Klingler A, et al. Single-shot live-attenuated chikungunya vaccine in healthy adults: a phase 1, randomised controlled trial. Lancet Infect Dis (2020) 20(10):1193–203. doi: 10.1016/S1473-3099(20)30238-3

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Chen GL, Coates EE, Plummer SH, Carter CA, Berkowitz N, Conan-Cibotti M, et al. Effect of a chikungunya virus-like particle vaccine on safety and tolerability outcomes: A randomized clinical trial. JAMA (2020) 323(14):1369–77. doi: 10.1001/jama.2020.2477

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Folegatti PM, Harrison K, Preciado-Llanes L, Lopez FR, Bittaye M, Kim YC, et al. A single dose of ChAdOx1 chik vaccine induces neutralizing antibodies against four chikungunya virus lineages in a phase 1 clinical trial. Nat Commun (2021) 12(1):4636. doi: 10.1038/s41467-021-24906-y

PubMed Abstract | CrossRef Full Text | Google Scholar

8. August A, Attarwala HZ, Himansu S, Kalidindi S, Lu S, Pajon R, et al. A phase 1 trial of lipid-encapsulated mRNA encoding a monoclonal antibody with neutralizing activity against chikungunya virus. Nat Med (2021) 27(12):2224–33. doi: 10.1038/s41591-021-01573-6

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Kovacikova K, van Hemert MJ. Small-molecule inhibitors of chikungunya virus: Mechanisms of action and antiviral drug resistance. Antimicrob Agents Chemother (2020) 64(12):e01788–20. doi: 10.1128/AAC.01788-20

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Chan YH, Lum FM, Ng LFP. Limitations of current in vivo mouse models for the study of chikungunya virus pathogenesis. Med Sci (Basel) (2015) 3(3):64–77. doi: 10.3390/medsci3030064

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Haese NN, Broeckel RM, Hawman DW, Heise MT, Morrison TE, Streblow DN. Animal models of chikungunya virus infection and disease. J Infect Dis (2016) 214(suppl 5):S482–S7. doi: 10.1093/infdis/jiw284

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Nguyen W, Nakayama E, Yan K, Tang B, Le TT, Liu L, et al. Arthritogenic alphavirus vaccines: Serogrouping versus cross-protection in mouse models. Vaccines (Basel) (2020) 8(2):209. doi: 10.3390/vaccines8020209

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Gardner J, Anraku I, Le TT, Larcher T, Major L, Roques P, et al. Chikungunya virus arthritis in adult wild-type mice. J Virol (2010) 84(16):8021–32. doi: 10.1128/jvi.02603-09

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Rawle DJ, Le TT, Dumenil T, Bishop C, Yan K, Nakayama E, et al. Widespread discrepancy in nnt genotypes and genetic backgrounds complicates granzyme a and other knockout mouse studies. Elife (2022) 11:e70207. doi: 10.7554/eLife.70207

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Wilson JA, Prow NA, Schroder WA, Ellis JJ, Cumming HE, Gearing LJ, et al. RNA-Seq analysis of chikungunya virus infection and identification of granzyme a as a major promoter of arthritic inflammation. PloS Pathog (2017) 13(2):e1006155. doi: 10.1371/journal.ppat.1006155

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Hallengard D, Kakoulidou M, Lulla A, Kummerer BM, Johansson DX, Mutso M, et al. Novel attenuated chikungunya vaccine candidates elicit protective immunity in C57BL/6 mice. J Virol (2014) 88(5):2858–66. doi: 10.1128/JVI.03453-13

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Garcia-Arriaza J, Cepeda V, Hallengard D, Sorzano CO, Kummerer BM, Liljestrom P, et al. A novel poxvirus-based vaccine, MVA-CHIKV, is highly immunogenic and protects mice against chikungunya infection. J Virol (2014) 88(6):3527–47. doi: 10.1128/JVI.03418-13

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Prow NA, Liu L, Nakayama E, Cooper TH, Yan K, Eldi P, et al. A vaccinia-based single vector construct multi-pathogen vaccine protects against both zika and chikungunya viruses. Nat Commun (2018) 9(1):1230. doi: 10.1038/s41467-018-03662-6

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Zhang YN, Zhang ZR, Li N, Pei XR, Li XD, Deng CL, et al. High-titer self-propagating capsidless chikungunya virus generated in vero cells as a strategy for alphavirus vaccine development. J Virol (2022) 96(6):e0148021. doi: 10.1128/JVI.01480-21

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Julander JG, Anderson N, Haese N, Andoh T, Streblow DN, Cortez P, et al. Therapeutic and prophylactic treatment with a virus-specific antibody is highly effective in rodent models of chikungunya infection and disease. Antiviral Res (2022) 202:105295. doi: 10.1016/j.antiviral.2022.105295

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Jeengar MK, Kurakula M, Patil P, More A, Sistla R, Parashar D. Antiviral activity of stearylamine against chikungunya virus. Chem Phys Lipids (2021) 235:105049. doi: 10.1016/j.chemphyslip.2021.105049

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Slifka DK, Raue HP, Weber WC, Andoh TF, Kreklywich CN, DeFilippis VR, et al. Development of a next-generation chikungunya virus vaccine based on the HydroVax platform. PloS Pathog (2022) 18(7):e1010695. doi: 10.1371/journal.ppat.1010695

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Voigt EA, Fuerte-Stone J, Granger B, Archer J, Van Hoeven N. Live-attenuated RNA hybrid vaccine technology provides single-dose protection against chikungunya virus. Mol Ther (2021) 29(9):2782–93. doi: 10.1016/j.ymthe.2021.05.018

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Rangel Margarita V, McAllister N, Dancel-Manning K, Noval Maria G, Silva Laurie A, Stapleford Kenneth A. Emerging chikungunya virus variants at the E1-E1 interglycoprotein spike interface impact virus attachment and inflammation. J Virol (2022) 96(4):e01586–21. doi: 10.1128/jvi.01586-21

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Fumagalli MJ, de Souza WM, de Castro-Jorge LA, de Carvalho RVH, Castro IA, de Almeida LGN, et al. Chikungunya virus exposure partially cross-protects against mayaro virus infection in mice. J Virol (2021) 95(23):e0112221. doi: 10.1128/JVI.01122-21

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Geng T, Lin T, Yang D, Harrison AG, Vella AT, Fikrig E, et al. A critical role for STING signaling in limiting pathogenesis of chikungunya virus. J Infect Dis (2021) 223(12):2186–96. doi: 10.1093/infdis/jiaa694

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Chen W, Foo SS, Zaid A, Teng TS, Herrero LJ, Wolf S, et al. Specific inhibition of NLRP3 in chikungunya disease reveals a role for inflammasomes in alphavirus-induced inflammation. Nat Microbiol (2017) 2(10):1435–45. doi: 10.1038/s41564-017-0015-4

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Ng LFP. Immunopathology of chikungunya virus infection: Lessons learned from patients and animal models. Annu Rev Virol (2017) 4(1):413–27. doi: 10.1146/annurev-virology-101416-041808

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Zhang Y, Yan H, Li X, Zhou D, Zhong M, Yang J, et al. A high-dose inoculum size results in persistent viral infection and arthritis in mice infected with chikungunya virus. PloS Negl Trop Dis (2022) 16(1):e0010149. doi: 10.1371/journal.pntd.0010149

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Rawle DJ, Dumenil T, Tang B, Bishop CR, Yan K, Le TT, et al. Microplastic consumption induces inflammatory signatures in the colon and prolongs a viral arthritis. Sci Total Environ (2022) 809:152212. doi: 10.1016/j.scitotenv.2021.152212

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Yue F, Cheng Y, Breschi A, Vierstra J, Wu W, Ryba T, et al. A comparative encyclopedia of DNA elements in the mouse genome. Nature (2014) 515(7527):355–64. doi: 10.1038/nature13992

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Lin S, Lin Y, Nery JR, Urich MA, Breschi A, Davis CA, et al. Comparison of the transcriptional landscapes between human and mouse tissues. Proc Natl Acad Sci U.S.A. (2014) 111(48):17224–9. doi: 10.1073/pnas.1413624111

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Takao K, Miyakawa T. Genomic responses in mouse models greatly mimic human inflammatory diseases. Proc Natl Acad Sci U.S.A. (2015) 112(4):1167–72. doi: 10.1073/pnas.1401965111

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Seok J, Warren HS, Cuenca AG, Mindrinos MN, Baker HV, Xu W, et al. Genomic responses in mouse models poorly mimic human inflammatory diseases. Proc Natl Acad Sci U.S.A. (2013) 110(9):3507–12. doi: 10.1073/pnas.1222878110

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Seeley RJ, MacDougald OA. Mice as experimental models for human physiology: when several degrees in housing temperature matter. Nat Metab (2021) 3(4):443–5. doi: 10.1038/s42255-021-00372-0

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Monaco G, van Dam S, Casal Novo Ribeiro JL, Larbi A, de Magalhaes JP. A comparison of human and mouse gene co-expression networks reveals conservation and divergence at the tissue, pathway and disease levels. BMC Evol Biol (2015) 15:259. doi: 10.1186/s12862-015-0534-7

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Godec J, Tan Y, Liberzon A, Tamayo P, Bhattacharya S, Butte AJ, et al. Compendium of immune signatures identifies conserved and species-specific biology in response to inflammation. Immunity (2016) 44(1):194–206. doi: 10.1016/j.immuni.2015.12.006

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Wolf J, Schlecht A, Rosmus DD, Boneva S, Agostini H, Schlunck G, et al. Comparative transcriptome analysis of human and murine choroidal neovascularization identifies fibroblast growth factor inducible-14 as phylogenetically conserved mediator of neovascular age-related macular degeneration. Biochim Biophys Acta Mol Basis Dis (2022) 1868(4):166340. doi: 10.1016/j.bbadis.2022.166340

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Bishop CR, Dumenil T, Rawle DJ, Le TT, Yan K, Tang B, et al. Mouse models of COVID-19 recapitulate inflammatory pathways rather than gene expression. Ploth Pathog (2022) 18(9):e1010867. doi: 10.1101/2022.02.24.481866

CrossRef Full Text | Google Scholar

40. Michlmayr D, Pak TR, Rahman AH, Amir ED, Kim EY, Kim-Schulze S, et al. Comprehensive innate immune profiling of chikungunya virus infection in pediatric cases. Mol Syst Biol (2018) 14(8):e7862. doi: 10.15252/msb.20177862

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Soares-Schanoski A, Baptista Cruz N, de Castro-Jorge LA, de Carvalho RVH, Santos CAD, Ros ND, et al. Systems analysis of subjects acutely infected with the chikungunya virus. PloS Pathog (2019) 15(6):e1007880. doi: 10.1371/journal.ppat.1007880

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Kumar R, Ahmed S, Parray HA, Das S. Chikungunya and arthritis: An overview. Travel Med Infect Dis (2021) 44:102168. doi: 10.1016/j.tmaid.2021.102168

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Simarmata D, Ng DC, Kam YW, Lee B, Sum MS, Her Z, et al. Early clearance of chikungunya virus in children is associated with a strong innate immune response. Sci Rep (2016) 6:26097. doi: 10.1038/srep26097

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Li S, Rouphael N, Duraisingham S, Romero-Steiner S, Presnell S, Davis C, et al. Molecular signatures of antibody responses derived from a systems biology study of five human vaccines. Nat Immunol (2014) 15(2):195–204. doi: 10.1038/ni.2789

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Labadie K, Larcher T, Joubert C, Mannioui A, Delache B, Brochard P, et al. Chikungunya disease in nonhuman primates involves long-term viral persistence in macrophages. J Clin Invest (2010) 120(3):894–906. doi: 10.1172/JCI40104

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Beddingfield BJ, Sugimoto C, Wang E, Weaver SC, Russell-Lodrigue KE, Killeen SZ, et al. Phenotypic and kinetic changes of myeloid lineage cells in innate response to chikungunya infection in cynomolgus macaques. Viral Immunol (2022) 35(3):192–9. doi: 10.1089/vim.2021.0171

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Roques P, Thiberville SD, Dupuis-Maguiraga L, Lum FM, Labadie K, Martinon F, et al. Paradoxical effect of chloroquine treatment in enhancing chikungunya virus infection. Viruses (2018) 10(5):268. doi: 10.3390/v10050268

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Suhrbier A, Jaffar-Bandjee MC, Gasque P. Arthritogenic alphaviruses-an overview. Nat Rev Rheumatol (2012) 8(7):420–9. doi: 10.1038/nrrheum.2012.64

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Chusri S, Siripaitoon P, Silpapojakul K, Hortiwakul T, Charernmak B, Chinnawirotpisan P, et al. Kinetics of chikungunya infections during an outbreak in southern Thailand, 2008-2009. Am J Trop Med Hyg (2014) 90(3):410–7. doi: 10.4269/ajtmh.12-0681

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Riswari SF, Ma'roef CN, Djauhari H, Kosasih H, Perkasa A, Yudhaputri FA, et al. Study of viremic profile in febrile specimens of chikungunya in bandung, Indonesia. J Clin Virol (2016) 74:61–5. doi: 10.1016/j.jcv.2015.11.017

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Her Z, Malleret B, Chan M, Ong EK, Wong SC, Kwek DJ, et al. Active infection of human blood monocytes by chikungunya virus triggers an innate immune response. J Immunol (2010) 184(10):5903–13. doi: 10.4049/jimmunol.0904181

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Ng KW, Chow A, Win MK, Dimatatac F, Neo HY, Lye DC, et al. Clinical features and epidemiology of chikungunya infection in Singapore. Singapore Med J (2009) 50(8):785–90.

PubMed Abstract | Google Scholar

53. Suhrbier A, Mahalingam S. The immunobiology of viral arthritides. Pharmacol Ther (2009) 124(3):301–8. doi: 10.1016/j.pharmthera.2009.09.005

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Prow NA, Hirata TDC, Tang B, Larcher T, Mukhopadhyay P, Alves TL, et al. Exacerbation of chikungunya virus rheumatic immunopathology by a high fiber diet and butyrate. Front Immunol (2019) 10:2736. doi: 10.3389/fimmu.2019.02736

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Poo YS, Nakaya H, Gardner J, Larcher T, Schroder WA, Le TT, et al. CCR2 deficiency promotes exacerbated chronic erosive neutrophil-dominated chikungunya virus arthritis. J Virol (2014) 88(12):6862–72. doi: 10.1128/JVI.03364-13

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Nakaya HI, Gardner J, Poo YS, Major L, Pulendran B, Suhrbier A. Gene profiling of chikungunya virus arthritis in a mouse model reveals significant overlap with rheumatoid arthritis. Arthritis Rheum (2012) 64(11):3553–63. doi: 10.1002/art.34631

PubMed Abstract | CrossRef Full Text | Google Scholar

57. Teo TH, Lum FM, Claser C, Lulla V, Lulla A, Merits A, et al. A pathogenic role for CD4+ T cells during chikungunya virus infection in mice. J Immunol (2013) 190(1):259–69. doi: 10.4049/jimmunol.1202177

PubMed Abstract | CrossRef Full Text | Google Scholar

58. Poo YS, Rudd PA, Gardner J, Wilson JA, Larcher T, Colle MA, et al. Multiple immune factors are involved in controlling acute and chronic chikungunya virus infection. PloS Negl Trop Dis (2014) 8(12):e3354. doi: 10.1371/journal.pntd.0003354

PubMed Abstract | CrossRef Full Text | Google Scholar

59. Teo TH, Her Z, Tan JJ, Lum FM, Lee WW, Chan YH, et al. Caribbean And la reunion chikungunya virus isolates differ in their capacity to induce proinflammatory Th1 and NK cell responses and acute joint pathology. J Virol (2015) 89(15):7955–69. doi: 10.1128/JVI.00909-15

PubMed Abstract | CrossRef Full Text | Google Scholar

60. Liu X, Poo YS, Alves JC, Almeida RP, Mostafavi H, Tang PCH, et al. Interleukin-17 contributes to chikungunya virus-induced disease. mBio (2022) 13(2):e0028922. doi: 10.1128/mbio.00289-22

PubMed Abstract | CrossRef Full Text | Google Scholar

61. Lee WW, Teo TH, Her Z, Lum FM, Kam YW, Haase D, et al. Expanding regulatory T cells alleviates chikungunya virus-induced pathology in mice. J Virol (2015) 89(15):7893–904. doi: 10.1128/JVI.00998-15

PubMed Abstract | CrossRef Full Text | Google Scholar

62. Gois BM, Peixoto RF, Guerra-Gomes IC, Palmeira PHS, Dias CNS, Araujo JMG, et al. Regulatory T cells in acute and chronic human chikungunya infection. Microbes Infect (2022) 24(3):104927. doi: 10.1016/j.micinf.2021.104927

PubMed Abstract | CrossRef Full Text | Google Scholar

63. Alves de Souza TM, Fernandes-Santos C, Araujo da Paixao de Oliveira J, Tome LCT, Fiestas-Solorzano VE, Nunes PCG, et al. Increased indoleamine 2,3-dioxygenase 1 (IDO-1) activity and inflammatory responses during chikungunya virus infection. Pathogens (2022) 11(4):444. doi: 10.3390/pathogens11040444

PubMed Abstract | CrossRef Full Text | Google Scholar

64. Rudd PA, Wilson J, Gardner J, Larcher T, Babarit C, Le TT, et al. Interferon response factors 3 and 7 protect against chikungunya virus hemorrhagic fever and shock. J Virol (2012) 86(18):9888–98. doi: 10.1128/JVI.00956-12

PubMed Abstract | CrossRef Full Text | Google Scholar

65. Schilte C, Couderc T, Chretien F, Sourisseau M, Gangneux N, Guivel-Benhassine F, et al. Type I IFN controls chikungunya virus via its action on nonhematopoietic cells. J Exp Med (2010) 207(2):429–42. doi: 10.1084/jem.20090851

PubMed Abstract | CrossRef Full Text | Google Scholar

66. Torres-Ruesta A, Teo TH, Chan YH, Renia L, Ng LFP. Pathogenic Th1 responses in CHIKV-induced inflammation and their modulation upon plasmodium parasites co-infection. Immunol Rev (2020) 294(1):80–91. doi: 10.1111/imr.12825

PubMed Abstract | CrossRef Full Text | Google Scholar

67. Ruiz Silva M, van der Ende-Metselaar H, Mulder HL, Smit JM, Rodenhuis-Zybert IA. Mechanism and role of MCP-1 upregulation upon chikungunya virus infection in human peripheral blood mononuclear cells. Sci Rep (2016) 6(1):32288. doi: 10.1038/srep32288

PubMed Abstract | CrossRef Full Text | Google Scholar

68. Wolf S, Taylor A, Zaid A, Freitas J, Herrero LJ, Rao S, et al. Inhibition of interleukin-1beta signaling by anakinra demonstrates a critical role of bone loss in experimental arthritogenic alphavirus infections. Arthritis Rheumatol (2019) 71(7):1185–90. doi: 10.1002/art.40856

PubMed Abstract | CrossRef Full Text | Google Scholar

69. Guaraldo L, Wakimoto MD, Ferreira H, Bressan C, Calvet GA, Pinheiro GC, et al. Treatment of chikungunya musculoskeletal disorders: a systematic review. Expert Rev Anti Infect Ther (2018) 16(4):333–44. doi: 10.1080/14787210.2018.1450629

PubMed Abstract | CrossRef Full Text | Google Scholar

70. Bautista-Vargas M, Puerta-Sarmiento G, Canas CA. Characteristics of chikungunya virus infection in patients with established rheumatoid arthritis. Clin Rheumatol (2020) 39(12):3639–42. doi: 10.1007/s10067-020-05198-x

PubMed Abstract | CrossRef Full Text | Google Scholar

71. de Brito CAA, Marques CDL, Franca RFO, Monteiro JR, de Brito MCM, Lopes F, et al. Reduced duration of postchikungunya musculoskeletal pain in rheumatological patients treated with biologicals. J Trop Med (2020) 2020:2071325. doi: 10.1155/2020/2071325

PubMed Abstract | CrossRef Full Text | Google Scholar

72. Zaid A, Rulli NE, Rolph MS, Suhrbier A, Mahalingam S. Disease exacerbation by etanercept in a mouse model of alphaviral arthritis and myositis. Arthritis Rheum (2011) 63(2):488–91. doi: 10.1002/art.30112

PubMed Abstract | CrossRef Full Text | Google Scholar

73. Miner JJ, Cook LE, Hong JP, Smith AM, Richner JM, Shimak RM, et al. Therapy with CTLA4-ig and an antiviral monoclonal antibody controls chikungunya virus arthritis. Sci Transl Med (2017) 9(375):eaah3438. doi: 10.1126/scitranslmed.aah3438

PubMed Abstract | CrossRef Full Text | Google Scholar

74. Chen W, Foo SS, Rulli NE, Taylor A, Sheng KC, Herrero LJ, et al. Arthritogenic alphaviral infection perturbs osteoblast function and triggers pathologic bone loss. Proc Natl Acad Sci U.S.A. (2014) 111(16):6040–5. doi: 10.1073/pnas.1318859111

PubMed Abstract | CrossRef Full Text | Google Scholar

75. Sugasti-Salazar M, Llamas-Gonzalez YY, Campos D, Gonzalez-Santamaria J. Inhibition of p38 mitogen-activated protein kinase impairs mayaro virus replication in human dermal fibroblasts and HeLa cells. Viruses (2021) 13(6):1156. doi: 10.3390/v13061156

PubMed Abstract | CrossRef Full Text | Google Scholar

76. Mazzon M, Castro C, Thaa B, Liu L, Mutso M, Liu X, et al. Alphavirus-induced hyperactivation of PI3K/AKT directs pro-viral metabolic changes. PloS Pathog (2018) 14(1):e1006835. doi: 10.1371/journal.ppat.1006835

PubMed Abstract | CrossRef Full Text | Google Scholar

77. Suchowiecki K, Reid SP, Simon GL, Firestein GS, Chang A. Persistent joint pain following arthropod virus infections. Curr Rheumatol Rep (2021) 23(4):26. doi: 10.1007/s11926-021-00987-y

PubMed Abstract | CrossRef Full Text | Google Scholar

78. Teo TH, Chan YH, Lee WW, Lum FM, Amrun SN, Her Z, et al. Fingolimod treatment abrogates chikungunya virus-induced arthralgia. Sci Transl Med (2017) 9(375):eaal1333. doi: 10.1126/scitranslmed.aal1333

PubMed Abstract | CrossRef Full Text | Google Scholar

79. Padmakumar B, Jayan JB, Menon RMR, Krishnankutty B, Payippallil R, Nisha RS. Comparative evaluation of four therapeutic regimes in chikungunya arthritis: a prospective randomized parallel-group study. Indian J Rheumatol (2009) 4(3):94–101. doi: 10.1016/S0973-3698(10)60189-6

CrossRef Full Text | Google Scholar

80. Mylonas AD, Harley D, Purdie DM, Pandeya N, Vecchio PC, Farmer JF, et al. Corticosteroid therapy in an alphaviral arthritis. J Clin Rheumatol (2004) 10(6):326–30. doi: 10.1097/01.rhu.0000147052.11190.36

PubMed Abstract | CrossRef Full Text | Google Scholar

81. Prow NA, Tang B, Gardner J, Le TT, Taylor A, Poo YS, et al. Lower temperatures reduce type I interferon activity and promote alphaviral arthritis. PloS Pathog (2017) 13(12):e1006788. doi: 10.1371/journal.ppat.1006788

PubMed Abstract | CrossRef Full Text | Google Scholar

82. Jones AC, Anderson D, Galbraith S, Fantino E, Gutierrez Cardenas D, Read JF, et al. Personalized transcriptomics reveals heterogeneous immunophenotypes in children with viral bronchiolitis. Am J Respir Crit Care Med (2019) 199(12):1537–49. doi: 10.1164/rccm.201804-0715OC

PubMed Abstract | CrossRef Full Text | Google Scholar

83. Giroux NS, Ding S, McClain MT, Burke TW, Petzold E, Chung HA, et al. Differential chromatin accessibility in peripheral blood mononuclear cells underlies COVID-19 disease severity prior to seroconversion. Sci Rep (2022) 12(1):11714. doi: 10.1038/s41598-022-15668-8

PubMed Abstract | CrossRef Full Text | Google Scholar

84. Sumitomo S, Nagafuchi Y, Tsuchida Y, Tsuchiya H, Ota M, Ishigaki K, et al. Transcriptome analysis of peripheral blood from patients with rheumatoid arthritis: a systematic review. Inflammation Regener (2018) 38:21. doi: 10.1186/s41232-018-0078-5

CrossRef Full Text | Google Scholar

86. Fontaine DA, Davis DB. Attention to background strain is essential for metabolic research: C57BL/6 and the international knockout mouse consortium. Diabetes (2016) 65(1):25–33. doi: 10.2337/db15-0982

PubMed Abstract | CrossRef Full Text | Google Scholar

87. Kraev A. Parallel universes of black six biology. Biol Direct (2014) 9:18. doi: 10.1186/1745-6150-9-18

PubMed Abstract | CrossRef Full Text | Google Scholar

88. Chopra A, Saluja M, Venugopalan A. Effectiveness of chloroquine and inflammatory cytokine response in patients with early persistent musculoskeletal pain and arthritis following chikungunya virus infection. Arthritis Rheumatol (2014) 66(2):319–26. doi: 10.1002/art.38221

PubMed Abstract | CrossRef Full Text | Google Scholar

89. Hoarau JJ, Jaffar Bandjee MC, Krejbich Trotot P, Das T, Li-Pat-Yuen G, Dassa B, et al. Persistent chronic inflammation and infection by chikungunya arthritogenic alphavirus in spite of a robust host immune response. J Immunol (2010) 184(10):5914–27. doi: 10.4049/jimmunol.0900255

PubMed Abstract | CrossRef Full Text | Google Scholar

90. Mylonas AD, Brown AM, Carthew TL, McGrath B, Purdie DM, Pandeya N, et al. Natural history of Ross river virus-induced epidemic polyarthritis. Med J Aust (2002) 177(7):356–60. doi: 10.5694/j.1326-5377.2002.tb04837.x

PubMed Abstract | CrossRef Full Text | Google Scholar

91. Manzoor KN, Javed F, Ejaz M, Ali M, Mujaddadi N, Khan AA, et al. The global emergence of chikungunya infection: An integrated view. Rev Med Virol (2022) 32(3):e2287. doi: 10.1002/rmv.2287

PubMed Abstract | CrossRef Full Text | Google Scholar

92. Robin S, Ramful D, Le Seach F, Jaffar-Bandjee M-C, Rigou G, Alessandri J-L. Neurologic manifestations of pediatric chikungunya infection. J Child Neurol (2008) 23(9):1028–35. doi: 10.1177/0883073808314151

PubMed Abstract | CrossRef Full Text | Google Scholar

93. Rafael KC, Preciado-Llanes L, Azar SR, Kim YC, Brandon O, Lopez-Camacho C, et al. Adenoviral-vectored mayaro and chikungunya virus vaccine candidates afford partial cross-protection from lethal challenge in A129 mouse model. Front Immunol (2020) 11:591885. doi: 10.3389/fimmu.2020.591885

PubMed Abstract | CrossRef Full Text | Google Scholar

94. Couderc T, Chretien F, Schilte C, Disson O, Brigitte M, Guivel-Benhassine F, et al. A mouse model for chikungunya: young age and inefficient type-I interferon signaling are risk factors for severe disease. PloS Pathog (2008) 4(2):e29. doi: 10.1371/journal.ppat.0040029

PubMed Abstract | CrossRef Full Text | Google Scholar

95. Baxter VK, Heise MT. Immunopathogenesis of alphaviruses. Adv Virus Res (2020) 107:315–82. doi: 10.1016/bs.aivir.2020.06.002

PubMed Abstract | CrossRef Full Text | Google Scholar

96. Fox JM, Diamond MS. Immune-mediated protection and pathogenesis of chikungunya virus. J Immunol (2016) 197(11):4210–8. doi: 10.4049/jimmunol.1601426

PubMed Abstract | CrossRef Full Text | Google Scholar

97. Farmer JF, Suhrbier A. Interpreting paired serology for Ross river virus and barmah forest virus diseases. Aust J Gen Pract (2019) 48(9):645–9. doi: 10.31128/AJGP-02-19-4845

PubMed Abstract | CrossRef Full Text | Google Scholar

98. Joubert PE, Werneke S, de la Calle C, Guivel-Benhassine F, Giodini A, Peduto L, et al. Chikungunya-induced cell death is limited by ER and oxidative stress-induced autophagy. Autophagy (2012) 8(8):1261–3. doi: 10.4161/auto.20751

PubMed Abstract | CrossRef Full Text | Google Scholar

99. Krejbich-Trotot P, Denizot M, Hoarau JJ, Jaffar-Bandjee MC, Das T, Gasque P. Chikungunya virus mobilizes the apoptotic machinery to invade host cell defenses. FASEB J (2011) 25(1):314–25. doi: 10.1096/fj.10-164178

PubMed Abstract | CrossRef Full Text | Google Scholar

100. Judith D, Mostowy S, Bourai M, Gangneux N, Lelek M, Lucas-Hourani M, et al. Species-specific impact of the autophagy machinery on chikungunya virus infection. EMBO Rep (2013) 14(6):534–44. doi: 10.1038/embor.2013.51

PubMed Abstract | CrossRef Full Text | Google Scholar

101. Kramnik I, Beamer G. Mouse models of human TB pathology: roles in the analysis of necrosis and the development of host-directed therapies. Semin Immunopathol (2016) 38(2):221–37. doi: 10.1007/s00281-015-0538-9

PubMed Abstract | CrossRef Full Text | Google Scholar

102. Perlman RL. Mouse models of human disease: An evolutionary perspective. Evol Med Public Health (2016) 2016(1):170–6. doi: 10.1093/emph/eow014

PubMed Abstract | CrossRef Full Text | Google Scholar

103. Andrews S. FastQC: a quality control tool for high throughput sequence data version 0.11.8 (2010). Available at: https://wwwbioinformaticsbabrahamacuk/projects/fastqc/.

Google Scholar

104. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for illumina sequence data. Bioinformatics (2014) 30(15):2114–20. doi: 10.1093/bioinformatics/btu170

PubMed Abstract | CrossRef Full Text | Google Scholar

105. Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics (2013) 29(1):15–21. doi: 10.1093/bioinformatics/bts635

PubMed Abstract | CrossRef Full Text | Google Scholar

106. Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The sequence Alignment/Map format and SAMtools. Bioinformatics (2009) 25(16):2078–9. doi: 10.1093/bioinformatics/btp352

PubMed Abstract | CrossRef Full Text | Google Scholar

107. Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-seq data with or without a reference genome. BMC Bioinf (2011) 12(1):323. doi: 10.1186/1471-2105-12-323

CrossRef Full Text | Google Scholar

108. Robinson MD, McCarthy DJ, Smyth GK. edgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics (2010) 26(1):139–40. doi: 10.1093/bioinformatics/btp616

PubMed Abstract | CrossRef Full Text | Google Scholar

109. The R Core Team. R: A language and environment for statistical computing version 4.1.0. r foundation for statistical computing (2021). Available at: https://www.R-project.org/.

Google Scholar

110. Durinck S, Spellman PT, Birney E, Huber W. Mapping identifiers for the integration of genomic datasets with the R/Bioconductor package biomaRt. Nat Protoc (2009) 4(8):1184–91. doi: 10.1038/nprot.2009.97

PubMed Abstract | CrossRef Full Text | Google Scholar

111. Larsson J. Eulerr: Area-proportional Euler and Venn diagrams with ellipses. In: R package version 6.1.1 (2021). Available at: https://CRAN.R-project.org/package=eulerr.

Google Scholar

112. Wickham H. ggplot2: elegant graphics for data analysis (use r!). New York: Springer (2009).

Google Scholar

113. Kassambara A. Ggpubr: 'ggplot2' based publication ready plots. In: R package version 0.4.0 (2020). Available at: https://CRAN.R-project.org/package=ggpubr.

Google Scholar

114. Krassowski M. Complex-upset. Create complex UpSet plots using ggplot2 components. In: R package version 1.1.3 (2020). doi: 10.5281/zenodo.3700590

CrossRef Full Text | Google Scholar

115. Blighe K, Lun A. PCAtools: Everything principal components analysis. In: R package version 2.4.0 (2021). Available at: https://github.com/kevinblighe/PCAtools.

Google Scholar

116. Neuwirth E. RColorBrewer: ColorBrewer palettes. In: R package version 1.1-2 (2014). Available at: https://CRAN.R-project.org/package=RColorBrewer.

Google Scholar

117. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U.S.A. (2005) 102(43):15545–50. doi: 10.1073/pnas.0506580102

PubMed Abstract | CrossRef Full Text | Google Scholar

118. Kolde R. Pheatmap: Pretty heatmaps. In: R package version 1.0.12 (2019). Available at: https://CRAN.R-project.org/package=pheatmap.

Google Scholar

Keywords: Chikungunya, C57BL/6J, mouse model, RNA-seq, bioinformatics, arthritis

Citation: Bishop CR, Caten FT, Nakaya HI and Suhrbier A (2022) Chikungunya patient transcriptional signatures faithfully recapitulated in a C57BL/6J mouse model. Front. Immunol. 13:1092370. doi: 10.3389/fimmu.2022.1092370

Received: 08 November 2022; Accepted: 25 November 2022;
Published: 12 December 2022.

Edited by:

Joel Henrique Ellwanger, Federal University of Rio Grande do Sul, Brazil

Reviewed by:

Penghua Wang, University of Connecticut Health Center, United States
Scott B. Biering, University of California, Berkeley, United States
Pierre Roques, Institut Pasteur de Guinée, Guinea

Copyright © 2022 Bishop, Caten, Nakaya and Suhrbier. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Helder I. Nakaya, aG5ha2F5YUBnbWFpbC5jb20=; Andreas Suhrbier, QW5kcmVhcy5TdWhyYmllckBxaW1yYmVyZ2hvZmVyLmVkdS5hdQ==

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.