Abstract
The microbiota plays a critical role in regulating organismal health and response to environmental stresses. Intermittent hypoxia and hypercapnia, a condition that represents the main hallmark of obstructive sleep apnea in humans, is known to induce significant alterations in the gut microbiome and metabolism, and promotes the progression of atherosclerosis in mouse models. To further understand the role of the microbiome in the cardiovascular response to intermittent hypoxia and hypercapnia, we developed a new rodent cage system that allows exposure of mice to controlled levels of O2 and CO2 under gnotobiotic conditions. Using this experimental setup, we determined the impact of the microbiome on the transcriptional response to intermittent hypoxia and hypercapnia in the left ventricle of the mouse heart. We identified significant changes in gene expression in both conventionally reared and germ-free mice. Following intermittent hypoxia and hypercapnia exposure, we detected 192 significant changes in conventionally reared mice (96 upregulated and 96 downregulated) and 161 significant changes (70 upregulated and 91 downregulated) in germ-free mice. Only 19 of these differentially expressed transcripts (∼10%) were common to conventionally reared and germ-free mice. Such distinct transcriptional responses imply that the host microbiota plays an important role in regulating the host transcriptional response to intermittent hypoxia and hypercapnia in the mouse heart.
Introduction
Obstructive sleep apnea (OSA) is the most common type of sleep disorder in humans, and OSA is very prevalent in adults and children. It is characterized by repeated episodes of complete or partial obstructions of the upper airways during sleep, which is associated with a reduction in blood oxygen saturation and an increase in blood CO2. It has been demonstrated by experimental and clinical studies that OSA may lead to a number of health consequences including cardiovascular diseases ranging from cardiometabolic disorders and arrhythmogenesis to heart failure (Al Lawati et al., 2009; Baguet et al., 2012; Levy et al., 2013; Drager et al., 2015; Lyons and Bradley, 2015; Hoyos et al., 2017; May et al., 2017).
In addition, evidence from both animal and human studies supports the notion that the gut microbiota plays a critical role in the onset and progression of various cardiovascular diseases (Karbach et al., 2016; Koch et al., 2017; Weis, 2018). Previous studies by others and us have shown that intermittent hypoxia with/without hypercapnia alters the gut microbiome (Moreno-Indias et al., 2015, 2016; Lucking et al., 2018; Tripathi et al., 2018, 2019). In turn, such alterations may induce developmental defects and functional dysregulations. For example, the human oral and intestinal microbiota has been recognized as a significant determinant of cardiovascular disease risk (Ordovas and Mooser, 2006), and suppression of the gut microbiome ameliorates age-dependent arterial dysfunction (Brunt et al., 2019). Recent studies also revealed a potential contribution of gut microbes to certain human cardiometabolic diseases (Vinje et al., 2014) and heart failure (Sandek et al., 2007; Nagatomo and Tang, 2015; Nemet et al., 2020).
Besides impacting cellular signaling, the microbiome has been demonstrated to modify the epigenetic landscape of the host genome, which may influence cellular differentiation and proliferation, leading to new phenotypic traits throughout the age spectrum (Takahashi et al., 2011; Pevsner-Fischer et al., 2017; Qin and Wade, 2018). For example, several studies have demonstrated a direct link between the gut microbiome and host epigenetic states in the context of inflammatory bowel disease, obesity, cancer, and behavioral abnormalities (Kellermayer et al., 2010; Manichanh et al., 2012; Dejea et al., 2014; Paul et al., 2015; Farhana et al., 2018; Sfanos et al., 2018; Chu et al., 2019).
Since gut microbiota plays an important role in regulating host gene expression (Patrignani et al., 2014; Chu et al., 2019; Richards et al., 2019; Weger et al., 2019), in the current study, we tested the hypothesis that the microbiome regulates the transcriptional response to intermittent hypoxia and hypercapnia (IHH) exposure in the heart. Using a newly developed experimental setup to impose IHH under gnotobiotic conditions, we performed 150-bp paired-end RNA-seq and utilized multiple bioinformatic approaches to analyze the transcriptional alterations following IHH treatment in the left ventricle (LV) of conventionally reared (CONV) and germ-free (GF) C57BL/6 mice. In addition, we compared our results with those reported for other organs in male GF mice of similar genetic background and age (Weger et al., 2019) (dataset GSE77221). We found significant and distinct differences in transcriptional responses between CONV and GF mice, demonstrating that, indeed, the microbiome plays an important role in regulating IHH response in the heart.
Materials and Methods
Experimental Animals
GF and CONV C57BL/6 mice that were 8–12 weeks old were used for this experiment. GF status was confirmed weekly by standard microbiological cultures of fecal extracts, using thioglycollate broth and potato dextrose broth for facultative anaerobes including fungi, and cooked meat broth, heart infusion broth, and Columbia blood agar for obligate anaerobes (Moody et al., 2019). All animal protocols were approved by the Animal Care Committee of the University of California, San Diego and followed the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health.
IHH Treatment
All mice were housed in specialized cages for the 2-week treatments (IHH or normoxia; GF or CONV) (Figure 1). These specialized cages were developed by us in the Gnotobiotic Animal Facility at the University of California, San Diego. Each sealed single cage (Sealed Positive Pressure System, Allentown, NJ, United States) was ventilated by a single, HEPA (high-efficiency particulate air)-filtered gas supply pump. The cages were opened and closed for stocking and servicing inside a Class IIa laminar-flow biosafety cabinet. The system allows complete control of gas inflow and outflow, as has been confirmed by measurements of exhaust gases using an O2/CO2 analyzer (Quantek Instruments, Grafton, MA, United States). The IHH condition was maintained in a computer-controlled atmosphere chamber system (OxyCycler, BioSpherix, Redfield, NY, United States) as previously described (Fan et al., 2005; Zhou et al., 2008). IHH exposure was introduced to the mice in short periods of ∼4 min of synchronized reduction in (O2) (from 21 to 8%) and increase in (CO2) (from ∼0.5 to 8%) separated by alternating ∼4-min periods of normoxic room air [(O2) = 21% and (CO2) ≈ 0.5%] with 2 ramp intervals (∼1 min each) for a total of 10 IHH events per hour for 10 h per day during the light cycle. Controls received a flow of room air throughout the entire experimental period. Prior works have shown that the oxygen saturation level (SaO2) ranged from 60 to 80% when mice were subjected to 5–10% hypoxic treatments (Sforza and Roche, 2016; Xue et al., 2017).
FIGURE 1
RNA Sequencing and Data Analysis
At the end of the IHH exposure, mice were sacrificed, and the heart was harvested and separated into the right and left ventricle under sterile conditions and then placed into separate tubes containing RNAlater Stabilization Solution (Thermo Fisher Scientific, Ashville, NC, United States). RNA from the LV was extracted using the RNeasy Mini Kit (QIAGEN, Hilden, Germany) and sent to Novogene (Chula Vista, CA, United States) for QC and sequencing (n = 3/group). A 150 nt paired-end sequencing on the Illumina HiSeq 4000 platform (Illumina, San Diego, CA, United States) was performed. Rigorous quality controls of the paired-end reads were assessed using FastQC tools, and the adapter sequences and low-quality bases were trimmed using Cutadapt (Didion et al., 2017). The resulting reads were processed using the Illumina BaseSpace software package using the STAR aligner (Dobin et al., 2013) to align reads against the Mus musculus UCSC mm9 reference genome. Differential expression was calculated using DESeq 2 with default settings (Anders and Huber, 2010), and differentially expressed genes (DEGs) were identified with a threshold of |Fold-Change| > 1.5 and adjusted p-value < 0.05. Bioinformatic data mining of the differentially expressed transcripts was carried out with the tools provided by the GWAS Catalog1 (Buniello et al., 2019), the Database for Annotation, Visualization, and Integrated Discovery (DAVID)2 (Dennis et al., 2003; Huang et al., 2007), the Ingenuity Pathway Analysis (IPA) package (QIAGEN, Redwood City, CA, United States), as well as KEGG3 (Ogata et al., 1999; Kanehisa and Goto, 2000) and the Reactome4 (Fabregat et al., 2017). Multiple testing correction was applied, and the FDR q-values were estimated to correct the p-values for the multiple testing issue.
Data Availability
The paired-end RNA-seq data of n = 12 samples are available at https://www.ncbi.nlm.nih.gov/bioproject/692551 (BioProject ID: PRJNA692551).
Results and Discussion
Experimental Design
The microbiome, particularly in the gut, plays an important role in maintaining systemic homeostasis and contributes to a number of human diseases [for selected reviews, see Lynch and Pedersen (2016), Gilbert et al. (2018), Dominguez-Bello et al. (2019), Harikrishnan (2019), Helmink et al. (2019), Team (2019)]. However, its potential involvement in regulating many physiologic and pathophysiologic responses under stressful conditions remains to be explored. In the current study, we used a prospective case-control study design to determine the impact of the microbiome on cardiac transcriptional responses to IHH in CONV and GF mice (Figure 1A). To maintain GF animals under controlled gas conditions, we developed a custom cage system with minimal dead space and programmable HEPA-filtered gas supply (Figure 1B). Analysis of the exhaust gases of the cage system revealed rapid (<20 s) equilibration of the cage environment to match the input gas mixture (Figure 1C). GF mice could be kept in the system for extended periods (>2 weeks) without compromising their GF state, demonstrating that the newly developed system was well suited to conduct short- and long-term experiments of constant or intermittent hypoxia with or without hypercapnia under gnotobiotic conditions.
GF and CONV mice of the same strain, sex, and age (kept in the same cage system as the GF mice) were either exposed to the IHH condition for 2 weeks to mimic sleep apnea in humans (Kanaan et al., 2006; Xue et al., 2017; Hunyor and Cook, 2018) or kept under normoxia/normocapnia throughout as controls. For a broad characterization of the cardiac gene expression responses, we extracted total RNA from the LV and applied paired-end RNA-seq as a genome-wide discovery approach. Approximately 38–50 million total reads were obtained per sample with over 94% of the reads uniquely mapped to the mouse reference genome (UCSC mm9) (Supplementary Table 1).
Microbiota Impacts Cardiac Gene Expression Under Normoxic Conditions
We first determined the impact of the microbiota on cardiac gene expression under normoxic and normocapnic conditions. A total of 117 DEGs were identified between CONV and GF mice. Of these, 73 were upregulated and 44 were downregulated (Fold change > 1.5 and adjusted p-value < 0.05) in the left ventricle of GF mice as compared to CONV controls (Figure 2A and Supplementary Table 2). The most upregulated transcripts included Amd1, Zfp950, Rnf125, and Khdc4, which are involved in polyamine biosynthesis, transcription regulation, ubiquitination, and subsequent proteasomal degradation of target proteins, and pre-mRNA splicing, respectively. In contrast, the most downregulated transcripts included genes encoding a transcription repressor (Zmynd11), a Calcium Voltage-Gated Channel Subunit (Cacna1d), and a serine/threonine kinase (Limk1) (Figure 2B). In addition, several functional implications were suggested by data mining and bioinformatic analyses using the DAVID and IPA tools. Results derived from DAVID functional annotation suggested alterations in transcription regulation, alpha-actinin binding, and sarcolemma organization, as well as pathways that regulate cardiac function and activity, including the GnRH pathway, the MAPK signaling pathway, the Oxytocin signaling pathway, and vascular muscle contraction (Supplementary Table 3A). In addition, Ingenuity analysis showed potential changes in several canonical pathways that are important in cardiac function, such as the mechanisms regulating synaptic long-term potentiation, Myc-mediated apoptosis signaling, and the role of NFAT in cardiac hypertrophy (Figure 3 and Supplementary Table 3B). The different alterations on the “Myc Mediated Apoptosis Signaling” and the “Role of NFAT in Cardiac Hypertrophy” suggested that cell death and proliferation signaling pathways were differently configured in the GF mice (Supplementary Table 3B). Since (1) Myc-mediated apoptosis plays a critical role in response to stress in the myocardium (Hou et al., 2016; Veeraveedu et al., 2017; Guan et al., 2019), (2) NFAT plays a major role in cardiac hypertrophy [for selected reviews see Clerk et al. (2007), Rohini et al. (2010), Nakayama et al. (2013)], and (3) previous studies also shown that gut microbiome regulates the activity of Myc in mouse liver (Kolodziejczyk et al., 2020) and NFAT in colorectal cancer in humans (Peuker et al., 2016), our current findings suggest that alterations of Myc and NFAT signaling in GF mice may modify cardiac activity and stress response.
FIGURE 2
FIGURE 3
Furthermore, we found that 36% of the DEGs (42 out of 117) in GF mice encode nuclear proteins, including transcription factors (e.g., Dach1 and Gtf2ird2), protein kinases (e.g., Clk1 and Limk1), DNA, and histone modification enzymes (e.g., Dnmt3b and Kmt5b). Such a significant enrichment of nuclear proteins suggested differences of transcriptional regulation mechanisms between GF and CONV mice, which may lead to distinct transcriptional responses to stress in the heart.
Distinct Transcriptional Alterations in Different Organs of GF Mice
Alterations in gene expression in GF mice have been reported in several organs (Weger et al., 2019; Mills et al., 2020). In order to determine common and distinct transcriptional responses between different organs in GF mice, we compared the reported changes (GSE77221, GEO, NCBI) with those in our study. As shown in Figure 4, each organ displayed a distinct transcriptional profile. Among the total 5,109 significantly altered transcripts (|Fold Change| > 1.5 and FDR < 0.05) in the heart, liver, white adipose tissue, duodenum, and ileum of GF mice, 3,386 (66%) were tissue specific. The largest number of unique changes was observed in the ileum (1,144 genes), followed by duodenum (1,096 genes), liver (587 genes), white fat (467 genes), and heart (92 genes), possibly reflecting the relative vicinity of the respective organs to the abundant gut microbiota. No genes were identified that were altered in all five organs, indicating a lack of common tissue responses to the microbiota. Although the vast majority of cardiac genes with significant expression changes were unique to the heart under GF conditions (79%), the most overlap was seen with the duodenum (5 genes, 4.3% of total DEGs in the heart) and the least with duodenum (or ileum) and white fat (1 gene each, <1% of total DEGs in the heart). Interestingly, the gene cytochrome b-245, beta chain (Cybb) was commonly altered in heart and gut, suggesting that the lack of a microbiome alters the activity of immune cells, especially neutrophils, which may lead to immune dysregulation in these organs.
FIGURE 4
GF Condition Alters Transcriptional Responses to IHH in the Heart
Previous studies have shown that stress can alter the composition of the gut microbiome, and in turn, such changes in gut microbiota may regulate stress response in various organs and tissues, including the central nervous system, the digestive system, and the immune system (Lyte et al., 2020; Park et al., 2020; Rengarajan et al., 2020; Tian et al., 2020; Vuong et al., 2020; Wang et al., 2020; Xu et al., 2020). The current study provides evidence that the microbiome may also affect stress responses in the heart through regulating transcription activities. Indeed, following 2 weeks of IHH exposure, we detected 192 significant changes in CONV mice (96 upregulated and 96 downregulated; | Fold Change| > 1.5 and FDR < 0.05) (Figure 5A and Supplementary Table 4) and 161 significant changes (70 upregulated and 91 downregulated; | Fold Change| > 1.5 and FDR < 0.05) in GF mice (Figure 5B and Supplementary Table 5). As shown in Table 1, the top DEGs in CONV mice include Rnf44, Cyb561a3, Cep250, Asph, Gm14296, and Epb41l1, whereas Kmt5a, Amd1, Tnnt2, Slc44a2, Golga5, and Cpeb3 were among the top changes in GF mice. Interestingly, only 19 transcripts were altered in both CONV and GF mice, representing around 10% of the changes in CONV and GF, respectively, demonstrating that IHH treatment induced distinct and specific transcriptional responses between CONV and GF mice. The majority of these common altered transcripts were downregulated, and the top three downregulated transcripts were Acta1 (Alpha actin 1), Rn45s (45S pre-ribosomal RNA), and Rhobtb1 (Rho-related BTB domain containing 1). On the other hand, Tfrc (transferrin receptor), Rgs2 (regulator of G-protein signaling 2), and Prkar1a (protein kinase, cAMP-dependent regulatory, type I, alpha) were the top three genes that were commonly upregulated in both CONV and GF mice. In contrast, two transcripts, Nadk2 (mitochondrial NAD kinase 2) and Ssbp1 (single-stranded DNA binding protein 1), were downregulated in CONV mice but upregulated in GF mice following IHH exposure (Figure 5C). Such dramatically different transcriptional responses suggested that IHH exposure may lead to different physiological and functional alterations in the heart of CONV and GF mice.
FIGURE 5
TABLE 1
| RefSeq_ID | Symbol | Normalized read counts | Fold change | Padj | |||||
| Air-1 | Air-2 | Air-3 | IHH-1 | IHH-2 | IHH-3 | ||||
| CONV Mice | |||||||||
| NR_027395 | Rnf44 | 272.4 | 260.7 | 211.7 | 0.0 | 0.0 | 0.0 | −1406.2 | 6.5E-15 |
| NM_201351 | Cyb561a3 | 158.4 | 127.1 | 150.7 | 0.0 | 0.0 | 0.0 | −827.7 | 1.8E-12 |
| NM_008383 | Cep250 | 141.4 | 111.9 | 138.3 | 0.0 | 0.0 | 0.0 | −740.1 | 6.2E-12 |
| NM_001286664 | Ssbp1 | 93.6 | 88.9 | 82.6 | 0.0 | 0.0 | 0.0 | −509.8 | 2.8E-10 |
| NM_001199248 | Kat5 | 123.6 | 107.7 | 17.6 | 0.0 | 0.0 | 0.0 | −471.9 | 8.8E-07 |
| NM_001163009 | Tnik | 62.4 | 73.3 | 97.8 | 0.0 | 0.0 | 0.0 | −443.1 | 2.6E-09 |
| NM_001204277 | Hdac7 | 81.2 | 89.8 | 54.1 | 0.0 | 0.0 | 0.0 | −422.7 | 5.2E-09 |
| NR_130970_1 | Gm14296 | 71.3 | 48.6 | 79.5 | 0.0 | 0.0 | 0.0 | −381.8 | 1.3E-08 |
| NM_008488 | Arhgef1 | 27.4 | 81.0 | 85.5 | 0.0 | 0.0 | 0.0 | −365.4 | 3.4E-07 |
| NM_008642 | Mttp | 61.8 | 46.2 | 74.7 | 0.0 | 0.0 | 0.0 | −349.2 | 3.0E-08 |
| NM_001347201 | Gucd1 | 0.0 | 0.0 | 0.0 | 14.4 | 28.5 | 43.9 | 150.7 | 2.3E-04 |
| NM_001012363 | Slc2a9 | 0.0 | 0.0 | 0.0 | 28.6 | 39.1 | 24.6 | 159.8 | 2.7E-05 |
| NM_001358355 | Zfx | 0.0 | 0.0 | 0.0 | 34.6 | 50.0 | 27.8 | 191.7 | 7.3E-06 |
| NM_178722 | Zfp438 | 0.0 | 0.0 | 0.0 | 38.2 | 42.6 | 55.3 | 234.8 | 6.2E-07 |
| NM_001167886 | Kmt5b | 0.0 | 0.0 | 0.0 | 48.3 | 70.0 | 21.8 | 241.4 | 7.3E-06 |
| NM_001284506 | Plxnb2 | 0.0 | 0.0 | 0.0 | 40.1 | 70.8 | 40.4 | 263.5 | 5.6E-07 |
| NM_001270495 | Tmem254a | 0.0 | 0.0 | 0.0 | 51.6 | 54.4 | 48.3 | 266.2 | 1.3E-07 |
| NM_001291122 | Epb41l1 | 0.0 | 0.0 | 0.0 | 100.8 | 47.7 | 33.0 | 314.3 | 8.4E-07 |
| NR_130970 | Gm14296 | 0.0 | 0.0 | 0.0 | 67.2 | 61.0 | 85.7 | 367.4 | 8.8E-09 |
| NM_001290367 | Asph | 0.0 | 0.0 | 0.0 | 168.1 | 41.5 | 83.9 | 502.3 | 3.9E-08 |
| GF Mice | |||||||||
| NM_001310725 | Kmt5a | 332.6 | 616.9 | 544.7 | 0.0 | 0.0 | 0.0 | −2582.6 | 1.9E-17 |
| NM_009665 | Amd1 | 321.6 | 421.3 | 295.8 | 0.0 | 0.0 | 0.0 | −1800.0 | 3.6E-16 |
| NM_001130180 | Tnnt2 | 0.0 | 69299.9 | 70470.5 | 61.6 | 13.2 | 12.4 | −1602.1 | 2.1E-03 |
| NM_134125 | Trip10 | 92.3 | 146.7 | 489.2 | 0.0 | 0.0 | 0.0 | −1251.8 | 1.0E-07 |
| NM_001102438 | Acbd5 | 107.8 | 132.9 | 183.7 | 0.0 | 0.0 | 0.0 | −734.7 | 1.1E-11 |
| NM_001360080 | Herc2 | 151.3 | 143.0 | 85.8 | 0.0 | 0.0 | 0.0 | −657.8 | 4.6E-11 |
| NM_001329143 | Adamts10 | 94.8 | 120.1 | 117.7 | 0.0 | 0.0 | 0.0 | −571.0 | 4.8E-11 |
| NM_001361696 | Slc38a10 | 170.4 | 58.5 | 88.2 | 10.3 | 29.6 | 12.4 | −6.1 | 9.7E-04 |
| NM_001039087 | Rapgef1 | 767.1 | 518.2 | 501.2 | 156.1 | 160.6 | 283.6 | −3.0 | 2.4E-06 |
| NM_001243837 | C7 | 290.4 | 258.7 | 163.4 | 80.0 | 95.4 | 69.2 | −2.9 | 3.2E-06 |
| NM_011638 | Tfrc | 1034.7 | 887.6 | 1268.0 | 2834.2 | 3037.6 | 3108.3 | 2.8 | 1.7E-28 |
| NM_001304810 | Kcnj3 | 37.4 | 11.6 | 29.9 | 79.4 | 115.0 | 66.8 | 3.3 | 1.7E-02 |
| NM_001277328 | Sox6 | 57.1 | 31.0 | 27.7 | 163.1 | 117.7 | 123.1 | 3.5 | 3.5E-05 |
| NM_001357298 | Tfrc | 11.3 | 65.4 | 25.6 | 136.5 | 133.8 | 147.1 | 4.1 | 1.3E-02 |
| NM_001313974 | Prkar1a | 26.0 | 24.4 | 23.2 | 256.7 | 24.7 | 170.0 | 6.1 | 1.1E-02 |
| NM_001355067 | Ctnnd1 | 0.0 | 0.0 | 0.0 | 145.9 | 152.4 | 18.6 | 600.8 | 1.7E-07 |
| NM_001199308 | Gm14440 | 0.0 | 0.0 | 0.0 | 133.6 | 156.4 | 63.1 | 670.1 | 2.0E-10 |
| NM_198300 | Cpeb3 | 0.0 | 0.0 | 0.0 | 185.3 | 184.8 | 41.5 | 785.4 | 1.5E-09 |
| NM_013747 | Golga5 | 0.0 | 0.0 | 0.0 | 151.5 | 128.8 | 200.1 | 912.2 | 9.2E-13 |
| NM_001199186 | Slc44a2 | 0.0 | 0.0 | 0.0 | 222.6 | 265.2 | 248.5 | 1400.1 | 3.3E-15 |
Top transcriptional alterations in CONV and GF mice following IHH exposure.
Distinct Alterations in Cell Signaling in Response to IHH Exposure Between CONV and GF Mice
We identified distinct cell signaling mechanisms in the hearts of CONV and GF mice in response to IHH. In the CONV mice, the top significantly altered transcripts encoding a group of nuclear proteins [66/192 (34%), FDR < 0.05] that are involved in “DNA-templated transcription” [GO:0006351, 42/192 (22%), FDR < 0.05] and “regulation of DNA-templated transcription” [GO:0006355, 46/192 (24%), FDR < 0.05] (Supplementary Table 6A). The top altered signaling mechanisms were Myc mediated apoptosis signaling, coronavirus pathogenesis pathway, and IL-9 signaling (Figure 6A and Supplementary Table 6B). In GF mice, on the other hand, the significantly altered genes were involved in the biological processes that regulate cell cycle (GO:0051726, including some cyclins, e.g., Cdkn1a, Ccnd1, and Cables1, FDR < 0.05), biological rhythms (FDR < 0.05), and MHC class II protein complex (GO:0042613, FDR < 0.05) (Supplementary Table 7A). The top representative pathways were “regulation of actin-based motility by Rho,” “integrin signaling,” “HGF signaling,” and “circadian rhythm signaling” (Figure 6B and Supplementary Table 7B).
FIGURE 6
Following IHH exposure, “regulation of DNA-templated transcription” was one of the most significantly altered biological processes represented by 24% of the significantly altered transcripts in CONV mice. These results suggested a fine-tuning of the transcriptome, at least in part, through epigenetic modifications and rebalance of transcription regulators. Indeed, as shown in Figure 7, transcripts encoding several epigenetic modifiers were significantly changed in the CONV mice, but not GF mice, following IHH treatment. In contrast, in the heart of the GF mouse, IHH-induced transcriptional alterations suggested an array of functional/activity adjustments, such as the circadian rhythm and cell cycle genes (e.g., Per2, Clock, Cdkn1a, and Stat3) that regulate the heart rate (Yamamura et al., 2010) or the metabolism, hypertrophy, and apoptosis of cardiac cells (Bonney et al., 2013; Hauck et al., 2016; Harhous et al., 2019).
FIGURE 7
Transcriptomic Indications of Potential Cardiac Disorders in GF and CONV Mice Following IHH Treatment
Substantial evidence has been reported that OSA significantly increases the risk of stroke, coronary heart disease, and heart failure in humans (McEvoy et al., 2016; Drager and Lee, 2018; Benjamin et al., 2019). Interestingly, we found that the IHH-induced transcriptional alterations in the GF and CONV mice suggested distinct risks of cardiac disorders, i.e., IHH exposure evoked different sets of well-characterized pathogenetic genes in GF and CONV mice. As shown in Table 2 and Figure 8, in CONV mice, IHH-induced transcriptional alterations mainly affect cardiac necrosis and cell death. In contrast, in the GF mice, the DEGs included several key regulators of cardiac hypertrophy. In the CONV mice, changes of both well-characterized cell death stimulators and inhibitors were detected, which include BCL2-like 1 (Bcl2l1), Crystallin alpha B (Cryab), Gelsolin (Gsn), Inhibitor of nuclear factor kappa B kinase subunit beta (Ikbkb), Nuclear receptor subfamily 4 group A member 1 (Nr4a1), Protein phosphatase 1 regulatory subunit 10 (Ppp1r10), and Signal transducer and activator of transcription 1 (Stat1), as well as the downregulation of Thioredoxin interacting protein (Txnip) [for selected reviews see Wagner and Siddiqui (2009), Li et al. (2012), Ganguly et al. (2014), Wang and Yoshioka (2017), Opferman and Kothari (2018), Herring et al. (2019), Tang et al. (2019)]. In the GF mice, however, most of the disease-related DEGs were cardiac hypertrophy regulators, such as the downregulation of Angiotensin I converting enzyme (Ace), Ankyrin repeat domain 1 (Ankrd1), Apelin receptor (Aplnr), Nitric oxide synthase 2 (Nos2), and Cardiac type Troponin T2 (Tnnt2), as well as the upregulation of Cyclin-dependent kinase inhibitor 1A (Cdkn1a), Four and a half LIM domains 2 (Fhl2), Regulator of G protein signaling 2 (Rgs2), and Signal transducer and activator of transcription 3 (Stat3) [for selected reviews, see Brooks et al. (1998), Drexler (1999), Takeda et al. (2005), Tsang et al. (2010), Haghikia et al. (2011), Liang et al. (2018), Yotti et al. (2019)]. Since several genes in this group are key regulators of hypertension (e.g., Ace, Nos2, and Rgs2) (Li and Forstermann, 2000; Prichard et al., 2001; Tsang et al., 2010), it is reasonable to hypothesize that dysregulation of blood pressure might be an important mechanism regulating IHH-induced hypertrophy in GF mice.
TABLE 2
![]() |
Transcriptome indication of potential IHH-induced cardiac disorders in GF and CONV mice.
Red font: Upregulation. Blue font: Downregulation.
FIGURE 8
In summary, our study demonstrated that the microbiome regulates the cardiac transcriptome at baseline and in response to IHH stress. A total of 117 DEGs were identified in the heart of GF mice under normoxic condition. These DEFs encode proteins involved in transcriptional regulation, alpha-actinin binding, and sarcolemma organization. Lack of microbiome also altered cardiac transcriptional responses to IHH stress, since we found 192 IHH-induced DEGs in CONV mice and 161 DEGs in GF mice with little overlap, suggesting that IHH exposure leads to different cardiac responses between CONV and GF mice. Indeed, 24% of the DEGs detected in CONV mice are involved in the “regulation of DNA-templated transcription,” while in GF mice, IHH-induced cardiac DEGs are mostly involved in functional or activity adjustments, such as circadian rhythm and cell cycle. Future functional studies will be needed to define the functional implications of the observed DEGs. For example, IHH-induced DEGs suggest distinct risks of cardiac disorders in GF and CONV mice, since disease-related DEGs mainly affected cardiac necrosis and cell death in CONV mice, while regulators of cardiac hypertrophy were mostly observed in GF mice. Additional studies are also warranted to explore the function of specific microbial species and their metabolites on regulating cardiac gene expression under baseline and stress conditions. Gnotobiotic mouse models are uniquely suited to unravel such microbial contributions to specific physiologic responses.
Conclusion
Even though it has been shown that OSA or IHH exposure may lead to the development of various types of cardiac disorders [for selected reviews, see Baguet et al. (2012), Dumitrascu et al. (2013), Levy et al. (2015), Turnbull (2018)], any contribution of the microbiome to this process is still largely uncharacterized. In the present study, we found that the microbiome plays an important role in the transcriptional responses to IHH stress and identified distinct responsive mechanisms in the heart of GF and CONV mice, suggesting that the microbiome modulates the cardiac function during stressful and disease conditions like OSA. Our study provides insight into the unique and common transcriptional response to IHH that is related to microbial colonization in mice. These transcriptional alterations suggest changes in cellular signaling and interaction networks regulating cardiac function, such as cardiac hypertrophy and cell death. Overall, our study is a step toward a better understanding of the role of microbiota in regulating cardiac health and disease.
Statements
Data availability statement
The paired-end RNA-seq data of n = 12 samples are available at https://www.ncbi.nlm.nih.gov/bioproject/692551 (BioProject ID: PRJNA692551).
Ethics statement
The animal study was reviewed and approved by the Animal Care Committee of the University of California, San Diego.
Author contributions
DZ, JX, LE, and GH conceptualized the study, designed the experiments, and wrote the manuscript. DZ, JX, YM, and OP performed the experiments. DZ and JX performed data curation and analysis. LE and GH acquired funding and facilitated access to resources for the study. GH supervised the study. All authors contributed to the article and approved the submitted version.
Funding
This work was partially supported by the NIH grant P30 DK120515 (LE).
Acknowledgments
We thank Iain Hartley and Travis Smith for technical support.
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2021.680275/full#supplementary-material
Abbreviations
- OSA
obstructive sleep apnea
- IHH
intermittent hypoxia and hypercapnia
- GF
germ-free
- CONV
conventional
- LV
left ventricle
- HEPA
high-efficiency particulate air
- DAVID
Database for Annotation Visualization and Integrated Discovery
- IPA
Ingenuity Pathway Analysis
- KEGG
Kyoto Encyclopedia of Genes and Genomes.
References
1
Al LawatiN. M.PatelS. R.AyasN. T. (2009). Epidemiology, risk factors, and consequences of obstructive sleep apnea and short sleep duration.Prog. Cardiovasc. Dis.51285–293. 10.1016/j.pcad.2008.08.001
2
AndersS.HuberW. (2010). Differential expression analysis for sequence count data.Genome Biol.11:R106. 10.1186/gb-2010-11-10-r106
3
BaguetJ. P.Barone-RochetteG.TamisierR.LevyP.PepinJ. L. (2012). Mechanisms of cardiac dysfunction in obstructive sleep apnea.Nat. Rev. Cardiol.9679–688. 10.1038/nrcardio.2012.141
4
BenjaminE. J.MuntnerP.AlonsoA.BittencourtM. S.CallawayC. W.CarsonA. P.et al (2019). Heart Disease and Stroke Statistics-2019 Update: a Report From the American Heart Association.Circulation139e56–e528. 10.1161/CIR.0000000000000659
5
BonneyS.KominskyD.BrodskyK.EltzschigH.WalkerL.EckleT. (2013). Cardiac Per2 functions as novel link between fatty acid metabolism and myocardial inflammation during ischemia and reperfusion injury of the heart.PLoS One8:e71493. 10.1371/journal.pone.0071493
6
BrooksG.PoolmanR. A.LiJ. M. (1998). Arresting developments in the cardiac myocyte cell cycle: role of cyclin-dependent kinase inhibitors.Cardiovasc. Res.39301–311. 10.1016/s0008-6363(98)00125-4
7
BruntV. E.Gioscia-RyanR. A.RicheyJ. J.ZiglerM. C.CuevasL. M.GonzalezA.et al (2019). Suppression of the gut microbiome ameliorates age-related arterial dysfunction and oxidative stress in mice.J. Physiol.5972361–2378. 10.1113/JP277336
8
BunielloA.MacarthurJ. A. L.CerezoM.HarrisL. W.HayhurstJ.MalangoneC.et al (2019). The NHGRI-EBI GWAS Catalog of published genome-wide association studies, targeted arrays and summary statistics 2019.Nucleic Acids Res.47D1005–D1012. 10.1093/nar/gky1120
9
ChuC.MurdockM. H.JingD.WonT. H.ChungH.KresselA. M.et al (2019). The microbiota regulate neuronal function and fear extinction learning.Nature574543–548. 10.1038/s41586-019-1644-y
10
ClerkA.CullingfordT. E.FullerS. J.GiraldoA.MarkouT.PikkarainenS.et al (2007). Signaling pathways mediating cardiac myocyte gene expression in physiological and stress responses.J. Cell Physiol.212311–322. 10.1002/jcp.21094
11
DejeaC. M.WickE. C.HechenbleiknerE. M.WhiteJ. R.Mark WelchJ. L.RossettiB. J.et al (2014). Microbiota organization is a distinct feature of proximal colorectal cancers.Proc. Natl. Acad. Sci. U. S. A.11118321–18326. 10.1073/pnas.1406199111
12
DennisG.Jr.ShermanB. T.HosackD. A.YangJ.GaoW.LaneH. C.et al (2003). DAVID: database for Annotation, Visualization, and Integrated Discovery.Genome Biol.4:P3.
13
DidionJ. P.MartinM.CollinsF. S. (2017). Atropos: specific, sensitive, and speedy trimming of sequencing reads.PeerJ5:e3720. 10.7717/peerj.3720
14
DobinA.DavisC. A.SchlesingerF.DrenkowJ.ZaleskiC.JhaS.et al (2013). STAR: ultrafast universal RNA-seq aligner.Bioinformatics2915–21. 10.1093/bioinformatics/bts635
15
Dominguez-BelloM. G.Godoy-VitorinoF.KnightR.BlaserM. J. (2019). Role of the microbiome in human development.Gut681108–1114. 10.1136/gutjnl-2018-317503
16
DragerL. F.LeeC. H. (2018). Treatment of obstructive sleep apnoea as primary or secondary prevention of cardiovascular disease: where do we stand now?Curr. Opin. Pulm. Med.24537–542. 10.1097/MCP.0000000000000523
17
DragerL. F.PolotskyV. Y.O’donnellC. P.CravoS. L.Lorenzi-FilhoG.MachadoB. H. (2015). Translational approaches to understanding metabolic dysfunction and cardiovascular consequences of obstructive sleep apnea.Am. J. Physiol. Heart Circ. Physiol.309H1101–H1111. 10.1152/ajpheart.00094.2015
18
DrexlerH. (1999). Nitric oxide synthases in the failing human heart: a doubled-edged sword?Circulation992972–2975. 10.1161/01.cir.99.23.2972
19
DumitrascuR.HeitmannJ.SeegerW.WeissmannN.SchulzR. (2013). Obstructive sleep apnea, oxidative stress and cardiovascular disease: lessons from animal studies.Oxid. Med. Cell. Longev.2013:234631. 10.1155/2013/234631
20
FabregatA.SidiropoulosK.ViteriG.FornerO.Marin-GarciaP.ArnauV.et al (2017). Reactome pathway analysis: a high-performance in-memory approach.BMC Bioinformatics18:142. 10.1186/s12859-017-1559-2
21
FanC.IacobasD. A.ZhouD.ChenQ.LaiJ. K.GavrialovO.et al (2005). Gene expression and phenotypic characterization of mouse heart after chronic constant or intermittent hypoxia.Physiol. Genomics22292–307. 10.1152/physiolgenomics.00217.2004
22
FarhanaL.BanerjeeH. N.VermaM.MajumdarA. P. N. (2018). Role of Microbiome in Carcinogenesis Process and Epigenetic Regulation of Colorectal Cancer.Methods Mol. Biol.185635–55. 10.1007/978-1-4939-8751-1_3
23
GangulyS.MitraA.SarkarS. (2014). Role of alpha-crystallin B in regulation of stress induced cardiomyocyte apoptosis.Cardiovasc. Hematol. Agents Med. Chem.1260–65. 10.2174/1871525713666150123151731
24
GilbertJ. A.BlaserM. J.CaporasoJ. G.JanssonJ. K.LynchS. V.KnightR. (2018). Current understanding of the human microbiome.Nat. Med.24392–400. 10.1038/nm.4517
25
GuanF.YangX.LiJ.DongW.ZhangX.LiuN.et al (2019). New Molecular Mechanism Underlying Myc-Mediated Cytochrome P450 2E1 Upregulation in Apoptosis and Energy Metabolism in the Myocardium.J. Am. Heart Assoc.8:e009871. 10.1161/JAHA.118.009871
26
HaghikiaA.StapelB.HochM.Hilfiker-KleinerD. (2011). STAT3 and cardiac remodeling.Heart Fail. Rev.1635–47. 10.1007/s10741-010-9170-x
27
HarhousZ.BoozG. W.OvizeM.BidauxG.KurdiM. (2019). An Update on the Multifaceted Roles of STAT3 in the Heart.Front. Cardiovasc. Med.6:150. 10.3389/fcvm.2019.00150
28
HarikrishnanS. (2019). Diet, the Gut Microbiome and Heart Failure.Card. Fail. Rev.5119–122. 10.15420/cfr.2018.39.2
29
HauckL.GrotheD.BilliaF. (2016). p21(CIP1/WAF1)-dependent inhibition of cardiac hypertrophy in response to Angiotensin II involves Akt/Myc and pRb signaling.Peptides8338–48. 10.1016/j.peptides.2016.07.003
30
HelminkB. A.KhanM. A. W.HermannA.GopalakrishnanV.WargoJ. A. (2019). The microbiome, cancer, and cancer therapy.Nat. Med.25377–388. 10.1038/s41591-019-0377-7
31
HerringJ. A.ElisonW. S.TessemJ. S. (2019). Function of Nr4a Orphan Nuclear Receptors in Proliferation, Apoptosis and Fuel Utilization Across Tissues.Cells8:1373. 10.3390/cells8111373
32
HouN.YeB.LiX.MarguliesK. B.XuH.WangX.et al (2016). Transcription Factor 7-like 2 Mediates Canonical Wnt/beta-Catenin Signaling and c-Myc Upregulation in Heart Failure.Circ. Heart Fail.9:e003010. 10.1161/CIRCHEARTFAILURE.116.003010
33
HoyosC. M.DragerL. F.PatelS. R. (2017). OSA and cardiometabolic risk: what’s the bottom line?Respirology22420–429. 10.1111/resp.12984
34
HuangD. W.ShermanB. T.TanQ.CollinsJ. R.AlvordW. G.RoayaeiJ.et al (2007). The DAVID Gene Functional Classification Tool: a novel biological module-centric algorithm to functionally analyze large gene lists.Genome Biol.8:R183. 10.1186/gb-2007-8-9-r183
35
HunyorI.CookK. M. (2018). Models of intermittent hypoxia and obstructive sleep apnea: molecular pathways and their contribution to cancer.Am. J. Physiol. Regul. Integr. Comp. Physiol.315R669–R687. 10.1152/ajpregu.00036.2018
36
KanaanA.FarahaniR.DouglasR. M.LamannaJ. C.HaddadG. G. (2006). Effect of chronic continuous or intermittent hypoxia and reoxygenation on cerebral capillary density and myelination.Am. J. Physiol. Regul. Integr. Comp. Physiol.290R1105–R1114. 10.1152/ajpregu.00535.2005
37
KanehisaM.GotoS. (2000). KEGG: kyoto encyclopedia of genes and genomes.Nucleic Acids Res.2827–30. 10.1093/nar/28.1.27
38
KarbachS. H.SchonfelderT.BrandaoI.WilmsE.HormannN.JackelS.et al (2016). Gut Microbiota Promote Angiotensin II-Induced Arterial Hypertension and Vascular Dysfunction.J. Am. Heart Assoc.5:e003698. 10.1161/JAHA.116.003698
39
KellermayerR.BalasaA.ZhangW.LeeS.MirzaS.ChakravartyA.et al (2010). Epigenetic maturation in colonic mucosa continues beyond infancy in mice.Hum. Mol. Genet.192168–2176. 10.1093/hmg/ddq095
40
KochC. D.GladwinM. T.FreemanB. A.LundbergJ. O.WeitzbergE.MorrisA. (2017). Enterosalivary nitrate metabolism and the microbiome: intersection of microbial metabolism, nitric oxide and diet in cardiac and pulmonary vascular health.Free Radic. Biol. Med.10548–67. 10.1016/j.freeradbiomed.2016.12.015
41
KolodziejczykA. A.FedericiS.ZmoraN.MohapatraG.Dori-BachashM.HornsteinS.et al (2020). Acute liver failure is regulated by MYC- and microbiome-dependent programs.Nat. Med.261899–1911. 10.1038/s41591-020-1102-2
42
LevyP.KohlerM.McnicholasW. T.BarbeF.McevoyR. D.SomersV. K.et al (2015). Obstructive sleep apnoea syndrome.Nat. Rev. Dis. Primers1:15015. 10.1038/nrdp.2015.15
43
LevyP.RyanS.OldenburgO.ParatiG. (2013). Sleep apnoea and the heart.Eur. Respir. Rev.22333–352. 10.1183/09059180.00004513
44
LiG. H.AroraP. D.ChenY.MccullochC. A.LiuP. (2012). Multifunctional roles of gelsolin in health and diseases.Med. Res. Rev.32999–1025. 10.1002/med.20231
45
LiH.ForstermannU. (2000). Nitric oxide in the pathogenesis of vascular disease.J. Pathol.190244–254. 10.1002/(SICI)1096-9896(200002)190:3<244::AID-PATH575<3.0.CO;2-8
46
LiangY.BradfordW. H.ZhangJ.SheikhF. (2018). Four and a half LIM domain protein signaling and cardiomyopathy.Biophys. Rev.101073–1085. 10.1007/s12551-018-0434-3
47
LuckingE. F.O’connorK. M.StrainC. R.FouhyF.BastiaanssenT. F. S.BurnsD. P.et al (2018). Chronic intermittent hypoxia disrupts cardiorespiratory homeostasis and gut microbiota composition in adult male guinea-pigs.EBioMedicine38191–205. 10.1016/j.ebiom.2018.11.010
48
LynchS. V.PedersenO. (2016). The Human Intestinal Microbiome in Health and Disease.N. Engl. J. Med.3752369–2379. 10.1056/NEJMra1600266
49
LyonsO. D.BradleyT. D. (2015). Heart Failure and Sleep Apnea.Can. J. Cardiol.31898–908. 10.1016/j.cjca.2015.04.017
50
LyteJ. M.GheorgheC. E.GoodsonM. S.Kelley-LoughnaneN.DinanT. G.CryanJ. F.et al (2020). Gut-brain axis serotonergic responses to acute stress exposure are microbiome-dependent.Neurogastroenterol. Motil.32:e13881. 10.1111/nmo.13881
51
ManichanhC.BorruelN.CasellasF.GuarnerF. (2012). The gut microbiota in IBD.Nat. Rev. Gastroenterol. Hepatol.9599–608. 10.1038/nrgastro.2012.152
52
MayA. M.Van WagonerD. R.MehraR. (2017). OSA and Cardiac Arrhythmogenesis: mechanistic Insights.Chest151225–241. 10.1016/j.chest.2016.09.014
53
McEvoyR. D.AnticN. A.HeeleyE.LuoY.OuQ.ZhangX.et al (2016). CPAP for Prevention of Cardiovascular Events in Obstructive Sleep Apnea.N. Engl. J. Med.375919–931. 10.1056/NEJMoa1606599
54
MillsR. H.WozniakJ. M.VrbanacA.CampeauA.ChassaingB.GewirtzA.et al (2020). Organ-level protein networks as a reference for the host effects of the microbiome.Genome Res.30276–286. 10.1101/gr.256875.119
55
MoodyL. V.MiyamotoY.AngJ.RichterP. J.EckmannL. (2019). Evaluation of Peroxides and Chlorine Oxides as Disinfectants for Chemical Sterilization of Gnotobiotic Rodent Isolators.J. Am. Assoc. Lab. Anim. Sci.58558–568. 10.30802/AALAS-JAALAS-18-000130
56
Moreno-IndiasI.TorresM.MontserratJ. M.Sanchez-AlcoholadoL.CardonaF.TinahonesF. J.et al (2015). Intermittent hypoxia alters gut microbiota diversity in a mouse model of sleep apnoea.Eur. Respir. J.451055–1065. 10.1183/09031936.00184314
57
Moreno-IndiasI.TorresM.Sanchez-AlcoholadoL.CardonaF.AlmendrosI.GozalD.et al (2016). Normoxic Recovery Mimicking Treatment of Sleep Apnea Does Not Reverse Intermittent Hypoxia-Induced Bacterial Dysbiosis and Low-Grade Endotoxemia in Mice.Sleep391891–1897. 10.5665/sleep.6176
58
NagatomoY.TangW. H. (2015). Intersections Between Microbiome and Heart Failure: revisiting the Gut Hypothesis.J. Card. Fail.21973–980. 10.1016/j.cardfail.2015.09.017
59
NakayamaH.FujioY.YamaguchiO. (2013). [Calcium dependent signaling in cardiac hypertrophy and cell death].Clin. Calcium23505–515. CliCa1304505515,
60
NemetI.SahaP. P.GuptaN.ZhuW.RomanoK. A.SkyeS. M.et al (2020). A Cardiovascular Disease-Linked Gut Microbial Metabolite Acts via Adrenergic Receptors.Cell180862–877.e22. 10.1016/j.cell.2020.02.016
61
OgataH.GotoS.SatoK.FujibuchiW.BonoH.KanehisaM. (1999). KEGG: kyoto Encyclopedia of Genes and Genomes.Nucleic Acids Res.2729–34. 10.1093/nar/27.1.29
62
OpfermanJ. T.KothariA. (2018). Anti-apoptotic BCL-2 family members in development.Cell Death Differ.2537–45. 10.1038/cdd.2017.170
63
OrdovasJ. M.MooserV. (2006). Metagenomics: the role of the microbiome in cardiovascular diseases.Curr. Opin. Lipidol.17157–161. 10.1097/01.mol.0000217897.75068.ba
64
ParkY. S.KimS. H.ParkJ. W.KhoY.SeokP. R.ShinJ. H.et al (2020). Melatonin in the colon modulates intestinal microbiota in response to stress and sleep deprivation.Intest. Res.18325–336. 10.5217/ir.2019.00093
65
PatrignaniP.TacconelliS.BrunoA. (2014). Gut microbiota, host gene expression, and aging.J. Clin. Gastroenterol.48S28–S31. 10.1097/MCG.0000000000000229
66
PaulB.BarnesS.Demark-WahnefriedW.MorrowC.SalvadorC.SkibolaC.et al (2015). Influences of diet and the gut microbiome on epigenetic modulation in cancer and other diseases.Clin. Epigenetics7:112. 10.1186/s13148-015-0144-7
67
PeukerK.MuffS.WangJ.KunzelS.BosseE.ZeissigY.et al (2016). Epithelial calcineurin controls microbiota-dependent intestinal tumor development.Nat. Med.22506–515. 10.1038/nm.4072
68
Pevsner-FischerM.ZmoraN.BravermanS.ElinavE. (2017). “Epigenetics and the Microbiome,” in Handbook of Nutrition, Diet, and Epigenetics, edsPatelV.PreedyV. (Cham: Springer International Publishing), 1–25. 10.1007/978-3-319-31143-2_35-1
69
PrichardB. N.CruickshankJ. M.GrahamB. R. (2001). Beta-adrenergic blocking drugs in the treatment of hypertension.Blood Press.10366–386. 10.1080/080370501753400665
70
QinY.WadeP. A. (2018). Crosstalk between the microbiome and epigenome: messages from bugs.J. Biochem.163105–112. 10.1093/jb/mvx080
71
RengarajanS.KnoopK. A.RengarajanA.ChaiJ. N.Grajales-ReyesJ. G.SamineniV. K.et al (2020). A Potential Role for Stress-Induced Microbial Alterations in IgA-Associated Irritable Bowel Syndrome with Diarrhea.Cell Rep. Med.1:100124. 10.1016/j.xcrm.2020.100124
72
RichardsA. L.MuehlbauerA. L.AlaziziA.BurnsM. B.FindleyA.MessinaF.et al (2019). Gut Microbiota Has a Widespread and Modifiable Effect on Host Gene Regulation.mSystems4:e00323. 10.1128/mSystems.00323-18
73
RohiniA.AgrawalN.KoyaniC. N.SinghR. (2010). Molecular targets and regulators of cardiac hypertrophy.Pharmacol. Res.61269–280. 10.1016/j.phrs.2009.11.012
74
SandekA.BauditzJ.SwidsinskiA.BuhnerS.Weber-EibelJ.Von HaehlingS.et al (2007). Altered intestinal function in patients with chronic heart failure.J. Am. Coll. Cardiol.501561–1569. 10.1016/j.jacc.2007.07.016
75
SfanosK. S.YegnasubramanianS.NelsonW. G.De MarzoA. M. (2018). The inflammatory microenvironment and microbiome in prostate cancer development.Nat. Rev. Urol.1511–24. 10.1038/nrurol.2017.167
76
SforzaE.RocheF. (2016). Chronic intermittent hypoxia and obstructive sleep apnea: an experimental and clinical approach.Hypoxia499–108. 10.2147/HP.S103091
77
TakahashiK.SugiY.NakanoK.TsudaM.KuriharaK.HosonoA.et al (2011). Epigenetic control of the host gene by commensal bacteria in large intestinal epithelial cells.J. Biol. Chem.28635755–35762. 10.1074/jbc.M111.271007
78
TakedaT.OtsuK.HoriM. (2005). [Calcium channel blocker and cardiac hypertrophy].Clin. Calcium151645–1649. CliCa051016451649,
79
TangD.KangR.BergheT. V.VandenabeeleP.KroemerG. (2019). The molecular machinery of regulated cell death.Cell Res.29347–364. 10.1038/s41422-019-0164-5
80
TeamN. I. H. H. M. P. A. (2019). A review of 10 years of human microbiome research activities at the US National Institutes of Health, Fiscal Years 2007-2016.Microbiome7:31. 10.1186/s40168-019-0620-y
81
TianP.O’riordanK. J.LeeY. K.WangG.ZhaoJ.ZhangH.et al (2020). Towards a psychobiotic therapy for depression: bifidobacterium breve CCFM1025 reverses chronic stress-induced depressive symptoms and gut microbial abnormalities in mice.Neurobiol. Stress12:100216. 10.1016/j.ynstr.2020.100216
82
TripathiA.MelnikA. V.XueJ.PoulsenO.MeehanM. J.HumphreyG.et al (2018). Intermittent Hypoxia and Hypercapnia, a Hallmark of Obstructive Sleep Apnea, Alters the Gut Microbiome and Metabolome.mSystems3e00020–18. 10.1128/mSystems.00020-18
83
TripathiA.XuZ. Z.XueJ.PoulsenO.GonzalezA.HumphreyG.et al (2019). Intermittent Hypoxia and Hypercapnia Reproducibly Change the Gut Microbiome and Metabolome across Rodent Model Systems.mSystems4e00058–19. 10.1128/mSystems.00058-19
84
TsangS.WooA. Y.ZhuW.XiaoR. P. (2010). Deregulation of RGS2 in cardiovascular diseases.Front. Biosci.2:547–557. 10.2741/s84
85
TurnbullC. D. (2018). Intermittent hypoxia, cardiovascular disease and obstructive sleep apnoea.J. Thorac. Dis.10S33–S39. 10.21037/jtd.2017.10.33
86
VeeraveeduP. T.SanadaS.OkudaK.FuH. Y.MatsuzakiT.ArakiR.et al (2017). Ablation of IL-33 gene exacerbate myocardial remodeling in mice with heart failure induced by mechanical stress.Biochem. Pharmacol.13873–80. 10.1016/j.bcp.2017.04.022
87
VinjeS.StroesE.NieuwdorpM.HazenS. L. (2014). The gut microbiome as novel cardio-metabolic target: the time has come!Eur. Heart J.35883–887. 10.1093/eurheartj/eht467
88
VuongH. E.PronovostG. N.WilliamsD. W.ColeyE. J. L.SieglerE. L.QiuA.et al (2020). The maternal microbiome modulates fetal neurodevelopment in mice.Nature586281–286. 10.1038/s41586-020-2745-3
89
WagnerM.SiddiquiM. A. (2009). Signaling networks regulating cardiac myocyte survival and death.Curr. Opin. Investig. Drugs10928–937.
90
WangB. F.YoshiokaJ. (2017). The Emerging Role of Thioredoxin-Interacting Protein in Myocardial Ischemia/Reperfusion Injury.J. Cardiovasc. Pharmacol. Ther.22219–229. 10.1177/1074248416675731
91
WangR. X.LeeJ. S.CampbellE. L.ColganS. P. (2020). Microbiota-derived butyrate dynamically regulates intestinal homeostasis through regulation of actin-associated protein synaptopodin.Proc. Natl. Acad. Sci. U. S. A.11711648–11657. 10.1073/pnas.1917597117
92
WegerB. D.GobetC.YeungJ.MartinE.JimenezS.BetriseyB.et al (2019). The Mouse Microbiome Is Required for Sex-Specific Diurnal Rhythms of Gene Expression and Metabolism.Cell Metab.29362–382.e8. 10.1016/j.cmet.2018.09.023
93
WeisM. (2018). Impact of the gut microbiome in cardiovascular and autoimmune diseases.Clin. Sci.1322387–2389. 10.1042/CS20180410
94
XuC.LeeS. K.ZhangD.FrenetteP. S. (2020). The Gut Microbiome Regulates Psychological-Stress-Induced Inflammation.Immunity53417–428.e4. 10.1016/j.immuni.2020.06.025
95
XueJ.ZhouD.PoulsenO.ImamuraT.HsiaoY. H.SmithT. H.et al (2017). Intermittent Hypoxia and Hypercapnia Accelerate Atherosclerosis, Partially via Trimethylamine-Oxide.Am. J. Respir. Cell Mol. Biol.57581–588. 10.1165/rcmb.2017-0086OC
96
YamamuraY.YanoI.KudoT.ShibataS. (2010). Time-dependent inhibitory effect of lipopolysaccharide injection on Per1 and Per2 gene expression in the mouse heart and liver.Chronobiol. Int.27213–232. 10.3109/07420521003769111
97
YottiR.SeidmanC. E.SeidmanJ. G. (2019). Advances in the Genetic Basis and Pathogenesis of Sarcomere Cardiomyopathies.Annu. Rev. Genomics Hum. Genet.20129–153. 10.1146/annurev-genom-083118-015306
98
ZhouD.WangJ.ZapalaM. A.XueJ.SchorkN. J.HaddadG. G. (2008). Gene expression in mouse brain following chronic hypoxia: role of sarcospan in glial cell death.Physiol. Genomics32370–379. 10.1152/physiolgenomics.00147.2007
Summary
Keywords
gnotobiotics, transcriptome, intermittent hypoxia and hypercapnia, heart, mice
Citation
Zhou D, Xue J, Miyamoto Y, Poulsen O, Eckmann L and Haddad GG (2021) Microbiota Modulates Cardiac Transcriptional Responses to Intermittent Hypoxia and Hypercapnia. Front. Physiol. 12:680275. doi: 10.3389/fphys.2021.680275
Received
13 March 2021
Accepted
20 May 2021
Published
25 June 2021
Volume
12 - 2021
Edited by
Rodrigo Iturriaga, Pontificia Universidad Católica de Chile, Chile
Reviewed by
Lena Lavie, Technion – Israel Institute of Technology, Israel; Bin Liu, University of Texas MD Anderson Cancer Center, United States
Updates
Copyright
© 2021 Zhou, Xue, Miyamoto, Poulsen, Eckmann and Haddad.
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: Gabriel G. Haddad, ghaddad@health.ucsd.edu
†These authors have contributed equally to this work and share first authorship
This article was submitted to Integrative Physiology, a section of the journal Frontiers in Physiology
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.
