Sepsis in preterm infants causes alterations in mucosal gene expression and microbiota profiles compared to non-septic twins

Sepsis is a life-threatening condition in preterm infants. Neonatal microbiota plays a pivotal role in the immune system maturation. Changes in gut microbiota have been associated to inflammatory disorders; however, a link with sepsis in the neonatal period has not yet been established. We aimed to analyze gut microbiota and mucosal gene expression using non-invasively obtained samples to provide with an integrative perspective of host-microbe interactions in neonatal sepsis. For this purpose, a prospective observational case-control study was conducted in septic preterm dizygotic twins and their non-septic twin controls. Fecal samples were used for both microbiota analysis and host genome-wide expression using exfoliated intestinal cells. Gene expression of exfoliated intestinal cells in septic preterm showed an induction of inflammatory and oxidative stress pathways in the gut and pro-oxidant profile that caused dysbiosis in the gut microbiota with predominance of Enterobacteria and reduction of Bacteroides and Bifidobacterium spp.in fecal samples, leading to a global reduction of beneficial anaerobic bacteria. Sepsis in preterm infants induced low-grade inflammation and oxidative stress in the gut mucosa, and also changes in the gut microbiota. This study highlights the role of inflammation and oxidative stress in neonatal sepsis on gut microbial profiles.

homeostasis and consequently the gut microbiota. Moreover, following a septic process preterm infants exhibit a distorted microbiome with predominance of Staphylococcus species and reduced diversity with no specific enrichment of potential pathogens 22,23 .
Genome-wide expression profiles can discriminate septic from non-septic preterm infants in the neonatal period 24 . Gene expression analysis of exfoliated intestinal cells (EIC) and the transcriptional information obtained could disclose non-invasively relevant information about the biologic situation of the intestinal epithelial tissue 16 . However, studies of gene expression in EIC and microbiota in septic preterm infants have not been yet conducted.
The aim of the present study was to get an insight into the processes taking place in the gut of preterm infants during sepsis compared to their non-septic twins searching for possible relationships between changes in the gut microbiota and gene expression of EIC.

Results
Population. Five pairs of preterm twins (≈ 30 weeks' gestation) were enrolled. Each pair included one twin who developed sepsis and a non-septic control. No other differences were observed between cases and controls ( Table 1). Two of the neonates with sepsis had a positive blood culture test. The causal agent was identified as coagulase-negative Staphylococcus strain.
Transcriptomic analysis of exfoliated epithelial cells. Total RNA from the fecal samples of the infants was hybridized with whole human genome microarrays (28,000 annotated genes). Three-dimensional unsupervised principal component analysis (PCA) showed two clustered groups that included the septic and non-septic control samples (Fig. 1a). Further analysis yielded 819 differentially expressed probe sets (ANOVA, p-value < 0.05) between septic and non-septic infants (Tables S1 and S2). Only 510 annotated genes were considered from which 343 genes were up-regulated and 167 were down-regulated. The unsupervised hierarchical clustering showed sensitive differences in gene expression between septic infants and non-septic controls (Fig. 1b). Table 2A shows the most up and down-regulated genes in both groups. Proteins encoded in the list reflect certain degree of inflammation and infiltration of immune cells in the gut mucosa (fecal up-regulation of CD40LG during sepsis). These data are in agreement with gene expression in peripheral blood cells of preterm infants with sepsis 24 (Table 2B).

Functional annotations, upstream regulators and signaling pathways. According to Pathway
Studio (PS) most of the induced gene expression changes were related to regulation of transcription (Table 3A). The most significant were DNA-dependent regulation of transcription (P = 2.53E-08) representing typical interaction modules with DNA and RNA (Table S3), and HDAC2 encoding for histone acetylase 2 25 . The annotation term metabolic process entailed a large number of genes encoding different enzymes involved in ATB-binding (Table S3) and oxidative stress ((glutathione peroxidase 1 (GPX1); glutathione S-transferase mu2 (GSTM2)) 26,27 .
Further, protein-protein networks and signaling pathways could be inferred from expression data. Annotated genes showing significant up-and down-regulation with greater fold changes (FC > 2 and FC < − 2) were selected for the PS analysis. They revealed interaction pathway related to Oxidative Stress (Fig. 2, Panel a) and its connection to the canonical NF-κ B inflammatory pathway (Fig. 2, Panel b). Upstream regulator analysis (Table 3B) identified FOXN4 and SP1 as most significant transcription regulators.
Ingenuity Pathway Analysis (IPA) provided a functional analysis with gene Networks conceptually complementary to the biological processes and master regulators described above (Table 3C). The most significant interaction network related a collection of genes to the molecular functions Cell Death and Survival, Cancer and Inflammatory Response (Fig. 3, Panel a). In agreement with PS, using the upstream regulator analysis tool of the IPA software a contribution to regulatory role could be again attributed to IL1-β , which predicts IL1-β to be activated as transcriptional regulator (z-score 2.368), although, due to the low number of activated genes in the pathway there was a low significance score (P = 0.18) (Fig. 3, Panel b).

Microbial diversity and composition determined by molecular methods. A primary analysis
of the main bacterial groups performed by qPCR of 16S rDNA showed differences between septic and control twins (Fig. 4a). Healthy controls exhibited higher amounts and prevalence of Bifidobacterium species (P = 0.050 and P = 0.055, respectively) and lower levels of Enterobacteriaceae (P = 0.093). Although Enterobacteriaceae, Enterococcus and Staphylococcus groups were common in the septic group, prevalence had low significance scores (P > 0.05). By pyrosequencing, we obtained an average of 10.265,5 (SD ± 4957.04) sequences per infant sample. The average length per sequence was 425 nucleotides (range 300-455). Table S4 shows the predominant phyla present in all samples. Thus while Lactobacillales and Enterobacteriales showed higher frequency in septic patients Bacillales were lower. Pseudomonadales and Vibrionales were found only in the septic group.
We observed that the septic group showed higher bacterial diversity than the control group although differences were not significant (Fig. 4b). Moreover, taxons of twin neonates suffering sepsis included predominantly aerobic or aerotolerant species with a decrease of strict anaerobes. Associations between microbes and neonatal intestinal gene expression. First we analyzed both the sets of genes involved in inflammation and oxidative stress (except SLC19A1 and OR2H2) to establish a correlation between these and the most prominent bacterial clades in septic patients (Fig. 2). Pearson's correlation of the expression of oxidative stress-related genes with significant clades revealed a significant inverse correlation between the order Bacteroidales, family Bacteroidaceae and genus Bacteroides and 8 of 21 genes in the oxidative stress pathway (11 of 21 genes if we include "Other Bacteroidales") ( Fig. 5a

Table 2. Up-and down-regulated genes in preterm infants with sepsis ordered by expression value (Panel A) and Comparison of gene expression (fold change) between data obtained from feces in the present study in preterm infants with sepsis and those from blood cells also in septic preterm newborn described in a previous study (Panel B). (see ref. 24).
Scientific RepoRts | 6:25497 | DOI: 10.1038/srep25497 confirmed with data obtained by qPCR, as a significant inverse correlation with the Bacteroides/Prevotella group was confirmed with 4 of 21 genes. Interestingly, qPCR data also indicated a significant inverse correlation between the calculated reads of the genus Bifidobacterium and the expression of oxidative stress genes (8 of 21) (Fig. 5b, Table S5A). Both genera Bifidobacterium and Bacteroides include strict anaerobic species. Other bacterial groups decreased significantly when signal transduction pathways controlled by the master regulator NF-κ B were active such as the order Mycoplasmatales and the genus Staphylococcus with significant correlation with 6 and 5 genes respectively out of 15 (Fig. 5b, Table S5B). Enterobacteria counts obtained by qPCR directly correlated with a moderate activation of the NF-κ B and IL-1β pathways (6 and 5 of 15 genes, respectively).

Discussion
The analysis of intestinal cells and gut microbiota in the present study offers an integrative perspective of the events occurring in the gut of preterm infants during sepsis considering that inter-individual differences were minimized by the selection of dizygotic twin pairs with and without sepsis. Gene expression differences between septic and non-septic infants were moderate, between 2 and − 2, because EIC integrate pools of various types of cells -epithelial cells, lymphocytes-at different physiological conditions, including anoikis. Despite this variability, our data demonstrated that EIC from septic patients significantly expressed more genes bound to inflammatory pathways and oxidative stress than non-septic controls. Of note, none of the analysis performed showed significant differences between culture positive samples and those with clinical diagnostic.  Some pro-inflammatory genes were up-regulated in EIC revealing the implication of inflammatory regulators; however, these genes do not participate in canonical responses to gut bacteria as pathogen-associated molecular patterns. The infiltration of lymphocytes in feces may partially explain finding of just some differentially expressed genes in exfoliated cells similar to those reported in blood of septic preterm infants 24 (Table 2) as expected in the case of inflammatory organ crosstalk. Genes showing a similar trend in EIC and peripheral blood -like OSCAR,  ARID5A SMYD2 and CD3e-may probably reveal the presence of circulating immune cells in feces. OSCAR encodes a collagen receptor that in monocytes and neutrophils plays a role in the release of pro-inflammatory cytokines and oxidative stress mediators 28 . ARID5A contributes to the stabilization of the mRNA encoding the pro-inflammatory cytokine IL-6 29,30 . SMYD2 encodes a mono-methyl-transferase that methylates acceptor lysine residues on histones 31,32 and other proteins like tumor suppressor proteins p53 and pRb relevant to tumor progression and/or undetermined epigenetic events [33][34][35] . Finally, CD3E down-regulation may indicate that antigen recognition may be uncoupled during sepsis in CD4 + T lymphocytes and T-cells 36 .  Functional annotation analysis offered intriguing clues on possible connections to epigenetic modifications. In septic infants we found activated biological processes like gene transcription that involved down-regulation of HDAC2 a histone acetylase 34,35 and PARP1 37 but also up-regulation of SMYD2 that methylates histones 31,32 , and DNA methyltransferase associated protein (DMAP1). Histone deacetylases can modulate gene expression through endogenous factors, but also through dietary components, synthetic inhibitors, and bacteria-derived signals 38,39 , thereby underlining the mutualistic relationship that exists between gut microbiota and host.
The finding of active oxidative stress pathways in exfoliated cells led us to reanalyze data obtained from peripheral blood cells of VLBW infants, and we found significant values for the oxidative stress pathway (p = 0.0015) that involved the genes BCL2, LCN2, MMP9, EPAS1, BTG1, TOR1 and VNN1 24 . The activation of the oxidative pathways in gut cells could be initiated by the circulating inflammatory mediators or by reactive oxygen (ROS) or nitrogen (RNS) species secreted by blood lymphocytes 40 . Metabolomic analysis has shown that sepsis in the neonatal period is associated to an increase in oxidative stress biomarkers 41 . ROS and RNS intermediates are powerful antimicrobials and major components of the innate defense frontline; however, beyond a threshold they have deleterious effects. Remarkably, oxidative stress has been related to changes in the microbial patterns in the gut 42,43 . Thus, we have observed a decrease of strict anaerobes (Bacteroides, Bifidobacterium) in the septic group, followed by the proliferation of aerobes (Enterobacteriaceae). The genera belonging to this family include species equipped with specialized enzymes such as tetrathionate reductase 44 and other versatile cytochrome related enzymes that allow the use of ROS-oxidized molecules as electron acceptors 45 . Other bacterial groups found in the sepsis samples can reduce nitrate and nitrite to NO, or have catalase or pseudo-catalase enzymes 46 . Altogether these properties would grant them competitive advantage. A general model can be inferred in septic preterm infants that would allow understanding of the likely relationship between blood lymphocyte inflammatory signals during sepsis, with expression of inflammatory and oxidative stress genes, as well as changes in the gut microbiota (Fig. 6). In blood cells, master regulators induce canonical inflammatory pathways 24 and oxidative stress elicited by bacterial septic invasion, leading to the secretion of abundant pro-inflammatory cytokines. Through inflammatory organ crosstalk that led to low but significant expression of inflammatory networks and oxidative stress genes in the gut mucosa -exfoliated epithelial cellsthat drives the secretion of ROS and NOS to the lumen, hence affecting the viability of sensitive bacterial species, like Bacteroides and Bifidobacterium species.
To summarize, this pilot study shows that sepsis in preterm infants induces inflammation in the gut mucosa and changes in gut microbiota extending inflammation-associated damage to sensitive organs. Our results also support the relevance of oxidative stress in sepsis, which may cause the reduction of Bifidobacterium spp. and other beneficial anaerobes in the gut lumen. These data will open new possibilities to reduce the neonatal sepsis using the knowledge and potential role of human breast milk acting as a "gut protector" and antioxidant 47 .

Methods
Patients and Design. We conducted a prospective, observational, cohort study in the NICU of the University and Polytechnic Hospital La Fe from April 2011 to September 2012. Eligible patients were twin pairs with a birth weight ≤ 1.500 g. When one of the twins exhibited clinical signs of sepsis, the other was simultaneously recruited as control if lacking signs of sepsis to homogenize the genetic background. Sepsis was considered clinical signs such as the presence of 3 or more of the following: (i) temperature instability (rectal temp. > 38 °C or < 36 °C); (ii) respiratory symptoms (respiratory distress, apnoea or cyanosis); (iii) cardiovascular symptoms: hypotension (blood pressure < 5th percentile for age), tachycardia (HR > 0180/ min), bradycardia (HR < 100/min) or poor perfusion; (iv) neurological symptoms: clinical or electrical seizures, hypotonia or lethargy; and (v) gastrointestinal symptoms: vomiting, poor feeding or feeding intolerance and/or abdominal distension, without identification of a bacterial pathogen from a sterile site. Exclusion criteria were: maternal immunodeficiency; Apgar 0 at 1 min; chromosomal abnormalities or major gastro-intestinal malformation; congenital infection; concomitant diagnosis of necrotizing enterocolitis.
The clinical interventions included in the study protocol were designed according to the Clinical Guidelines approved by the Spanish Neonatal Society (SENeo) and the Ethical Guidelines of the Comité de Ética e Investigación Médica (CEIM) of the Health Research Institute La Fe. The study protocol was approved by the CEIM. The inform consent was signed by parents of all the enrolled patients.

Diagnosis of sepsis.
A standard blood culture was performed (BacT/Alert ® PF; Biomérieux ® , Durham, NC; USA) only in patients with suspected sepsis. Microbiologic positive diagnosis was considered when a microorganism was isolated from blood culture together with clinical signs and risk factors 48 . Two positive blood cultures were required to diagnose CONS sepsis. Cultures obtained from other sterile sites were not considered for diagnosis of sepsis.
Samples. Stool samples from both twins were collected immediately after sepsis was suspected. Samples were frozen and stored at − 80 °C for later analysis.

DNA extraction and microbial composition by quantitative PCR (qPCR). DNA was extracted from
200-300 mg fecal samples using a modified Qiagen stool DNA extraction kit (QIAgen, Hilden, Germany), and qPCR amplification and detection were performed in a LightCycler ® 480 Real-Time PCR System (Roche) as previously described 50 . The bacterial concentration in each sample was calculated by comparing the Ct values obtained from standard curves. These were created using serial 10-fold dilution of pure culture-specific DNA fragments corresponding to 10 to 10 9 number of gene copies/ml. For statistical analysis, SPSS 15.0 (SPSS Inc, Chicago, IL) software was used. Due to non-normal distribution, microbial data are expressed as medians with interquartile ranges (IQR). The Kruskall-Wallis test was used to compare data of more than two groups of babies, while Mann-Whitney U test was used for comparing data of two groups. The χ -square test was used to establish differences in the prevalence of bacteria between the groups. A p < 0.050 was considered statistically significant. The possible correlation between variables was studied by applying the Pearson correlation coefficient and significance was established at 0.5%.
Microbial diversity and composition estimated by pyrosequencing total 16S rRNA genes in feces. A barcoded primer set based on universal primers 27F and 533R was used to amplify 500 bps of the 16S rRNA genes encompassing the V1-V3 region. PCR was carried out using a high-fidelity KAPA-HiFi polymerase (Kappa Biosystems, US) with an annealing temperature of 52 °C and 30 cycles. Amplicons were checked and measured using the Agilent High Sensitivity DNA assay in Agilent 2100 Expert. Purified PCR products were pyrosequenced from the forward primer end only using a GS-Junior sequencer with Titanium chemistry (Roche) at the Scientific Platform unit SCSIE (University of Valencia, Spain). From the resulting raw data set, low quality sequences or sequences with a length less than 150 nucleotides were discarded and also, chimeric sequences were removed. The QIIME pipeline (version 1.8.0 and green genes data base 13_5) was used for identifying representative sequences for each operational taxonomic unit (OTU) generated from complete linkage clustering with a 97% similarity. Alpha diversity indices were determined from rarefied tables using the Shannon-Wiener index for diversity and the Chao1 index for species richness; Observed Species (number of unique OTUs) and Phylogenetic Distance (PD_whole) were also determined. An OTU was a cluster of 16S rRNA sequences that were > 95% identical, a conservative estimate for the boundary between species, established at 97% (species) and 95% (genus) for full-length 16S gene sequences. A beta diversity distance matrix was computed from the previously constructed OTU table using UniFrac analysis. Unweighted (presence/ absence matrix) and weighted (presence/absence/abundance matrix) UniFrac distances were used to construct two-and three-dimensional Principal Coordinates Analysis (PCoA) plots. Biplots were generated as part of the beta diversity analysis in QIIME, using genus level OTU tables showing principle coordinate sample clustering Scientific RepoRts | 6:25497 | DOI: 10.1038/srep25497 alongside weighted taxonomic group data. Data on assigned sequences at genus level shared between samples were used to generate a Venn diagram [51][52][53][54][55][56][57][58][59][60] .