Skip to main content

Cytomegalovirus infection disrupts the influence of short-chain fatty acid producers on Treg/Th17 balance



Both the gut microbiota and chronic viral infections have profound effects on host immunity, but interactions between these influences have been only superficially explored. Cytomegalovirus (CMV), for example, infects approximately 80% of people globally and drives significant changes in immune cells. Similarly, certain gut-resident bacteria affect T-cell development in mice and nonhuman primates. It is unknown if changes imposed by CMV on the intestinal microbiome contribute to immunologic effects of the infection.


We show that rhesus cytomegalovirus (RhCMV) infection is associated with specific differences in gut microbiota composition, including decreased abundance of Firmicutes, and that the extent of microbial change was associated with immunologic changes including the proliferation, differentiation, and cytokine production of CD8+ T cells. Furthermore, RhCMV infection disrupted the relationship between short-chain fatty acid producers and Treg/Th17 balance observed in seronegative animals, showing that some immunologic effects of CMV are due to disruption of previously existing host-microbe relationships.


Gut microbes have an important influence on health and disease. Diet is known to shape the microbiota, but the influence of concomitant chronic viral infections is unclear. We found that CMV influences gut microbiota composition to an extent that is correlated with immunologic changes in the host. Additionally, pre-existing correlations between immunophenotypes and gut microbes can be subverted by CMV infection. Immunologic effects of CMV infection on the host may therefore be mediated by two different mechanisms involving gut microbiota.

Video Abstract


Cytomegalovirus (CMV) is a β-herpesvirus that infects an estimated 80% of the global population [1] and drives significant changes in immune-cell phenotypes and functions [2]. Despite an arsenal of host immune responses against CMV, the infection is never cleared; instead, the virus establishes latency and can reactivate and cause disease when the host is immunocompromised, such as in cases of congenital infection, organ transplants, or acquired immunodeficiency syndrome (AIDS) [3]. Throughout the course of asymptomatic and predominantly latent infection, periodic viral reactivation can occur and stimulate the immune system, resulting in up to 10% of both CD4+ and CD8+ memory T-cell populations having specificity for CMV [4]. A large effect of persistent viral infection on the adaptive immune compartments might be expected, but surprisingly, CMV has equally substantial, durable effects on innate immune cells. Despite relative infrequency of CMV infection among antigen-presenting cells (APC), for example, the phenotype of these cells is transformed after infection [5]. Similarly, natural killer (NK) cells are impacted, and the imposed changes appear lifelong [6]. The virus’s natural ability to promote large expansions of adaptive immune cells prompted investigations into using rhesus CMV (RhCMV) as a vaccine vector to protect against the simian immunodeficiency virus (SIV), a model for human immunodeficiency virus (HIV) [7, 8]. We showed that vaccination with RhCMV-vectored vaccines is associated with IL-15-dependent expansion of innate-memory cells with SIV killing function [9].

Other environmental factors such as the gut microbiome have been shown to interact with the host to modify cytokine production capacity [10], potentially modifying vaccine responses. We demonstrated that infant diet and associated gut microbiota differences significantly impact immune development over the first 3–5 years of life [11, 12]. Additionally, type-1 diabetes onset in genetically predisposed infants was associated with decreased alpha diversity of the gut microbiome [13]. Inoculation of germ-free mice with segmented filamentous bacteria was shown to be critical for the development of Th17 cells [14], which are important for maintenance of the intestinal barrier and robust mucosal immunity [15]. Given the importance of microbes in shaping the host immune system, it is not surprising that certain constituents of the gut microbiota, including members of the phyla Actinobacteria and Firmicutes, have been associated with better vaccine responses, while Proteobacteria and Bacteroidetes have been associated with poor responses [16].

Interactions between viral infections and the commensal microbiota further complicate their relationships with host immunity. The gut microbiota can promote viral infection in some contexts but suppress viral infection in others [17]. For example, short-chain fatty acids (SCFAs) produced by commensal bacteria have been shown to reactivate latent herpesviruses such as Epstein-Barr virus and Kaposi’s sarcoma-associated herpesvirus [18], but lactic acid and Lactobacillus cell-wall components have been shown to inhibit herpes simplex virus activity [19, 20]. The persistence of murine norovirus differs per gut microbiota composition [21]. Another mouse study found that previous infection with murine CMV altered responses to yellow fever vaccination [22]. Moreover, mice latently infected with murine CMV exhibited greater resistance to sublethal doses of Listeria monocytogenes and reduced Yersinia pestis replication and spread [23]. Subclinical infections with RhCMV in adult macaques have been shown to induce changes in the gut microbiota and result in reduced immune responses to influenza A vaccination [24]. Since both CMV and the gut microbiota induce signaling cascades that control immune responses, it is reasonable that the affected pathways may overlap and interact [25, 26].

Nonetheless, interactions between the gut microbiota and chronic viral infections have only been superficially explored. Rhesus macaques represent a particularly useful model for studying how chronic viral infections alter relationships between gut microbiota and host [27]. Rhesus immune cells and many features of the adaptive response are similar to those found in humans [28]. Studies of RhCMV have demonstrated effects on host immunity similar to those apparent after human CMV (HCMV) infection [29]. The pathogenesis of fetal infections with RhCMV and HCMV are similar, making RhCMV an ideal model for study of HCMV pathogenesis, vaccines, and effects on the immune system [30]. To investigate the impact of RhCMV on gut microbe-host relationships, we studied the gut microbiotas and immune systems of RhCMV-seropositive and -seronegative infant macaques. We found that RhCMV infection had a direct effect on abundance of certain bacterial taxa in the gut and altered relationships between gut microbial taxa and immune-cell subsets.


Study design

Rectal swabs and blood samples were collected from 5- to 11-month-old infant rhesus macaques that were seropositive (RhCMV+, n = 29) or seronegative (RhCMV−, n = 38) for RhCMV. Animals were co-housed in outdoor corrals; samples were collected from all available animals in the correct age range at the time of sampling, without selection. Rectal swabs were stored in RNAlater™ Storage Solution (Sigma-Aldrich) at −70 °C until DNA extraction. Peripheral blood mononuclear cells (PBMCs) were isolated by gradient density purification using Lymphocyte Separation Medium (MP Biomedicals, LLC), and then washed in medium containing fetal bovine serum and 10% dimethyl sulfoxide (DMSO) before cryopreservation in liquid nitrogen prior to analysis. Stool samples were also collected from a separate longitudinal study in which 24 adult female indoor-housed RhCMV-negative rhesus macaques were vaccinated with a RhCMV 68-1-based vaccine vector [31]. Baseline samples were collected 2 weeks prior to vaccination, and post-vaccination samples were collected 3 weeks post vaccination.

DNA extraction and 16S rDNA sequencing

DNA from rectal swabs and stools was extracted using the MoBio PowerSoil kit (Qiagen). Amplicon libraries were generated by amplifying the V3-V4 or V4 variable region of 16S rRNA genes using primers 319F and 806R or 515F and 806R, respectively [32]. Both forward and reverse primers contained a unique 8-nt barcode, a primer pad, a linker sequence, and the Illumina adaptor sequences. Each sample was barcoded with a unique forward and reverse barcode combination. PCR reactions contained 1 unit KAPA2G Robust HotStart Polymerase (Kapa Biosystems), 1.5 mM MgCl2, 10 μmol of each primer, 10 mM dNTPs, and 1 ul of DNA. PCR conditions were as follows: an initial incubation at 95 °C for 2 min; 30 cycles of 95 °C for 15 s, 50 °C for 20 s; 72 °C for 20 s, and a final extension at 72°C for 3 min. The final product was quantified on the Qubit instrument using the Qubit High Sensitivity DNA kit (Invitrogen), and individual amplicon libraries were pooled, cleaned by AMPure XP beads (Beckman Coulter), and sequenced using a 250 bp paired-end method on an Illumina MiSeq instrument in the Genome Center DNA Technologies Core, University of California, Davis. The quality of sequencing reads was checked using FastQC. Sequences were trimmed and annotated to the genus level using the DADA2 package [33] to Greengenes database version 13_8 [34] within the R 4.1.1 software [35] using RStudio [36]. The numbers of reads per sample per bacterial feature were stored in a matrix and used in downstream statistical analysis.

Immune-cell phenotyping by flow cytometry

Immune cells were stained and analyzed as described in previously published paper using samples from the same animal cohort [9]. Distribution of innate cells and T-cell subsets in peripheral blood samples and their activation status were determined by flow cytometry using freshly isolated PBMC samples. The following antibodies were used: anti-CD3-Alexa 700, anti-CD95-APC (clone DX2), anti-CD28-APC-H7 (clone CD28.2), anti-CD8–PE-Cy5.5 (clone 3B5), anti-CD4–BV650 (clone L200), anti-HLADR-ECD (clone Immu-357), anti-CD14-Qdot 605 (clone TüK4), anti-CD16-PacBlue (clone 3G8), anti-CD20-ECD (clone B9E9), anti-CD11c-AF700 (clone 3.9), anti-CD123-PerCP-Cy™5.5 (clone 7G3), anti-CD80-FITC (clone L307.4), anti-CD83-PE (clone HB15e), and anti-CD86-APC (clone FUN-1). A cell viability dye (Invitrogen Aqua LIVE/DEAD Fixable Dead Cell Stain) was included to discriminate live from dead cells. Cells were washed and permeabilized using a Fix/Perm kit (BioLegend) according to the manufacturer’s instructions, intracellularly stained with anti-Ki67–Alexa 488 (clone B56) and fixed in phosphate-buffered saline containing 1% paraformaldehyde. Data were acquired on Aria or Fortessa cytometers (BD Biosciences) and analyzed using FlowJo software version 10.3 (BD Life Sciences).

Intracellular cytokine immunostaining

Immune cells were stained and analyzed as described in previously published paper using samples from the same animal cohort [9]. To measure the level of cytokine production in response to mitogenic stimulation, PBMC (1 million cells) were incubated for 4 h at 37 °C with phorbol 12-myristate 13-acetate (50 ng/ml) and ionomycin (1 μg/ml) in complete RPMI 1640 medium and with GolgiPlug (5 μg/ml). Cells were washed and immunostained with anti-CD3-PacBlue (clone SP34-2), anti-CD8–PE-Cy5.5 (clone 3B5), anti-CD4–BV650 (clone L200), anti-CD95-APC (clone DX2), and anti-CD28-APC-H7 (clone CD28.2) and a stain reagent (Invitrogen Aqua LIVE/DEAD Fixable Dead Cell Stain) to exclude dead cells. Cells were washed and permeabilized using a Cytofix/Cytoperm kit (BD Biosciences) according to the manufacturer’s instructions and intracellularly immunostained with anti-IL-17–PE (clone eBio64DEC17), anti–IFN-γ-PE-Cy7 (clone B27), and anti-TNFα-Alexa Fluor 700 (clone Mab11). Cells were finally washed and fixed in phosphate-buffered saline containing 1% paraformaldehyde. Data were acquired on Aria or Fortessa cytometers (BD Biosciences) and analyzed using FlowJo software version 10.3 (BD Life Sciences).

Statistical analysis

R-4.1.1 in RStudio was used for all statistical analysis [35, 36]. R packages phyloseq [37], vegan [38], limma [39], glmnet [40], and pROC [41] were used to filter taxa with < 5% prevalence, calculate distance matrices, perform differential abundance analysis, perform elastic-net logistic regression, and assess robustness of elastic-net algorithm, respectively. When assessing relationships between bacterial taxa and immune parameters (Table 2), residuals of each model were checked for homoscedasticity (Breusch-Pagan test) and normal distribution (Shapiro-Wilk test). P-values were adjusted for multiple testing using false discovery rate estimation with the qvalue package [42]. All the data were analyzed, without elimination of any outliers.


Gut microbial communities of RhCMV+ and RhCMV− animals are similar at the phylum level but cluster separately at the genus level

To test the impact of RhCMV infection on gut microbiota, we performed 16S ribosomal RNA gene sequencing on rectal swabs from 67 macaques aged 5 to 11 months, housed outdoors, and screened for anti-RhCMV antibodies. 29 seropositive animals and 38 seronegative animals were sampled (Table 1). An average of 81,000 raw reads were acquired per sample; an average of 65% ± 8% were retained after processing with the DADA2 pipeline [33], resulting in a mean of 52,000 analyzed reads per sample. We used Greengenes database version 13_8 [34] for taxonomic assignment. Bacterial sequence counts were agglomerated to the genus level, and reads not classified to the genus level were assigned to the lowest taxonomic assignment. Taxa that were not present in at least 5% of the samples were filtered, resulting in 92 bacterial taxa for downstream analysis. Bacteroidetes and Firmicutes were the most abundant phyla in all samples, accounting for 49% and 47% of reads, respectively (Fig. 1A). Alpha diversities within RhCMV+ and RhCMV− animals were assessed using Shannon’s diversity index, evenness, and richness, and no significant differences were seen (Wilcoxon rank-sum test P > 0.05). Beta diversity at the genus level was assessed using complete-linking clustering and principal component analysis (PCA) with Aitchison distances. In this analysis, animals were shown to cluster by RhCMV serostatus (Fig. 1B, P = 0.01).

Table 1 Study groups
Fig. 1
figure 1

Gut microbial communities of RhCMV+ and RhCMV− animals are similar at the phylum level but cluster separately at the genus level. A Relative abundance of microbiota for each animal at the phylum level. Samples are ordered by complete-linkage clustering based on Aitchison distance. Animals are color coded as RhCMV− (red) and RhCMV+ (blue) on top of the bar graph. B PCA plot of all samples analyzed at the genus level, based on Aitchison distance

RhCMV-seropositive animals and those experimentally infected with RhCMV68-1 vaccine vector had significantly decreased abundances of bacteria from phylum Firmicutes

To detect specific differences in microbial communities associated with RhCMV infection, we analyzed genus-level data (86% of total reads) using linear modeling of abundances with variances moderated by an empirical Bayes procedure (limma-voom, ref. [39, 43]). We found that 11 of 62 detected bacterial genera, six of which were Firmicutes, had different abundance in animals infected with RhCMV (Fig. 2A, unadjusted P < 0.05). Four genera were determined to be significantly changed (adjusted P < 0.1) after adjustment for multiple comparisons. Butyrivibrio, Sarcina, and Blautia, all known short-chain fatty acid (SCFA) producers [44,45,46], were less abundant in RhCMV+ animals; Streptococcus was more abundant.

Fig. 2
figure 2

RhCMV-infected animals had significantly decreased abundance of bacteria from the order Clostridia. A Differentially abundant bacterial genera analyzed using the limma-voom pipeline with P < 0.05 (*adjusted P < 0.1). B Features selected by elastic-net regression to differentiate between RhCMV− and RhCMV+ animals. C PCA plot generated using log-transformed counts of genera selected by elastic net within the training set. D ROC curve to assess the robustness of elastic-net output. E Summary score for each animal generated by the elastic net. F PCA plot generated using log-transformed counts of genera selected by elastic-net within the testing set (RhCMV-vectored vaccine recipients)

We used elastic-net logistic regression to attempt better modeling of the RhCMV-associated microbiota despite multicollinearity in the dataset [47, 48]. A 10× cross-validation was performed on a grid of alpha values to determine the optimal alpha (mixture between ridge and lasso) and lambda. Eighteen bacterial genera were found to contribute to an optimal ensemble associated with RhCMV infection (Fig. 2B). Streptococcus, Aggregatibacter, Lachnospira, YRC22, Treponema, RFN20, Fibrobacter, and SMB53 had positive coefficients, signifying that they have positive correlations with RhCMV infection, while Campylobacter, Bacteroides, Parabacteroides, Dialister, Blautia, Ruminococcus, Lactobacillus, Succinivibrio, Faecalibacterium, and Butyrivibrio had negative coefficients. As expected, when abundances of these 18 genera were reduced to two principal components, plotting demonstrated substantial separation of the RhCMV-seronegative vs. -seropositive animals (Fig. 2C). To determine if RhCMV infection itself was the cause of such associations, we performed 16S rRNA sequencing of stool samples taken 2 weeks before vs. 3 weeks after administration of an RhCMV68-1 vaccine vector [7]. The elastic-net model that had been trained using samples from animals with natural infection (“training set,” above) was applied to the cohort of vaccinated animals (“testing set”) to determine if the 18 selected genera were good identifiers of previous RhCMV exposure via vaccination. The area under the receiver operating characteristic (ROC) curve was 0.87 with a confidence interval of 0.766–0.973, showing that the method was robust in differentiating RhCMV naïve vs. RhCMV-infected or -vaccinated animals (Fig. 2D). The application of elastic net resulted in a summary score calculated by the model to describe the degree of “RhCMVness” in the microbiota of individual animals (Fig. 2E), which showed clear separation between RhCMV− vs. RhCMV+ and pre- vs. post-vaccinated animals. The results were further confirmed by examination of a principal-component plot of the log2-transformed counts of important taxa as determined by elastic-net regression, in which pre- and post-vaccinated animals (“testing set,” Fig. 2F; PERMANOVA P = 0.001) clustered separately from each other.

Bacterial features were correlated with immune phenotypes

We described immunologic changes resulting from CMV infection in a previously published paper whose experiments were performed in parallel with the microbiome investigation [9]. Animals were seen to cluster according to RhCMV serostatus on a PCA plot summarizing immune-cell frequencies, demonstrating the large and consistent immunologic change imposed by this infection (Fig. 3A, PERMANOVA, P < 0.05; see ref. [2]). Specific differences in immune-cell populations between RhCMV+ and RhCMV− animals were consistent with those previously reported, including lower frequency of naïve and greater frequency of memory/effector CD4+ and CD8+ T cells (Fig. 3A; ref [23]). Proliferation of T cells, marked by Ki-67 expression, was also highly correlated with RhCMV infection (Fig. 3A).

Fig. 3
figure 3

Immune cell types correlated with microbiome constituents. A PCA plot of immune markers. Loadings of important differences associated with RhCMV infection are shown. B Immune cell subsets correlated with the CMV-microbe score shown in Fig. 2E

To test for a possible association between immunologic and microbial changes, we used the CMV-microbe score from Fig. 2E as a summary score for microbial change and then searched for associations between these scores and circulating immune cells. Frequencies of 48 out of 84 (57%) immune subsets examined were significantly associated with the CMV-microbe score (Spearman rank test, adjusted P < 0.05; Table S1 and examples in Fig. 3B), suggesting an important relationship between CMV serostatus, the gut microbiota, and immune markers.

Cytomegalovirus subverts relationships between gut bacteria and immunophenotypes

We hypothesized that the immunologic impact of RhCMV infection would be strong enough to swamp many preexisting effects of commensal gut bacteria on immunity, subverting existing gut microbiota-host relationships. We searched for features of the host immune system whose relationship to specific bacterial taxa was modified by RhCMV infection (indicated by a significant interaction term in a multivariate analysis accounting for age). To avoid spurious significance due to sparse read counts, we further filtered taxa not represented in at least 50% of samples, resulting in 2905 total comparisons between 35 bacterial genera and 83 immune markers. Sixty-seven comparisons were statistically significant for both the overall F-test (P < 0.05) and the interaction term (P < 0.05) while also passing basic model checks (Table 2; see “Methods”). Most of the significantly altered relationships detected (44 of 63) reflected significant correlations in the RhCMV-seronegative animals that were not observed in the seropositive cohort, suggesting destruction of a preexisting relationship by RhCMV, confirming our hypothesis (Table 2, marked “RhCMV− only” at right, and examples in Fig. 4A, top section). Many such cases involved CD4+ and CD8+ T-cell populations known to be highly impacted by RhCMV. Specifically, among RhCMV-seronegative animals, cytokine-producing T-cell subsets were negatively correlated with genera from the Lachnospiraceae family (Roseburia, Oribacterium, Coprococcus, Lachnospira, and Dorea), often short-chain fatty acid producers [44, 45], but these relationships were diminished in seropositive animals (Table 2 and examples in Fig. 4A, top section). Additionally, in seronegative macaques, SCFA producers were associated with Th17/Treg balance through positive correlations with Tregs (marked “CD25+, CD127lo, CD4”) and negative correlations with Th17 cells (“IL17+, CD4”; Fig. 4A, top section, and Fig. 4B). In rarer cases, relationships not seen in seronegative macaques were present in seropositive macaques (16 of 67 interactions; Fig. 4A, middle section). Genera from the Veillonellaceae family (Veillonella, Megasphaera, and Anaerovibrio) were seen to have more impact on the immune system in seropositive macaques (Table 2 and examples in Fig. 4A, middle section), suggesting that in CMV infection, there is a shift within class Clostridia from dominant effects of Lachnospiraceae on immunity to greater impact of Veillonellaceae. Seven of 67 interactions were significantly correlated in both RhCMV− and RhCMV+ animals, but with opposite polarity (Table 2; Fig. 4A, bottom section; and Fig. 4 B–C). Oribacterium, Roseburia, and Faecalibacterium are all known short-chain fatty acid producers. As expected, these genera have negative correlations with cytokine-producing CD8+ T cells in RhCMV-seronegative animals; however, the opposite is true for seropositive animals (Fig. 4 B–C).

Table 2 Results of multiple linear regression resulting in significant interaction due to RhCMV infection
Fig. 4
figure 4

RhCMV subverts relationships between immune-cell subsets and gut bacteria. A Significant associations between immune cells and gut bacterial abundance that were changed by RhCMV infection, as indicated by a significant interaction term in regression. The significance of interaction term was indicated by different sizes (largest have P < 0.05, medium have 0.05 ≤ P < 0.1, smallest have P ≥ 0.1). Top panels show significant bacterial-immune correlations within RhCMV− animals only; middle panels show significant bacterial-immune correlations within RhCMV+ animals only; bottom panels show significant bacterial-immune correlations in both RhCMV+ and RhCMV− animals. Red signifies negative correlations, while blue signifies positive correlations. aGram-positive and bgram-negative bacterial genera known to produce SCFA. B Significant correlations between Treg (CD25+CD127lo CD4+ T cells) or Th17 (IL17+ CD4+ T cells) and known SCFA producers Bulleidia, Dialister, Oribacterium, and Faecalibacterium in RhCMV− animals but not RhCMV+ animals, also shown in the top panel of Fig. 4A. C Significant immune correlations detected in both RhCMV− and RhCMV+ animals with Oribacterium, Roseburia, and Faecalibacterium, also shown in the bottom panel of Fig. 4A


We found that CMV infection is associated with alterations in both gut microbes and relationships between microbes and host immunity. A substantial literature demonstrates the effects of both human and rhesus cytomegaloviruses on host immune function [4, 9, 24, 29, 49, 50]. Brodin and colleagues showed that HCMV infection impacted more than 50% of the immune parameters examined [2]. Our data suggest the possibility that altered host sensitivity to microbial constituents forms part of the mechanism of these CMV effects—that is, CMV may act partly via changed sensitivity of the host to its own microbiota. Thus, understanding the effects of CMV on host immunity requires full understanding of interactions between chronic viral infections, microbes, and the immune system.

While previous studies did not find drastically changed host microbiota during RhCMV infection as assessed by beta diversity measures [24, 51, 52], we found significant differences between seronegative and seropositive microbiotas using methods for compositional data analysis [53], albeit with some overlap between the two groups (Fig. 1B). RhCMV-seropositive animals demonstrated increased abundance of Streptococcus and decreased abundances of Blautia, Butyrivibrio, and Sarcina—consistent with a previous report of decreased Firmicutes abundance in CMV-infected human infants [52]. We used different methods to identify differentially abundant individual taxa (limma-voom) vs. compositionally altered microbial communities (elastic net). Elastic-net regression identified an additional 15 genera that were associated with RhCMV seropositivity, while agreeing with many individual results from the limma-voom analysis, e.g., Streptococcus and Butyrivibrio. A separate longitudinal cohort of macaques administered RhCMV-vectored vaccines was used to test the machine-learning results; these experimentally infected macaques demonstrated similar changes in their microbiotas, indicating a causal relationship between RhCMV infection and microbiota changes.

All taxa that declined in abundance during RhCMV infection are SCFA-producing genera [44,45,46]. Changes in abundances of SCFA producers in the gut are commonly associated with altered diets, with higher fiber intake resulting in greater abundance of SCFA producers and increased host health [54, 55]. All animals in our study were fed the same diet, however, likely increasing sensitivity of the experiment to RhCMV-imposed changes. Previous research has demonstrated reduced abundance of SFCA producers in people suffering from inflammatory bowel disease, especially Crohn’s disease [56, 57], suggesting that reduced abundance of these taxa is a correlate of systemic inflammation. In a parallel line of experimentation, we demonstrated that RhCMV infection induces outgrowth of innate-memory CD8+ cells through an IL-15-dependent pathway [9]. While IL-15 may be important for protection against SIV, the cytokine has also been shown to promote intestinal dysbiosis, manifesting as reduction in SCFA producers [58], similar to the finding in this study (Fig. 2 A–B). We hypothesize that the changes in bacterial abundance reported above are due to altered host-gene expression, including IL-15 expression, in gut tissue that is caused by the immune response to RhCMV.

Interruptions of microbial influence on the host immune system due to RhCMV infection were observed mostly in T-cell subsets. Previously, the microbiome has been implicated in the development of mucosal Th17 cells and Tregs [59, 60]. Th17 cells and Tregs, present ubiquitously at the mucosal surface, are both induced by TGF-β signaling during activation. In the presence of IL-6, naïve CD4+ T cells commit to the Th17 lineage, while in the absence of IL-6 and other pro-inflammatory molecules, to the Treg lineage [61]. A study in mice found that supplementation with Faecalibacterium significantly decreased both IL-17A expression and levels of IL-17 in the plasma [62]. We similarly observed inverse correlations between abundances of Faecalibacterium, Dialister, and Oribacterium and the frequency of circulating Th17 cells among RhCMV-seronegative (but not -seropositive) animals. Another SCFA-producing genus, Bulleidia [63], positively correlated with Tregs in RhCMV-seronegative animals but not RhCMV-seropositive animals (Fig. 4 A–B). RhCMV infection therefore disrupts the relationship between SCFA producers and Th17/Treg balance, supplanting a homeostatic mechanism that controls this balance in seronegative macaques.

While in this example RhCMV infection subverts a preexisting gut microbe-immune system relationship, in other instances, RhCMV appears to sensitize its host to microbial influence: Anaerovibrio, Megasphaera, and Veillonella are important correlates of certain immunophenotypes only in RhCMV-seropositive macaques. While previously classified in the Clostridiales order based on 16S rRNA gene sequences and metabolic characteristics (all are SCFA producers) [46], these genera are gram negative in contrast to other members of Clostridiales and have been reclassified to a novel bacterial order, Negativicutes [64]. Among RhCMV-seropositive macaques, Anaerovibrio positively correlates with T-cell proliferation, while Megasphaera and Veillonella positively correlate with T-cell effector functions. Lipopolysaccharide (LPS), a key component in the cell wall of gram-negative bacteria, stimulates T cells via engagement of toll-like receptor 4 (TLR4) within innate immune cells [65, 66]. In addition, studies have shown that LPS promotes immediate-early gene expression of HCMV [67] and reactivation of latent CMV in mice [68]. Thus, CMV and the altered microbiota may synergistically promote T-cell activation, which in turn promotes further CMV replication.

Surprisingly, RhCMV infection inverts the relationship between abundance of SCFA producers and effector CD8+ T cells that is seen in RhCMV-naïve macaques. Oribacterium, Roseburia, and Faecalibacterium negatively correlate with the frequency of effector CD8+ T cells in seronegative macaques but positively correlate in RhCMV-seropositive animals (Fig. 4A). SCFA are thought to reduce inflammation via blockade of NF-κB activation [69,70,71,72] and inhibition of histone deacetylase [73,74,75,76,77]. However, SCFA can also cause inflammation when other TLR agonists are present [78] and cause inflammation when interacting with different G protein-coupled receptors [79, 80], suggesting that the net effect of SCFA partly depends on microenvironmental factors. In fact, SCFA has been shown to enhance the susceptibility to and induce the replication of CMV in human cell lines [81,82,83,84]. Therefore, CMV infection may compromise the anti-inflammatory signaling normally generated by the host’s endogenous microbiota and harness the SCFA produced by these bacteria to enhance its own survival.

In summary, we found that RhCMV infection was associated with a profound change in the relationship of an infected host to its microbiota. The methods we employed likely limited sensitivity of our study to other examples of altered host-microbe relationships, as 16S rRNA sequencing surveys bacterial taxonomy but not bacterial functions. A future study using metagenomic techniques would better resolve functions between bacteria and would possibly reveal more bacteria-immune system relationships that were impacted by RhCMV. Bacteria within a genus, or even within a species, maybe have different functions and consequently different effects on the host [85]. Additionally, gut-specific immune and transcriptomic data would provide mechanistic insights into the local interactions occurring between microbiome and host during RhCMV infection. Previous studies have shown that depleting Tregs in the salivary gland but not in the spleen resulted in CMV reactivation at one site but not the other [86], so profiling the gut immune system would likely reveal other interactions with CMV that are not seen in circulating cells. Nevertheless, our work provides a framework for discovering interactions between chronic viral infections and gut commensal microorganisms in vivo. Our future studies will use this framework to understand the importance of RhCMV for adaptive immune responses to infection or vaccination.


RhCMV infection is associated with a profound change in the relationship of an infected host to its microbiota. SCFA-producing genera are found in lower abundance in RhCMV-infected macaques. Furthermore, relationships between gut bacteria and host immune functions are disrupted. We observed inverse correlations between abundances of SFCA producers and the frequency of circulating Th17 cells among RhCMV-seronegative—but not -seropositive—animals. Another SCFA-producing genus, Bulleidia [63], positively correlated with Tregs in RhCMV-seronegative but not RhCMV-seropositive animals. In addition, RhCMV infection inverts the relationship between abundance of SCFA producers and effector CD8+ T cells that is seen in RhCMV-uninfected macaques. Oribacterium, Roseburia, and Faecalibacterium negatively correlated with the frequency of effector CD8+ T cells in seronegative macaques but positively correlate in RhCMV-seropositive animals. RhCMV infection thus compromises the anti-inflammatory signaling normally generated by the host’s endogenous microbiota and harnesses the SCFA produced by these bacteria to enhance its own survival. Thus, part of the mechanism of vast CMV effects on host immunity is alteration of the host-microbiome relationship.

Availability of data and materials

The datasets generated and analyzed during the current study are available in the SRA under BioProject ID PRJNA789338 and GitHub (


  1. Zuhair M, Smit GSA, Wallis G, Jabbar F, Smith C, Devleesschauwer B, et al. Estimation of the worldwide seroprevalence of cytomegalovirus: A systematic review and meta-analysis. Rev Med Virol. 2019;29:e2034.

  2. Brodin P, Jojic V, Gao T, Bhattacharya S, Angel CJL, Furman D, et al. Variation in the human immune system is largely driven by non-heritable influences. Cell. 2015;160:37–47.

  3. Dupont L, Reeves MB. Cytomegalovirus latency and reactivation: recent insights into an age old problem. Rev Med Virol. 2016;26:75–89.

    Article  PubMed  Google Scholar 

  4. Sylwester AW, Mitchell BL, Edgar JB, Taormina C, Pelte C, Ruchti F, et al. Broadly targeted human cytomegalovirus-specific CD4+ and CD8+ T cells dominate the memory compartments of exposed subjects. J Exp Med. 2005;202:673–85.

  5. Sinclair J, Reeves M. The intimate relationship between human cytomegalovirus and the dendritic cell lineage. Front Microbiol. 2014;5:1–14.

    Article  CAS  Google Scholar 

  6. Lee J, Zhang T, Hwang I, Kim A, Nitschke L, Kim M, et al. Epigenetic modification and antibody-dependent expansion of memory-like NK cells in human cytomegalovirus-infected individuals. Immunity. 2015;42:431–42.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Hansen SG, Ford JC, Lewis MS, Ventura AB, Hughes CM, Coyne-Johnson L, et al. Profound early control of highly pathogenic SIV by an effector memory T-cell vaccine. Nature. 2011;473:523–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Früh K, Picker L. CD8+ T cell programming by cytomegalovirus vectors: applications in prophylactic and therapeutic vaccination. Curr Opin Immunol. 2017;47:52–6.

    Article  Google Scholar 

  9. Méndez-Lagares G, Chin N, Chang WLW, Lee J, Rosás-Umbert M, Kieu HT, et al. Cytomegalovirus mediates expansion of IL-15–responsive innate-memory cells with SIV killing function. J Clin Invest. 2021;131.

  10. Schirmer M, Smeekens SP, Vlamakis H, Jaeger M, Oosting M, Franzosa EA, et al. Linking the human gut microbiome to inflammatory cytokine production capacity. Cell. 2016;167:1125–1136.e8.

    Article  CAS  Google Scholar 

  11. Ardeshir A, Narayan NR, Méndez-Lagares G, Lu D, Rauch M, Huang Y, et al. Breast-fed and bottle-fed infant rhesus macaques develop distinct gut microbiotas and immune systems. Sci Transl Med. 2014;6.

  12. Narayan NR, Méndez-Lagares G, Ardeshir A, Lu D, Van Rompay KKA, Hartigan-O’Connor DJ. Persistent effects of early infant diet and associated microbiota on the juvenile immune system. Gut Microbes. 2015;6:284–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Kostic AD, Gevers D, Siljander H, Vatanen T, Hyötyläinen T, Hämäläinen A-MM, et al. The dynamics of the human infant gut microbiome in development and in progression toward type 1 diabetes. Cell Host Microbe. 2015;17:260–73.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Ivanov II, Atarashi K, Manel N, Brodie EL, Shima T, Karaoz U, et al. Induction of intestinal Th17 cells by segmented filamentous bacteria. Cell. 2009;139:485–98.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Hartigan-O’Connor DJ, Abel K, Van Rompay KKA, Kanwar B, McCune JM. SIV Replication in the infected rhesus macaque is limited by the size of the preexisting T H 17 cell compartment. Sci Transl Med. 2012;4:136ra69.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Zimmermann P, Curtis N. The influence of the intestinal microbiome on vaccine responses. Vaccine. 2018;36:4433–9.

    Article  CAS  PubMed  Google Scholar 

  17. Li N, Ma W-T, Pang M, Fan Q-L, Hua J-L. The commensal microbiota and viral infection: a comprehensive review. Front Immunol. 2019:1–16.

  18. Gorres KL, Daigle D, Mohanram S, Miller G. Activation and repression of Epstein-Barr virus and Kaposi’s sarcoma-associated herpesvirus lytic cycles by short- and medium-chain fatty acids. J Virol. 2014;88:8028–44.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Conti C, Malacrino C, Mastromarino P. Inhibition of herpes simplex virus type 2 by vaginal lactobacilli. J Physiol Pharmacol. 2009;60(SUPPL.6):19–26.

    PubMed  Google Scholar 

  20. Mastromarino P, Cacciotti F, Masci A, Mosca L. Antiviral activity of Lactobacillus brevis towards herpes simplex virus type 2: role of cell wall associated components. Anaerobe. 2011;17:334–6.

    Article  PubMed  Google Scholar 

  21. Baldridge MT, Nice TJ, McCune BT, Yokoyama CC, Kambal A, Wheadon M, et al. Commensal microbes and interferon-λ determine persistence of enteric murine norovirus infection. Science. 2015;347:266–9.

    Article  CAS  Google Scholar 

  22. Reese TA, Bi K, Kambal A, Filali-Mouhim A, Beura LK, Bürger MC, et al. Sequential infection with common pathogens promotes human-like immune gene expression and altered vaccine response. Cell Host Microbe. 2016;19:713–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Barton ES, White DW, Cathelyn JS, Brett-McClellan KA, Engle M, Diamond MS, et al. Herpesvirus latency confers symbiotic protection from bacterial infection. Nature. 2007;447:326–9.

    Article  CAS  PubMed  Google Scholar 

  24. Santos Rocha C, Hirao LA, Weber MG, Méndez-Lagares G, Chang WLW, Jiang G, et al. Subclinical cytomegalovirus infection is associated with altered host immunity, gut microbiota, and vaccine responses. J Virol. 2018;92:e00167–18.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Noriega VM, Haye KK, Kraus TA, Kowalsky SR, Ge Y, Moran TM, et al. Human cytomegalovirus modulates monocyte-mediated innate immune responses during short-term experimental latency in vitro. J Virol. 2014;88:9391–405.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Larsson E, Tremaroli V, Lee YS, Koren O, Nookaew I, Fricker A, et al. Analysis of gut microbial regulation of host gene expression along the length of the gut and regulation of gut microbial ecology through MyD88. Gut. 2012;61:1124–31.

    Article  CAS  PubMed  Google Scholar 

  27. McKenna P, Hoffmann C, Minkah N, Aye PP, Lackner A, Liu Z, et al. The macaque gut microbiome in health, lentiviral infection, and chronic enterocolitis. PLoS Pathog. 2008;4:e20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Messaoudi I, Estep R, Robinson B, Wong SW. Nonhuman primate models of human immunology. Antioxid Redox Signal. 2011;14:261–73.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Oxford KL, dela Pena-Ponce MGA, Jensen K, Eberhardt MK, Spinner A, Van Rompay KKA, et al. The interplay between immune maturation, age, chronic viral infection and environment. Immun Ageing. 2015;12:3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Itell HL, Kaur A, Deere JD, Barry PA, Permar SR. Rhesus monkeys for a nonhuman primate model of cytomegalovirus infections. Curr Opin Virol. 2017;25 cCMV:126–33.

    Article  Google Scholar 

  31. Chang WLW, Deere JD, Kieu HT, Castillo LD, Machmach K, Shen X, et al. RhCMV serostatus and vaccine adjuvant impact immunogenicity of RhCMV/SIV vaccines. Sci Rep. 2020;10:14056.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Caporaso JG, Lauber CL, Walters WA, Berg-Lyons D, Lozupone CA, Turnbaugh PJ, et al. Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. Proc Natl Acad Sci U S A. 2011;108(SUPPL. 1):4516–22.

    Article  CAS  Google Scholar 

  33. Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJA, Holmes SP. DADA2: high-resolution sample inference from Illumina amplicon data. Nat Methods. 2016;13:581–3.

    Article  CAS  Google Scholar 

  34. DeSantis TZ, Hugenholtz P, Larsen N, Rojas M, Brodie EL, Keller K, et al. Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl Environ Microbiol. 2006;72:5069–72.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. R Core Team. R: A Language and Environment for Statistical Computing. 2021.

  36. RStudio Team. RStudio: Integrated Development for R. 2020.

  37. McMurdie PJ, Holmes S, Watson M. phyloseq: An R Package for Reproducible Interactive Analysis and Graphics of Microbiome Census Data. PLoS ONE. 2013;8(4):e61217.

  38. Oksanen J, Blanchet GF, Friendly M, Kindt R, Legendre P, McGlinn D, et al. vegan: Community Ecology Package. 2019.

  39. Law CW, Chen Y, Shi W, Smyth GK. Voom: precision weights unlock linear model analysis tools for RNA-seq read counts. Genome Biol. 2014;15:R29.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Friedman J, Hastie T, Tibshirani R. Regularization Paths for Generalized Linear Models via Coordinate Descent. J Stat Softw. 2010;33(1).

  41. Robin X, Turck N, Hainard A, Tiberti N, Lisacek F, Sanchez J, Müller M. pROC: an open-source package for R and S+ to analyze and compare ROC curves. BMC Bioinformatics. 2021;12(1):77.

  42. Storey JD, Bass AJ, Dabney A, Robinson D. qvalue: Q-value estimation for false discovery rate control. 2019.

  43. Zhang X, Nieuwdorp M, Groen AK, Zwinderman AH. Statistical evaluation of diet-microbe associations. BMC Microbiol. 2019;19:90.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Biddle A, Stewart L, Blanchard J, Leschine S. Untangling the genetic basis of fibrolytic specialization by Lachnospiraceae and Ruminococcaceae in diverse gut communities. Diversity. 2013;5:627–40.

    Article  Google Scholar 

  45. Markowiak-Kopeć P, Śliżewska K. The effect of probiotics on the production of short-chain fatty acids by human intestinal microbiome. Nutrients. 2020;12:1107.

    Article  CAS  PubMed Central  Google Scholar 

  46. Ríos-Covián D, Ruas-Madiedo P, Margolles A, Gueimonde M, de los Reyes-Gavilán CG, Salazar N. Intestinal short chain fatty acids and their link with diet and human health. Front Microbiol. 2016;7:FEB:1–9.

    Article  Google Scholar 

  47. Gonzales GB, De Saeger S. Elastic net regularized regression for time-series analysis of plasma metabolome stability under sub-optimal freezing condition. Sci Rep. 2018;8:3659.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Torang A, Gupta P, Klinke DJ. An elastic-net logistic regression approach to generate classifiers and gene signatures for types of immune cells and T helper cell subsets. BMC Bioinformatics. 2019;20:433.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Chidrawar S, Khan N, Wei W, McLarnon A, Smith N, Nayak L, et al. Cytomegalovirus-seropositivity has a profound influence on the magnitude of major lymphoid subsets within healthy individuals. Clin Exp Immunol. 2009;155:423–32.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Čičin-Šain L, Sylwester AW, Hagen SI, Siess DC, Currier N, Legasse AW, et al. Cytomegalovirus-specific T cell immunity is maintained in immunosenescent rhesus macaques. J Immunol. 2011;187:1722–32.

    Article  Google Scholar 

  51. Gianella S, Chaillon A, Mutlu EA, Engen PA, Voigt RM, Keshavarzian A, et al. Effect of cytomegalovirus and Epstein–Barr virus replication on intestinal mucosal gene expression and microbiome composition of HIV-infected and uninfected individuals. AIDS. 2017;31:2059–67.

    Article  PubMed  Google Scholar 

  52. Simmons KE. Early-life human cytomegalovirus infections in Canadians : implications for a healthy development: University of British Columbia; 2019.

    Book  Google Scholar 

  53. Gloor GB, Macklaim JM, Pawlowsky-Glahn V, Egozcue JJ. Microbiome datasets are compositional: and this is not optional. Front Microbiol. 2017:1–6.

  54. Baxter NT, Schmidt AW, Venkataraman A, Kim KS, Waldron C, Schmidt TM. Dynamics of human gut microbiota and short-chain fatty acids in response to dietary interventions with three fermentable fibers. MBio. 2019;10.

  55. Medawar E, Haange S-B, Rolle-Kampczyk U, Engelmann B, Dietrich A, Thieleking R, et al. Gut microbiota link dietary fiber intake and short-chain fatty acid metabolism with eating behavior. Transl Psychiatry. 2021;11:500.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Takahashi K, Nishida A, Fujimoto T, Fujii M, Shioya M, Imaeda H, et al. Reduced abundance of butyrate-producing bacteria species in the fecal microbial community in Crohn’s disease. Digestion. 2016;93:59–65.

    Article  CAS  PubMed  Google Scholar 

  57. Joossens M, Huys G, Cnockaert M, De Preter V, Verbeke K, Rutgeerts P, et al. Dysbiosis of the faecal microbiota in patients with Crohn’s disease and their unaffected relatives. Gut. 2011;60:631–7.

    Article  PubMed  Google Scholar 

  58. Meisel M, Mayassi T, Fehlner-Peach H, Koval JC, O’Brien SL, Hinterleitner R, et al. Interleukin-15 promotes intestinal dysbiosis with butyrate deficiency associated with increased susceptibility to colitis. ISME J. 2017;11:15–30.

    Article  CAS  PubMed  Google Scholar 

  59. Omenetti S, Pizarro TT. The Treg/Th17 axis: a dynamic balance regulated by the gut microbiome. Front Immunol. 2015;6:DEC:1–8.

    Article  CAS  Google Scholar 

  60. Pandiyan P, Bhaskaran N, Zou M, Schneider E, Jayaraman S, Huehn J. Microbiome dependent regulation of Tregs and Th17 cells in mucosa. Front Immunol. 2019;10:1–17.

    Article  CAS  Google Scholar 

  61. Kimura A, Kishimoto T. IL-6: regulator of Treg/Th17 balance. Eur J Immunol. 2010;40:1830–5.

    Article  CAS  PubMed  Google Scholar 

  62. Huang X-L. Faecalibacterium prausnitzii supernatant ameliorates dextran sulfate sodium induced colitis by regulating Th17 cell differentiation. World J Gastroenterol. 2016;22:5201.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Downes J, Olsvik B, Hiom SJ, Spratt DA, Cheeseman SL, Olsen I, et al. Bulleidia extructa gen. nov., sp. nov., isolated from the oral cavity. Int J Syst Evol Microbiol. 2000;50:979–83.

    Article  CAS  Google Scholar 

  64. Marchandin H, Teyssier C, Campos J, Jean-Pierre H, Roger F, Gay B, et al. Negativicoccus succinicivorans gen. nov., sp. nov., isolated from human clinical samples, emended description of the family Veillonellaceae and description of Negativicutes classis nov., Selenomonadales ord. nov. and Acidaminococcaceae fam. nov. in the ba. Int J Syst Evol Microbiol. 2010;60:1271–9.

    Article  CAS  Google Scholar 

  65. Tough DF, Sun S, Sprent J. T cell stimulation in vivo by lipopolysaccharide (LPS). J Exp Med. 1997;185:2089–94.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  66. Ulmer AJ, Flad H, Rietschel T, Mattern T. Induction of proliferation and cytokine production in human T lymphocytes by lipopolysaccharide (LPS). Toxicology. 2000;152:37–45.

    Article  CAS  PubMed  Google Scholar 

  67. Kew VG, Wills MR, Reeves MB. LPS promotes a monocyte phenotype permissive for human cytomegalovirus immediate-early gene expression upon infection but not reactivation from latency. Sci Rep. 2017;7:810.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  68. Cook CH, Trgovcich J, Zimmerman PD, Zhang Y, Sedmak DD. Lipopolysaccharide, tumor necrosis factor alpha, or interleukin-1β triggers reactivation of latent cytomegalovirus in immunocompetent mice. J Virol. 2006;80:9151–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  69. Lakhdari O, Tap J, Béguet-Crespel F, Le Roux K, de Wouters T, Cultrone A, et al. Identification of NF-κB modulation capabilities within human intestinal commensal bacteria. J Biomed Biotechnol. 2011;2011:1–9.

    Article  CAS  Google Scholar 

  70. Sokol H, Pigneur B, Watterlot L, Lakhdari O, Bermudez-Humaran LG, Gratadoux J-J, et al. Faecalibacterium prausnitzii is an anti-inflammatory commensal bacterium identified by gut microbiota analysis of Crohn disease patients. Proc Natl Acad Sci. 2008;105:16731–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  71. Pryde S. The microbiology of butyrate formation in the human colon. FEMS Microbiol Lett. 2002;217:133–9.

    Article  CAS  PubMed  Google Scholar 

  72. Koh A, De Vadder F, Kovatcheva-Datchary P, Bäckhed F. From dietary fiber to host physiology: short-chain fatty acids as key bacterial metabolites. Cell. 2016;165:1332–45.

    Article  CAS  PubMed  Google Scholar 

  73. Waldecker M, Kautenburger T, Daumann H, Busch C, Schrenk D. Inhibition of histone-deacetylase activity by short-chain fatty acids and some polyphenol metabolites formed in the colon. J Nutr Biochem. 2008;19:587–93.

    Article  CAS  PubMed  Google Scholar 

  74. Fellows R, Denizot J, Stellato C, Cuomo A, Jain P, Stoyanova E, et al. Microbiota derived short chain fatty acids promote histone crotonylation in the colon through histone deacetylases. Nat Commun. 2018;9:105.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  75. Davie JR. Inhibition of histone deacetylase activity by butyrate. J Nutr. 2003;133:2485S–93S.

    Article  CAS  PubMed  Google Scholar 

  76. Hamer HM, Jonkers D, Venema K, Vanhoutvin S, Troost FJ, Brummer R-J. Review article: the role of butyrate on colonic function. Aliment Pharmacol Ther. 2007;27:104–19.

    Article  CAS  PubMed  Google Scholar 

  77. Silva LG, Ferguson BS, Avila AS, Faciola AP. Sodium propionate and sodium butyrate effects on histone deacetylase (HDAC) activity, histone acetylation, and inflammatory gene expression in bovine mammary epithelial cells. J Anim Sci. 2018;96:5244–52.

    PubMed  PubMed Central  Google Scholar 

  78. Mirmonsef P, Zariffard MR, Gilbert D, Makinde H, Landay AL, Spear GT. Short-chain fatty acids induce pro-inflammatory cytokine production alone and in combination with toll-like receptor ligands. Am J Reprod Immunol. 2012;67:391–400.

    Article  CAS  PubMed  Google Scholar 

  79. Park J, Wang Q, Wu Q, Mao-Draayer Y, Kim CH. Bidirectional regulatory potentials of short-chain fatty acids and their G-protein-coupled receptors in autoimmune neuroinflammation. Sci Rep. 2019;9:8837.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  80. Li M, van Esch BCAM, Wagenaar GTM, Garssen J, Folkerts G, Henricks PAJ. Pro- and anti-inflammatory effects of short chain fatty acids on immune and endothelial cells. Eur J Pharmacol. 2018;831:52–9.

    Article  CAS  Google Scholar 

  81. Radsak K, Fuhrmann R, Franke RP, Schneider D, Kollert A, Brücher KH, et al. Induction by sodium butyrate of cytomegalovirus replication in human endothelial cells. Arch Virol. 1989;107:151–8.

    Article  CAS  PubMed  Google Scholar 

  82. Tanaka J, Sadanari H, Sato H, Fukuda S. Sodium butyrate-inducible replication of human cytomegalovirus in a human epithelial cell line. Virology. 1991;185:271–80.

    Article  CAS  PubMed  Google Scholar 

  83. Wu QH, Ascensao J, Almeida G, Forman SJ, Shanley JD. The effect of short-chain fatty acids on the susceptibility of human umbilical vein endothelial cells to human cytomegalovirus infection. J Virol Methods. 1994;47:37–50.

    Article  CAS  PubMed  Google Scholar 

  84. Vysochan A, Sengupta A, Weljie AM, Alwine JC, Yu Y. ACSS2-mediated acetyl-CoA synthesis from acetate is necessary for human cytomegalovirus infection. Proc Natl Acad Sci. 2017;114:E1528–35.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  85. Greenblum S, Carr R, Borenstein E. Extensive strain-level copy-number variation across human gut microbiome species. Cell. 2015;160:583–94.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  86. Almanan M, Raynor J, Sholl A, Wang M, Chougnet C, Cardin RD, et al. Tissue-specific control of latent CMV reactivation by regulatory T cells. PLoS Pathog. 2017;13:e1006507.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references


We thank Blythe Durbin-Johnson and Matt Settles for advice on statistical analysis. We thank the veterinary, Colony Research Clinical Laboratories, and pathology faculty and staff at CNPRC for expert technical assistance. The UC Davis Host Microbe Systems Biology Core and the CNPRC Flow Cytometry Core provided services supporting this work.


This study was funded by National Institutes of Health grants R01 AI150554, R01 AI143554, P01 AI131568, and R01 AI118451 (to DJHO) and National Institutes of Health grant OD P51 OD011107 (to CNPRC).

Author information

Authors and Affiliations



Conceptualization, GML, AA, and DJHO. Methodology, NC, NRN, and DJHO. Software, NC. Validation, NC and GML. Formal analysis, NC, NRN, and DJHO. Investigation, NC, NRN, GML, JHF, CC, HTK, and WL. Resources, GML, AA, WLWC, JDD, PAB, EES, and DJHO. Data curation, NC, NRN, GML, and DJHO. Writing — original draft, NC and NRN. Writing — review and editing, NC, PAB, EES, and DJHO. Visualization, NC and DJHO. Supervision, DJHO. Project administration, PAB, EES, and DJHO. Funding acquisition, DJHO. The authors read and approved the final manuscript.

Corresponding author

Correspondence to Dennis J. Hartigan-O’Connor.

Ethics declarations

Ethics approval and consent to participate

This study was performed in compliance with the NIH Guide for the Care and Use of Laboratory Animals. Established policies of the Institutional Animal Care and Use Committee (IACUC) of the University of California, Davis, were followed for sample collections, housing, and medical care. The study was performed at the California National Primate Research Center (CNPRC), which is one of seven centers supported by the National Institutes of Health, Office of the Director (NIH/OD), and is accredited by the Association for the Assessment and Accreditation of Laboratory Animal Care International (AAALAC). The CNPRC houses nonhuman primates, most of which are rhesus macaques (Macaca mulatta), in outdoor corrals or indoor cages.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

Spearman correlations between immune cell subsets and CMV-microbe score.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Chin, N., Narayan, N.R., Méndez-Lagares, G. et al. Cytomegalovirus infection disrupts the influence of short-chain fatty acid producers on Treg/Th17 balance. Microbiome 10, 168 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Host-microbe interactions
  • Microbiome
  • Cytomegalovirus infection
  • Immunophenotype
  • Elastic net
  • Rhesus macaque
  • 16S analysis