|
|
||||||||
1 Istituto di Anestesiologia e Rianimazione, Università degli Studi, Polo San Paolo, via Di Rudinì 8, I-20142 Milano, Italy
2 Centro Ricerche Biotecnologiche, Università Cattolica del Sacro Cuore, Via Milano 24, 26100 Cremona, Italy
3 Istituto di Microbiologia, Università Cattolica del Sacro Cuore, Via E. Parmense 84, 29100 Piacenza, Italy
Correspondence
Lorenzo Morelli
lorenzo.morelli{at}unicatt.it
Received 9 May 2007
Accepted 7 April 2008
Abbreviations: cLOS, critical length of stay; DGGE, denaturing gradient gel electrophoresis; ICU, intensive care unit; LOS, length of stay; SAPSII, simplified acute physiology score II; SOFA, sequential organ failure analysis; TGGE, temperature gradient gel electrophoresis.
| INTRODUCTION |
|---|
|
|
|---|
A number of factors that predispose the gut ecology to alterations can occur during the care of critically ill patients. Those that have been described include: luminal hypoxia and hypercarbia; the use of agents to neutralize gastric secretions; vasoactive agents; sedation and analgesia that may impair intestinal motility; nutrient scarcity within the colon due to parenteral or enteral diets poor in non-digestible carbohydrates (Schneider et al., 2000); and severe sepsis/injury or digestive surgical stress (Marshall, 1999; Alverdy et al., 2003). It is thus highly possible that disruptions to the ecological balance of the intestinal microbial community could occur in critically ill patients undergoing antibiotic treatments.
It is difficult to draw a full and accurate picture of the microbiota in the intestinal tract, as the gut contains more than 600 different bacterial species, most of which cannot be cultured under laboratory conditions (Guarner, 2006). However, culture-independent techniques to investigate microbial gut ecology are now available (Tannock, 2002) and allow better monitoring of these intestinal communities and also of population shifts. Despite the potential relevance of the disruption of the gut microbial ecology to the final clinical outcomes in critically ill patients, there are few studies on this topic (Alverdy et al., 2003; Shimizu et al., 2006). The aim of our current observational study was to assess, by means of molecular techniques such as denaturing gradient gel electrophoresis (DGGE), faecal microbiota changes in high-risk, critically ill patients during their intensive care unit (ICU) stay and treatment with antibiotics, according to the standard ICU protocols.
| METHODS |
|---|
|
|
|---|
Microbiota analysis. For each patient, faeces were collected at ICU admission (T0) and then weekly (T1, T2, etc.) from a rectal sinus, thus avoiding collection by enema, and were stored at –20 °C until DNA analysis. Patients were excluded if the T0 or T1 samples could not be harvested. The presence of Clostridium difficile toxin in the diarrhoeal stools was also assayed using ImmunoCard Toxins A&B (Meridian Bioscience). Total bacterial DNA extractions were performed with 200 mg faecal material using the PSP Spin Stool DNA kit (Invitek) according to the manufacturer's instructions. Target DNA was amplified using commercially available Megamix (Microzone) and a Perkin Elmer thermal cycler 9700. PCR amplifications of total bacterial DNA were obtained using the primer pair HDA1-GC and HDA2, and the thermocycling program as described previously (Walter et al., 2000).
DGGE analysis was performed with a DGene System (Bio-Rad Hercules) using a 0.8 mm 6 % polyacrylamide gel with a ratio of acrylamide to bisacrylamide of 37.5 : 1 and containing a 35–55 % gradient of urea and formamide, as reported previously (Walter et al., 2000). The electrophoretic run was performed in 1x Tris/acetate/EDTA buffer (TAE) at 120 V for 5 h. DGGE gels were stained for 30 min with SYBR Green I (Roche Applied Science) in 1x TAE buffer. The DGGE profiles obtained for each sampling were then compared within each gel using Fingerprint II software (Bio-Rad Hercules). These analyses included the calculation of the number of bands for each pattern and comparisons between patterns using the Pearson coefficient as a measure of the degree of similarity (Seksik et al., 2003). Two identical profiles create a similarity value of 100 % whereas completely different profiles result in a similarity value of 0 %.
To identify the bacteria producing the DNA banding pattern by DGGE at the species level, bands were excised from the stained gels with a sterile scalpel and eluted overnight in 20 µl sterile water at 4 °C, as described by Ben Omar & Ampe (2000). One microlitre of the eluted DNA from each band was then used as the template for a PCR reaction under the conditions described above. The resulting PCR products were purified using a Wizard SV Gel and PCR clean-up system (Promega). The final eluents, containing the purified PCR fragments, were sequenced by BMR Bio Molecular Research (CRIBI, Padova, Italy). The closest known relatives to the partial 16S rRNA gene sequences obtained were identified using the BLAST and RDP programs.
The percentage similarities between the DGGE profiles at T0 and at the subsequent sampling times for each patient were used as an index of microbiota modifications. A value equal or higher than the median value of similarity for all of the T1 profiles, compared with the baseline, was defined as the index of microbiota biodiversity preservation. This index value, coupled with the appearance in some samples of a new dominant DNA band (which was identified as Enterococcus spp.), was then evaluated against clinical variables.
Quantification of Enterococcus by real-time PCR. We quantified the presence of Enterococcus genomic DNA in the patient faecal samples by real-time PCR. The reactions were performed using the LightCycler System and LightCycler FastStart DNA MasterPLUS SYBR Green I (Roche Applied Science). The fragment of the 16S rRNA gene specific for the Enterococcus genus was amplified using the primers ECF (AGAAATTCCAAACGAACTTG) and ECR (CAGTGCTCTACCTTCCATCATT) (He & Jiang, 2005). The 20 µl total reaction volumes contained 4 µl LightCycler FastStart DNA MasterPLUS SYBR Green I (Roche Diagnostics), 0.5 µM forward and reverse primers and 5 µl extracted DNA. The amplification program included an initial denaturation step at 95 °C for 10 min, followed by 40 cycles of denaturation for 15 s at 95 °C, annealing for 20 s at 63 °C, and extension for 30 s at 72 °C. The temperature transition rate was 20 °C s–1.
For melting curve analysis, the samples were denatured at 65 °C for 15 s, and then slowly heated to 95 °C with a temperature transition rate of 0.1 °C s–1 and a continuous monitoring of fluorescence. Enterococcus faecium DNA (BAA-472D-5) was purchased from the American Type Culture Collection and used to optimize the real-time reactions and generate quantification standards. For the construction of standard curves, a 10-fold dilution series ranging from 0.01 pg to 10 ng of the target species genomic DNA was used. The standard curves had correlation coefficient values of between 0.99 and 1.0, whereas the amplification slope was between –3.35 and –3.38. The standard curves for the real-time PCR assays were used to quantify the target bacterial DNA from faecal DNA preparations and PCR results were converted to a mean estimate of the target bacterial genomes present in 1 g faeces (wet wt). The estimated mean genome size for Enterococcus that was used for converting the results obtained by real-time PCR was 3 Mb. Three parallel PCR reactions for each sample were analysed. The mean amounts and standard deviations of the target genomes in 1 g faeces, calculated by PCR, are shown in Table 3
.
|
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
|
|
In our current study, we performed DGGE analysis of faecal samples in all of our patients at ICU admission and found, except for patient 14, a complex banding pattern (Fig. 1a, b
) which included a number of bands previously described in healthy adults using the same technique (Seksik et al., 2003). The similarity between the DGGE profiles in each patient when comparing the T1 versus T0 samples was 62.06±27.11 %, median 68 % (Table 2
). By univariate analysis, the variation in the faecal microbiota at T1 was related to the degree of organ failure (SOFA: R2=0.43, P=0.011) and complexity of ICU treatment (cLOS: R2=0.38, P=0.015). The multivariate equation that we used to describe these relationships was: percentage similarity of DGGE profiles=139.9–3.8 SOFA (T1)–8.8 cLOS (T1) (R2=0.60, P=0.007). The median of the similarity of the DGGE profiles after one week was used to divide the cases into two groups with remarkably different clinical outcomes; the seven patients showing major modifications of their DGGE profiles had a higher SOFA score at T1 (5.4±2.6 vs 2.5±2.2; P=0.029), longer period of invasive ICU support (cLOS: 24.6±25.2, 16, 10–28 vs 7.9±4.1, 8, 4.5–11; P=0.0174) and ICU LOS (25.1±25.9, 16, 11–28 vs 10.3±3.3, 10.5, 7.5–12.0; P=0.0315), coupled with a higher mortality (4 out of 7 vs 0 out of 8; P=0.026). Moreover, although not statistically significant, the patients with greater DGGE profile changes received more anti-anaerobe drugs (71.4 % vs 37.5 %).
|
The second interesting finding from these analyses was that only 1 week of ICU treatment could cause the appearance of a very intense DGGE band which was faint or absent in the T0 profiles of the same subjects (Fig. 1
). Sequence analysis of these bands, excised from the gels for all patients showing this variation, subsequently revealed a significant (more than 95 %) similarity with the sequences of the 16S rRNA gene of the Enterococcus genus. At T2, this Enterococcus band was still present in five patients out of the seven sampled.
As DGGE can only provide a qualitative estimate of the levels of a specific bacterial group, we used RT-PCR to perform more quantitative analysis of bacteria of interest (based on the DNA bands that were detected). The same approach was used previously by Rinttila et al. (2004), who showed that RT-PCR could result in an overestimation when compared with plate counting. The results that we present herein probably suffer from the same overestimation rate, but the relevant outcomes are the comparisons between the sampling points. Rinttila et al. (2004) employed RT-PCR for the quantification of Enterococcus from faecal samples of healthy volunteers and found that the number of bacterial genomes in 1 g faeces ranged from 106 to 107. The levels of Enterococcus target genomes (g faeces)–1 calculated in our present study are presented in Table 3
. Nine of the patients considered in this study exceeded the upper limit.
The increase in the presence of Enterococcus bacteria among the gut faecal microbiota of our patients during their whole ICU stay could only be associated with the overall disruption to the microbiota by univariate analysis (% of similarity of DGGE profiles vs T0 <68 %, P=0.006) and with the clindamycin treatments (P=0.006) that were administered in cases of ab ingestis pneumonia prophylaxis, peritonitis and severe cellulitis. The role of clindamycin in promoting overgrowth of enterococci has been previously shown by Jernberg et al. (2005). It should be noted also that Enterococcus is an emerging genus in ICU epidemiology (Weber & Gold, 2003).
The loss of intestinal bacteria, particularly anaerobes, can induce an overgrowth of Enterococcus, a normal component of the human microbiota (Stiefel et al., 2004; Donskey, 2004). Moreover, a decrease in the number of gut anaerobes impairs the colonic degradation of undigested carbohydrates, provided via enteral nutrition, to absorbable lactic acid and short-chain fatty acids. This leads to an increased osmotic load. Reduced availability of butyrate, an important source of energy for the colonic mucosa (Sheppach et al., 1997), causes functional disorders (Topping & Clifton, 2001) resulting in diarrhoea (Young & Schmidt, 2004). Interestingly, as an indirect sign of severe alteration of bacterial gut ecology, 80 % of the patients in our current cohort experienced at least one diarrhoea episode during enteral nutrition and antibiotic treatment.
Since an evaluation of this group of patients as a whole would make it difficult to draw any conclusions, we adopted a case by case approach and focused our attention on individual patients whose index of microbiota biodiversity preservation was under the cut-off (Tables 1
, 2
and 3
; Fig. 1
). We considered patients 5, 6 and 10 to be characterized by a minor ecological alteration (50–63 % similarity, DGGE profile vs basal profile) after a 1 week ICU stay. Bands later identified as Enterococcus were found to be present in each of these cases at a later stage in the ICU. Patients 5 and 6 had been admitted for severe respiratory failure due to pneumonia and later died. Both were responsive to intensive treatment as shown by their clinical condition coupled with a stable or improved SOFA score at T1. After T1, patient 5 worsened (SOFA score increased to 7) and died of septic shock before T2. Of note in this regard, enterococci were shown to have massively accumulated from T0 to T1 in this patient. Patient 6 remained clinically and ecologically stable until T2, and then had a brisk septic episode (positive blood culture for Staphylococcus aureus) with huge irreversible impairment of the lung and the development of cardiovascular conditions (SOFA 11). The elevated levels of Enterococcus, as detected by RT-PCR, occurred from T0 to T1 and approached a 2-log increase.
In contrast, patient 10 was admitted for a severe diffuse bronchitis and an initial cardiac failure but survived. This individual showed a discrepancy between a very high severity and low organ failure score. This is characteristic of the bias that can arise with severity scores that give excess weight to old age and comorbidities. This male patient showed greater ecological alterations than the two deceased patients (a 50 % similarity and just over the upper normal limit for Enterococcus content), but also showed consistent improvements in his clinical condition (SOFA T1=2); he was transferred to the main wards of the hospital within 2 weeks of admittance to the ICU and before the collection of the T2 sample.
A further scenario that we investigated involved patients 2, 7, 11 and 14, who showed huge biota alterations with the percentage similarity of profiles observed well below the median value, and a massive Enterococcus content after 1 week in the ICU. Patients 2 and 7, who were admitted for pneumonia with severe respiratory failure, had high severity and organ failure scores and subsequently died. The clinical course of the illness in both cases was a slow progressive worsening, as shown by their high SOFA scores at the sampling times and just before death (7 for both). Accordingly, the alterations in the DGGE profiles and enterococci appearance remained very large, even during the second week for patient 7, who died after T2.
Patients 11 and 14 showed the same reduced biota profile at T1 as patients 2 and 7, but survived. Patient 11 was a young man who developed septic shock as a result of extensive cellulitis from the right knee to the axilla. At admission, he showed a high severity (young, no comorbidities) and a huge organ failure score (14, about a 60 % predicted probability of ICU death). The clinical condition of this patient progressively improved, however, and the SOFA decreased. In addition, the number of detectable bands on the DGGE gel increased from T1 to T2, and the enterococci levels showed an increase from the basal values. This patient was transferred to the main wards after 18 days in the ICU. Patient 14 was a 66-year-old woman admitted to our ICU for septic shock due to bilateral Pneumococcus pneumonia after 7 days of neglected illness at home. This was associated with a reduction in the number of DGGE bands at ICU admission, suggesting that her gut microbiota was already altered, and the enterococci levels were calculated to be 8x107 genomes (g faeces)–1. The severe respiratory insufficiency of this patient was further complicated by mediastinitis, which resulted from a tracheal lesion during a difficult rescue intubation. The extensive lung consolidation and a long period of pulmonary septic capillary leak resulted in a refractory pulmonary hypertension that required 156 days of intensive/subintensive stay before a referral to another hospital for rehabilitation (Giacomini et al., 2007). The combination of severe lung infection, mediastinitis with systemic sepsis and full antibiotic treatment over 1 week was associated with an almost total disappearance of gut microbiota, which to a large extent was due to the outgrowth of enterococci. This critical, septic phase lasted for many of the ICU days for this patient and we collected faecal samples for 5 weeks and then on day 58, when the clinical conditions seemed to advise a positive outcome. However, the DGGE profile of the gut microbiota of patient 14 showed a slow progressive increase in the number of bands from a T1 value of similarity of 3.8 % to T3 (17.9 %: SOFA 7), T4 (18.5 %: SOFA 7), day 39 (16.8 %: SOFA 5) and day 58 (43.0 %: SOFA 7). The enterococci counts remained in the high range at T1 and T2 (6.6E+9, 1.7E+10, 1.1E+10, 2.0E+10).
We conclude from our analyses that clinical improvements correspond to a progressive, if slow, enrichment of the faecal microbiota, even if the timing of the sampling may have obscured the evidence of this phenomenon. The corollary of this is that the microbiota remains reduced in number if clinical conditions do not improve for the patient.
The small sample size is a clear limitation of this study. However, for the first time, to our knowledge, we have shown that ecological alterations in gut microbiota in critically ill patients can be detected using molecular tools and that these data have clinical significance. Our findings underline: (a) the impact of even a single week of critical illness on the maintenance of the gut ecosystem and its negative influence on clinical outcome, and (b) the specific risk of Enterococcus outgrowth following both gut microbiota changes and clindamycin treatments. Moreover, these alterations, coupled with the effects of increased enterococci levels, seem to be associated with organ failure, the duration of intensive care, as well as the negative outcomes for the four patients in our current study cohort who died in the ICU.
| REFERENCES |
|---|
|
|
|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| INT J SYST EVOL MICROBIOL | J MED MICROBIOL | MICROBIOLOGY | J GEN VIROL | ALL SGM JOURNALS |