Microbes in the neonatal intensive care unit resemble those found in the gut of premature infants
© Brooks et al.; licensee BioMed Central Ltd. 2014
Received: 27 August 2013
Accepted: 28 November 2013
Published: 28 January 2014
The source inoculum of gastrointestinal tract (GIT) microbes is largely influenced by delivery mode in full-term infants, but these influences may be decoupled in very low birth weight (VLBW, <1,500 g) neonates via conventional broad-spectrum antibiotic treatment. We hypothesize the built environment (BE), specifically room surfaces frequently touched by humans, is a predominant source of colonizing microbes in the gut of premature VLBW infants. Here, we present the first matched fecal-BE time series analysis of two preterm VLBW neonates housed in a neonatal intensive care unit (NICU) over the first month of life.
Fresh fecal samples were collected every 3 days and metagenomes sequenced on an Illumina HiSeq2000 device. For each fecal sample, approximately 33 swabs were collected from each NICU room from 6 specified areas: sink, feeding and intubation tubing, hands of healthcare providers and parents, general surfaces, and nurse station electronics (keyboard, mouse, and cell phone). Swabs were processed using a recently developed ‘expectation maximization iterative reconstruction of genes from the environment’ (EMIRGE) amplicon pipeline in which full-length 16S rRNA amplicons were sheared and sequenced using an Illumina platform, and short reads reassembled into full-length genes. Over 24,000 full-length 16S rRNA sequences were produced, generating an average of approximately 12,000 operational taxonomic units (OTUs) (clustered at 97% nucleotide identity) per room-infant pair. Dominant gut taxa, including Staphylococcus epidermidis, Klebsiella pneumoniae, Bacteroides fragilis, and Escherichia coli, were widely distributed throughout the room environment with many gut colonizers detected in more than half of samples. Reconstructed genomes from infant gut colonizers revealed a suite of genes that confer resistance to antibiotics (for example, tetracycline, fluoroquinolone, and aminoglycoside) and sterilizing agents, which likely offer a competitive advantage in the NICU environment.
We have developed a high-throughput culture-independent approach that integrates room surveys based on full-length 16S rRNA gene sequences with metagenomic analysis of fecal samples collected from infants in the room. The approach enabled identification of discrete ICU reservoirs of microbes that also colonized the infant gut and provided evidence for the presence of certain organisms in the room prior to their detection in the gut.
From birth to death, humans spend approximately 90% of their time indoors . This realization, coupled with advancements in DNA sequencing technologies, has spawned a new interest in studying buildings as ecosystems. Pioneering efforts have revealed a built environment (BE), a term used here to collectively describe both the biotic and abiotic features of a building structure, that is far more complex than originally imagined [2, 3]. Diverse microbial communities have been uncovered in a variety of BEs  and surprisingly, from sites engineered to be sterile or near sterile, such as NASA clean rooms [5, 6] and high-risk hospital wards [7–10]. Additionally, recent studies characterizing different building types have revealed general trends suggesting a room’s function or architecture dictates the BE’s microbiome [8, 11]. Intrabuilding experiments in hospitals have corroborated this notion, showing general use areas, such as waiting rooms and lobbies, have a markedly different microbial community compared to more restrictive hospital zones such as intensive care units . The exchange between the BE microbiome and the human microbiome communities remains unclear; however, the observation that human pathogens are enriched for in hospital settings is of obvious concern . Here, we aimed to characterize the interaction between the BE’s microbiome and the human microbiome through study of very low birth weight (VLBW, <1,500 g) infants housed in a neonatal intensive care unit (NICU) as our model system.
Infants housed in a NICU are well suited to studies that aim to characterize interactions between the BE and occupants. In utero, infants are canonically thought to exist in a sterile or near-sterile environment . Acquisition of the microbiome starts at birth and is strongly influenced by mode of delivery . Patterns of colonization in full-term infants tend to follow a well documented trajectory affected by diet, host genotype, and a limited set of other variables, with the infant gut converging on an adult-like state around 2.5 years of life [14, 15]. In VLBW infants, early gut succession is characterized by extremely limited diversity, chaotic flux in community composition, and an abundance of opportunistic pathogens [16–19]. It is possible that a high rate of caesarean deliveries and the routine use of broad-spectrum antibiotics during the first week of life serve to decouple VLBW infants from source inoculum introduced during the birthing process. These influences likely render premature infant microbiomes especially susceptible to environmental influences.
There is strong evidence suggesting that the ICU serves as a reservoir of clinically relevant pathogens. ‘Outbreaks’ of disease in ICUs are relatively common, and a recent study estimated at least 38% of all ICU outbreaks could be attributed to microbial sources within the ICU environment, such as equipment, or personnel . In addition, upward of 63% of extremely preterm infants develop life-threatening infections . Epidemiologic investigations indicate environmental sources of infective agents in air , infant incubators [23, 24], sink drains , soap dispensers , thermometers , and baby toys . Clearly there is a growing need for comprehensive ecological surveys of the hospital BE to better understand the overall process of microbe migration and establishment on and in the body of occupants. Here, we performed the first matched time series characterization of the NICU and infant gut. Our analysis used metagenomic sequencing of microbial community DNA extracted from fecal samples to evaluate the metabolic potential of gut colonizing microorganisms and a recently developed ‘expectation maximization iterative reconstruction of genes from the environment’ (EMIRGE) amplicon protocol to profile the microbial community composition of BE samples collected from six environment types . Our protocol was aimed at addressing the hypothesis that the BE, specifically room surfaces frequently touched by humans, is a predominant source of colonizing microbes in the GI tract of premature infants.
Health profile of premature infant cohort
26 3/7 weeks
28 2/7 weeks
Day of life (DOL) 1 to 7 antibiotics
DOL 14 to 16, vancomycin, cefotaxime
DOL 3, maternal milk
DOL 8, artificial formula
Survive to discharge
All samples were obtained from a private-style NICU at Magee-Womens Hospital of the University of Pittsburgh Medical Center. Room samples were collected concurrently with fecal samples and spanned four timepoints on days of collection (9:00, 12:00, 13:00, and 16:00). Most frequently touched surfaces were determined by visual observation and health care provider interviews in the weeks leading up to sample collection. Microbial cells were removed from surfaces using foam tipped swabs (BBL CultureSwab EZ Collection and Transport System, Franklin Lakes, NJ, USA) and a sampling buffer of 0.15 M NaCl and 0.1% Tween20. Six frequently touched areas were processed per infant room: sink, feeding and intubation tubing, hands of healthcare providers and parents, general surfaces, access knobs on the incubator, and nurse station electronics (keyboard, mouse, and cell phone). All samples were placed in a sterile transport tube and stored within 30 minutes at -80°C until further processing.
DNA extraction and PCR amplification
Frozen fecal samples were thawed on ice and 0.25 g of thawed sample added to tubes with prewarmed (65°C) lysis solution from the PowerSoil DNA Isolation Kit (MoBio Laboratories, Carlsbad, CA, USA). The incubation was conducted for 5 minutes and the manufacturer’s protocol followed thereafter. Swab heads followed the same procedure, except heads were cut with sterilized scissors into the extraction tube before starting the protocol.
DNA extracted from swabs was pooled such that the four timepoints sampled in 1 day per environment were consolidated into one sample. Pooled DNA was used as template for amplification of the full-length 16S rRNA gene with 27 F (5’-AGAGTTTGATCCTGGCTCAG-3’) and 1492R (5’-GGTTACCTTGTTACGACTT-3’) primers . To limit PCR bias, gradient PCR was performed with 5 units/μL of TaKaRa Ex Taq™ (Takara Bio Inc., Otsu, Japan) across 7 different annealing temperatures with the following reaction: 1 minute at 94°C; 35 cycles of 1 minute at 94°C, 30 s at 48°C to 58°C (7°C temperature gradient) and 1 minute at 72°C; and a final extension for 7 minutes at 72°C. Amplicons were combined across gradients and cleaned with the QIAquick PCR Purification Kit (Qiagen, Hilden, Germany) as directed by the manufacturer. Cleaned amplicons were quantified via Qubit (Life Technologies, Carlsbad, CA, USA) and input into an Illumina library preparation pipeline.
Sequencing preparation and sequencing
Illumina library construction followed standard protocols at the University of California Davis DNA Technologies Core Facility (http://dnatech.genomecenter.ucdavis.edu) as previously described . Briefly, amplicons were fragmented to an average size of 225 bp using the Bioruptor NGS (Diagenode, Seraing, Belgium), and sheared fragments were used in a robotic library preparation protocol using the Appollo 324 robot (Integenx, Pleasanton, CA, USA) following the manufacturer’s instructions. Each sample was tagged with unique barcodes consisting of six nucleotides internal to the adapter read as a separate indexing read, and ligated to each fragment. There were 12 cycles of PCR enriched for adapter-ligated fragments before library quantification and validation. Fecal samples underwent the same preparation with two exceptions: (1) genomic DNA was used and (2) DNA was fragmented to 550 bp. Libraries were added, in equimolar amounts, to the Illumina HiSeq 2000 platform. Paired-end sequences were obtained with 100 cycles and the data processed with Casava version 1.8.2. Raw read data has been deposited in the NCBI Short Read Archive (accession number SRP033353).
EMIRGE assembly of full-length 16S rRNA gene amplicons
EMIRGE is an iterative template-guided assembler that relies on a database of 16S rRNA gene sequences to probabilistically generate full-length 16S rRNA gene sequences and provide the relative abundance of these sequences in the assayed consortia . For the reference database, we used version 108 of the SILVA SSU database, filtered to exclude sequences <1,200 bp and <1,900 bp . To remove closely related sequences, we clustered the database at 97% identity with USEARCH . A total of 1 million paired-end reads from each barcoded library were sampled randomly without replacement to accommodate computational restrictions associated with use of the full dataset. Reads from the subsample from each library were stringently trimmed using Sickle  for quality scores <30 and length <60 bp. Trimmed reads were input into an amplicon-optimized version of EMIRGE  for assembly using default parameters. A total of 80 iterations were performed for each subsample. EMIRGE-reconstructed sequences without Ns and with an estimated abundance of 0.01% or greater were kept for analysis. Putative chimeras were removed by using the intersection between two chimera detection programs, DECIPHER  and UCHIME v6.0  searched against the 2011 Greengenes database . Finally, reconstructed sequences from a spike-in control experiment (data not shown) were removed for downstream analysis. Sequences used in the analysis are publicly available as a project attachment at http://ggkbase.berkeley.edu/NICU-Micro/.
Metagenomic EMIRGE assembly of 16S rRNA gene
Metagenomic sequencing of 16 fecal samples on 1 lane of an Illumina HiSeq 2000 produced approximately 350 Mbp of 101 bp paired-end reads. Trimmed reads were input into EMIRGE and default parameters run for 80 iterations using the aforementioned database. After the final iteration, 153,980 reads, spanning all samples, were used in reconstructing fecal 16S rRNA sequences. Downstream filtering and analysis of reconstructed 16S rRNA gene sequences from fecal samples followed that of the room samples.
Community analysis of room and fecal samples
For community analysis, EMIRGE-reconstructed sequences were input into the standard QIIME 1.5.0 workflow . For presence/absence analyses, representative operational taxonomic units (OTUs) were clustered at the <97% identity level using USEARCH  and an OTU table was constructed using QIIME’s pick_otus_through_otu_table.py script. An adjusted OTU table that incorporated EMIRGE generated abundances was constructed using an in-house script  and is publicly available as a project attachment at http://ggkbase.berkeley.edu/NICU-Micro/. OTUs were aligned to the Greengenes  reference alignment (gg_97_otus_4feb2011.fasta) using the PyNAST aligner  and a phylogenetic tree built using FastTree v.2.1.3  with default parameters. Beta diversity was calculated from similar trees using Fast UniFrac scores and visualized with principle coordinates analysis (PCoA) . Taxonomy was assigned to each OTU at the genera and/or species level using the Ribosomal Database Project (RDP) classifier  at a confidence interval of 0.8 and trained with the same Greengenes database. OTUs were visualized across room-infant pairs in a spring-weighted, edge-embedded network plot by using QIIME’s make_otu_network.py script  with the modified OTU table as input.
Metagenomic assembly and gene prediction
Assemblies were constructed using idba_ud  and an iterative implementation of Velvet [45, 46]. For idba_ud assemblies, trimmed reads were assembled using default parameters. For the Velvet assemblies, sequence coverage bins representing major genomes in the dataset were identified by first running the program with permissive parameters in which the k-mer size covered the whole range of observed coverages. We summed the k-mer coverages for all contigs generated by this assembly to define the coverage bins (each of which contains one or more genomes). This provided bin-specific expected coverage, k-mer size, coverage cutoff, and coverage collection threshold parameters for the iterative assembly. After each iteration targeting a specific bin, the bin-specific reads were removed from the dataset.
Time-series-coverage-based emergent self-organizing maps (ESOMs) were used to bin scaffolds generated by metagenomic assembly . Genes were predicted and translated into protein sequences using Prodigal . Functional annotation was added with an in-house pipeline . Genome completeness was determined based on the number of single-copy genes and other conserved genes [49, 50] identified in each bin. The relative abundance of each organism in each sample was calculated by mapping reads to unique regions on the assembled genomes. Metagenomic assemblies along with their annotations are publicly available at http://ggkbase.berkeley.edu/NICU-Micro/.
Enterococcus faecalis concatenated ribosomal protein phylogeny
For phylogenetic resolution beyond the 16S rRNA gene, 32 highly conserved, single copy ribosomal proteins were used from infant 1 and 2’s assemblies (RpL10, 13, 14, 16, 17, 18, 19, 2, 20, 21, 22, 24, 27, 29, 3, 30, 4, 5, and RpS10, 11, 12, 13, 15, 16, 17, 18, 19, 20, 5, 6, 7, 8). The same genes from recently sequenced E. faecalis genomes, in addition to genes from more distantly related taxa, were obtained from the JGI IMG database. Together, each gene set was aligned using MUSCLE 3.8.31 [51, 52] and manually curated to remove ambiguously aligned regions and end gaps . The curated alignments were concatenated to form a 32-gene, 39-taxa, 4,101-position alignment. A maximum likelihood phylogeny for the concatenated alignment was conducted using PhyML under the LG + α + γ model of evolution with 100 bootstrap replicates.
Stability of NICU room samples over time and space
Sample collection summary and summary of the number of 16S rRNA genes assembled
No. of samples
No. of EMIRGE sequences
No. of OTUs
No. of unique OTUs
Alpha diversity indexes from neonatal intensive care unit (NICU) room and fecal samples
Time-series characterization of fecal samples
More than 94% of the reads from infant 1’s samples mapped to scaffolds generated by the idba_ud assembly. Consequently, this assembly was accepted for further analysis. In comparison, the initial idba_ud assembly of metagenomic data from infant 2 was highly fragmented, and less than 40% of reads could be mapped to the assembled scaffolds. Subsequent reassembly of metagenomic data from infant 2’s samples using the iterative Velvet-based assembly approach  generated a significantly better result. As <90% of reads could be mapped to the scaffolds generated by the Velvet assembly, this assembly was chosen for further analysis.
Staphylococcus epidermidis plasmid
Actinomyces neuii strain 1
Actinomyces neuii strain 2
Enterococcus faecalis phage
Enterococcus faecalis plasmid
Finegoldia plasmid 1
Finegoldia plasmid 2
Staphylococcus epidermidis mobile
Staphylococcus epidermidis phage 2
Staphylococcus epidermidis strain
Staphylococcus phage 1
Staphylococcus phage 2
Staphylococcus phage 3 and plasmid
Infant 1 and infant 2’s gastrointestinal tract (GIT) microbial communities are distinctly different. Infant 1’s colonization pattern echoes the canonical observation in infant GIT succession that facultative anaerobes dominate early phase colonization whereas late stage colonizers are primarily obligate anaerobes . This shift is observed on day of life 12 in infant 1, but is not observed in infant 2, in whom facultative anaerobes were observed throughout the study period. The metagenomic EMIRGE analyses corroborated the binning-based compositional analyses in that no sequences for new taxa were assembled for scaffolds included in the ESOM. Some 16S rRNA genes were identified in the metagenomic assemblies and match EMIRGE generated sequences with approximately 100% identity. The E. faecalis sequence from infant 1 was not identified by EMIRGE due to low abundance, but was extracted from the assembly using RNAmmer for the phylogenetic analysis .
Highly connected BE microbes
The NICU as a reservoir for gut colonists
Shared gut colonizers
To further explore similarity of shared strains, reads from infant 1 were mapped to infant 2’s assembled contigs. Infant 1’s reads covered 95% of the length of infant 2’s assembly at an average of 4.66X coverage. Read mapping revealed two distinct SNP profiles for infant 1’s reads, a major strain divergent from infant 2’s assembly and a minor strain identical to the strain in infant 2. In all, 77% of the length of infant 2’s E. faecalis assembly is covered by infant 1’s reads mapped as mate pairs with no mismatches. This suggests that infant 1’s E. faecalis minor strain is the same strain dominating infant 2’s gut. Pheromone-responsive plasmids were found in both infants. The plasmid from infant 2 occurs in low abundance in infant 1 (as expected based on the low representation of E. faecalis in infant 1), but with high sequence identity.
Genes relevant to adaptation to the NICU environment
Analysis of reconstructed genomes for gut microorganisms can lend clues as to how organisms detected in the GIT and room environment are able to persist in the NICU, which is subjected to regular cleaning/sterilization. Numerous antibiotic resistance genes were found in genomes of microorganisms in fecal samples of both infants. A large portion of these were efflux pumps, with representatives from all four families of multidrug transporters: major facilitator superfamily (MFS), small multidrug resistant (SMR), resistance-nodulation-cell division (RND), and multidrug and toxic compound extrusion (MATE) proteins . Particularly interesting are genes encoding the QacA/B MFS, SugE SMR, and MexA/B RND proteins, which are a growing concern in hospitals due to coselection through the practice of combining two or more types of antibiotic treatments . Resistance to multiple types of antibiotics can arise from a single resistance mechanism such as efflux pumping . In addition to antibiotics, these pumps can expel quaternary ammonium compounds (QACs), the active biocide in the detergent used to clean hospital surfaces during the study. Other notable observations were the presence of biofilm forming genes in most colonizers, which can be induced by exposure to aminoglycosides , a suite of genes that confer resistance to starvation, and the presence of antibiotic resistance genes encoded on several phage and plasmid genomes, as well as microbial genomes.
Increasing throughput, decreasing cost, and rapid development of informatics and sequencing pipelines has reshaped the field of microbial ecology, allowing researchers to survey a breadth of new environments [34, 61–63]. Recently, the first ICU survey to utilize next generation sequencing technology was published  and showed a surprising amount of bacterial diversity for an environment under constant attack via aggressive sanitation and antibiotic treatment efforts. The consortia were generally diverse, but some consortia contained a high representation of members of the family Enterobacteriaceae, typically considered to be gut microbes. Shortly after this publication, a study characterizing a snapshot of surfaces and sinks in two NICU rooms corroborated high proportions of fecal coliform bacteria on surface samples . Certainly the NICU has the capacity to retain enteric microbes, but their propensity to migrate to the gut remains unclear.
Next-generation sequencing surveys in the ICU have reported high levels of community diversity. Poza et al. found 1,145 distinct OTUs in an ICU in Spain  and subsequent studies reported 1,621 and 3,925 OTUs in a NICU in the US and in an Austrian ICU, respectively [9, 10]. While comparing these studies is difficult due to differences in sample size and protocols, we can begin to appreciate the need to better understand why so many types of bacteria can be found in a regularly cleaned environment. Our study, the first time series survey of an ICU using next-generation sequencing technologies, unveiled over 20,000 OTUs across 2 NICU rooms occupied by different infants with partial time overlap. Our study is distinct from prior NICU surveys in that it used amplicon-EMIRGE, a 16S rRNA gene assembly software which can be more sensitive in OTU detection  and provide increased confidence when making lower taxonomic level classifications . The increase in OTUs from study to study might be attributed to increases in sequencing read lengths and, in this study, increased information from reassembled, full-length genes, but the biological relevance of this increase is unclear. Notably, of the over 20,000 OTUs characterized here, only 984 were found in 2 or more samples. Further surveys are needed, integrating time-series sampling and samples from multiple surface types from different hospitals, to better characterize the expected number of OTUs in an ICU and the implications of this number for ICU occupants.
The increased sensitivity provided by EMIRGE was helpful when evaluating temporal patterns, especially pertaining to source-sink characterization. Similarly, our source-sink analyses benefited from the increased number of samples and timepoints relative to prior studies [8–10], which did not attempt to identify source-sink relationships. The SourceTracker results suggest the most probable room reservoir for gut colonists is tubing followed by surfaces, incubators, and hands (Figure 6). The tubing area sampled, the hub of the silastic nasogastric feeding tube, is the closest in proximity to the infant and, since SourceTracker is not bidirectional, it is difficult to tease out the directionality in this exchange . Incubators from both infants also appear to mirror successional patterns in the infant’s GIT, but without finer scale temporal sampling it is difficult to determine the true source and sink. The observation that hands tend to show a variable amount of potential fecal colonist is likely due to the variability in sampling and hand hygiene, as hand samples were taken both before and after infants received care from healthcare providers. A good example of this is infant 1’s DOL 27 hand sample in which the large spike in Escherichia likely came from a swab collected directly after contact with the infant (Figure 5).
Given the large inventory of sequences and the time-series dataset, it was possible to identify likely reservoirs of microorganisms in the room environment, prior to their appearance in the GIT (for example, the asterisked OTUs in Figure 4). Many of these sequences had perfect or near perfect identity between room and GIT 16S rRNA genes. Two notable examples include the Klebsiella pneumoniae in infant 1 and Finegoldia magna in infant 2, whose fecal to room sequence best hits averaged 99.4% and 99.6% identity respectively. Infant 1’s K. pneumoniae is first detectable in the gut on DOL 9, but NICU samples first detect the organism on electronic and sink samples starting at DOL 3, our earliest sampling point. Interestingly, the K. pneumoniae is outcompeted in the gut, yet is reintroduced on two separate occasions. This could be a byproduct of our detection limits, but the relatively high abundance of Klebsiella in many NICU samples and its availability at all timepoints, suggests the opportunity for reinoculation from multiple room reservoirs. The F. magna in infant 2’s samples exhibit similar patterns in that it is initially a high-ranking taxa that is out competed by other Firmicutes, but is reintroduced later in the time series.
If the environment is a reservoir for gut colonizing microbes in our cohort, then it is likely infants housed in close proximity will share the same strain. The 16S rRNA gene survey shows the availability of reservoirs of colonizing populations (likely with multiple strain variants) in the infant’s immediate environment. However, it cannot discriminate at the strain level, so the mere existence of a phylotype in the room prior to gut colonization is not a direct measure of BE to infant transfer. The current work resolves this, by using extensive genome sequence comparison of E. faecalis from the gut of two infants housed in the same ward to establish that environment to room occupant transfer occurs in the NICU. The mode of acquisition of infant 2’s abundant strain by infant 1 is unclear, but nosocomial infection by enterococci is not uncommon.
Enterococci are particularly difficult to classify due the plasticity of their genomes. Upwards of 25% of E. faecalis genomes may be comprised of mobile or acquired elements . Recent experiments attribute this genome flexibility partially to the ability to produce transconjugant hybrid strains in which several 100 kb fragments can be transferred between donor and recipient strain . Transfer of these genome fragments is dependent on pheromone-responsive plasmids, which were found in all strains studied here. The ability to form hybrids not only confounds the ability to confirm identical strains, unless the entire genome has been recovered, it also provides a competitive advantage in the hospital BE where enterococci have been problematic for decades [65, 67]. Enterococci are notoriously hardy and are able to persist on medical equipment and hospital surfaces for long periods of time [65, 68]. They are able to withstand chlorine, heat, some alcohol treatments, and possibly most concerning, several types of antibiotics . Their genome plasticity and ability to easily acquire new genes from other strains make them particularly well suited to thrive in the hospital environment.
Gut colonists must withstand selective pressures both inside and out of the gut. Two obvious forms of selection in the NICU come from hospital cleaning and the broad use of antibiotics. All rooms were cleaned daily using wet solutions containing QACs and all infants were administered multiple types of antibiotics. Incorrect administration of biocides, through misuse or unintended mixing with existing fluids (that is, water from sink samples or removing sanitizing agents via water rinsing), could enrich for resistance genes . Even if used to factory standards, if surface-dried cells or biofilms remain, biocide activity could be ineffective and contribute to cross resistance to biocides and antibiotics . Biofilm forming communities can be upwards of 1,000 times more resistant to QACs than their planktonic forms  and biofilm formation can be triggered by the types of antibiotics administered in this study . This may be a contributing factor as to why a recent study found enteric microbial communities to be relatively unaltered before and after routine NICU surface cleaning . Certain types of biofilms in many Enterobacteriaceae, including those studied here, contain amyloid fibers, called curli. Curli have been implicated in adhesion to abiotic surfaces, such as polystyrene, Teflon, and stainless steel, and contribute to adhesion to host epithelial cells and invasion by Escherichia coli in the gut . This type of dual-purpose adaptation may allow enteric organisms to persist on NICU surfaces until transmission to a more favorable environment such as the gut. Efflux pumps are another multipurpose adaptation conferring competitive advantages inside and out of the gut. Numerous pumps from every major class of efflux pump were identified here and, collectively, can function to pump out QACs and administered antibiotics. Previous studies have positively correlated high QAC minimum inhibitory concentrations (MICs) with increased antibiotic resistance markers in enteric microbes , indicating biocide efflux may be an important function for microbes in the ICU. Efflux and biofilm formation are two of many possible explanations as to how colonizers combat both biocides administered during NICU cleaning and host-administered antibiotics.
Through a time series analysis using full-length rRNA gene sequences, we have established that organisms that appear in the GI tract in the early phase of colonization have reservoirs in the room environment. The findings point to a scenario in which gut microbes are introduced from room sources, thrive in the gut, and are disseminated to the immediate environment, creating a cycle of room to infant colonization. The research also highlights the value of extensive genome comparisons to link colonists from different individuals, an approach that in the future may also target populations sampled directly from room reservoirs.
Day of life
Emergent self-organizing map
Intensive care unit
Multidrug and toxic compound extrusion
Major facilitator superfamily
Minimum inhibitory concentration
Neonatal intensive care unit
Operational taxonomic unit
Principle coordinates analysis
Quaternary ammonium compound
Small multidrug resistant
Very low birth weight.
Funding was provided through the Alfred P Sloan Foundation and the National Science Foundation’s Graduate Research Fellowship Program.
- Klepeis NE, Nelson WC, Ott WR, Robinson JP, Tsang AM, Switzer P, Behar JV, Hern SC, Engelmann WH: The National Human Activity Pattern Survey (NHAPS): a resource for assessing exposure to environmental pollutants. J Expo Anal Environ Epidemiol. 2001, 11: 231-252. 10.1038/sj.jea.7500165.View ArticlePubMedGoogle Scholar
- Tringe SG, Zhang T, Liu X, Yu Y, Lee WH, Yap J, Yao F, Suan ST, Ing SK, Haynes M, Rohwer F, Wei CL, Tan P, Bristow J, Rubin EM, Ruan Y: The airborne metagenome in an indoor urban environment. PLoS One. 2008, 3: e1862-10.1371/journal.pone.0001862.View ArticlePubMedPubMed CentralGoogle Scholar
- Rintala H, Pitkäranta M, Toivola M, Paulin L, Nevalainen A: Diversity and seasonal dynamics of bacterial community in indoor environment. BMC Microbiol. 2008, 8: 56-10.1186/1471-2180-8-56.View ArticlePubMedPubMed CentralGoogle Scholar
- Kelley ST, Gilbert JA: Studying the microbiology of the indoor environment. Genome Biol. 2013, 14: 202-10.1186/gb-2013-14-2-202.View ArticlePubMedPubMed CentralGoogle Scholar
- La Duc MT, Dekas A, Osman S, Moissl C, Newcombe D, Venkateswaran K: Isolation and characterization of bacteria capable of tolerating the extreme conditions of clean room environments. Appl Environ Microbiol. 2007, 73: 2600-2611. 10.1128/AEM.03007-06.View ArticlePubMedPubMed CentralGoogle Scholar
- La Duc MT, Vaishampayan P, Nilsson HR, Torok T, Venkateswaran K: Pyrosequencing-derived bacterial, archaeal, and fungal diversity of spacecraft hardware destined for Mars. Appl Environ Microbiol. 2012, 78: 5912-5922. 10.1128/AEM.01435-12.View ArticlePubMedPubMed CentralGoogle Scholar
- Perkins SD, Mayfield J, Fraser V, Angenent LT: Potentially pathogenic bacteria in shower water and air of a stem cell transplant unit. Appl Environ Microbiol. 2009, 75: 5363-5372. 10.1128/AEM.00658-09.View ArticlePubMedPubMed CentralGoogle Scholar
- Poza M, Gayoso C, Gómez MJ, Rumbo-Feal S, Tomás M, Aranda J, Fernández A, Bou G: Exploring bacterial diversity in hospital environments by GS-FLX titanium pyrosequencing. PLoS One. 2012, 7: e44105-10.1371/journal.pone.0044105.View ArticlePubMedPubMed CentralGoogle Scholar
- Oberauner L, Zachow C, Lackner S, Högenauer C, Smolle K-H, Berg G: The ignored diversity: complex bacterial communities in intensive care units revealed by 16S pyrosequencing. Sci Rep. 2013, 3: 1-12.View ArticleGoogle Scholar
- Hewitt KM, Mannino FL, Gonzalez A, Chase JH, Caporaso JG, Knight R, Kelley ST: Bacterial diversity in two neonatal intensive care units (NICUs). PLoS One. 2013, 8: e54703-10.1371/journal.pone.0054703.View ArticlePubMedPubMed CentralGoogle Scholar
- Kembel SW, Jones E, Kline J, Northcutt D, Stenson J, Womack AM, Bohannan BJ, Brown GZ, Green JL: Architectural design influences the diversity and structure of the built environment microbiome. ISME J. 2012, 6: 1469-1479. 10.1038/ismej.2011.211.View ArticlePubMedPubMed CentralGoogle Scholar
- Penders J, Thijs C, Vink C, Stelma FF, Snijders B, Kummeling I, Van Den Brandt PA, Stobberingh EE: Factors influencing the composition of the intestinal microbiota in early infancy. Pediatrics. 2006, 118: 511-521. 10.1542/peds.2005-2824.View ArticlePubMedGoogle Scholar
- Dominguez-Bello MG, Costello EK, Contreras M, Magris M, Hidalgo G, Fierer N, Knight R: Delivery mode shapes the acquisition and structure of the initial microbiota across multiple body habitats in newborns. Proc Natl Acad Sci U S A. 2010, 107: 11971-11975. 10.1073/pnas.1002601107.View ArticlePubMedPubMed CentralGoogle Scholar
- Palmer C, Bik EM, DiGiulio DB, Relman DA, Brown PO: Development of the human infant intestinal microbiota. PLoS Biol. 2007, 5: e177-10.1371/journal.pbio.0050177.View ArticlePubMedPubMed CentralGoogle Scholar
- Trosvik P, Stenseth NC, Rudi K: Convergent temporal dynamics of the human infant gut microbiota. ISME J. 2010, 4: 151-158. 10.1038/ismej.2009.96.View ArticlePubMedGoogle Scholar
- Morowitz MJ, Denef VJ, Costello EK, Thomas BC, Poroyko V, Relman DA, Banfield JF: Strain-resolved community genomic analysis of gut microbial colonization in a premature infant. Proc Natl Acad Sci U S A. 2011, 108: 1128-1133. 10.1073/pnas.1010992108.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang Y, Hoenig JD, Malin KJ, Qamar S, Petrof EO, Sun J, Antonopoulos DA, Chang EB, Claud EC: 16S rRNA gene-based analysis of fecal microbiota from preterm infants with and without necrotizing enterocolitis. ISME J. 2009, 3: 944-954. 10.1038/ismej.2009.37.View ArticlePubMedPubMed CentralGoogle Scholar
- Morowitz MJ, Poroyko V, Caplan M, Alverdy J, Liu DC: Redefining the role of intestinal microbes in the pathogenesis of necrotizing enterocolitis. Pediatrics. 2010, 125: 777-785. 10.1542/peds.2009-3149.View ArticlePubMedGoogle Scholar
- Mshvildadze M, Neu J, Shuster J, Theriaque D, Li N, Mai V: Intestinal microbial ecology in premature infants assessed with non-culture-based techniques. J Pediatr. 2010, 156: 20-25. 10.1016/j.jpeds.2009.06.063.View ArticlePubMedPubMed CentralGoogle Scholar
- Gastmeier P, Loui A, Stamm-Balderjahn S, Hansen S, Zuschneid I, Sohr D, Behnke M, Obladen M, Vonberg R-P, Rüden H: Outbreaks in neonatal intensive care units - they are not like others. Am J Infect Control. 2007, 35: 172-176. 10.1016/j.ajic.2006.07.007.View ArticlePubMedGoogle Scholar
- Stoll BJ, Hansen NI, Bell EF, Shankaran S, Laptook AR, Walsh MC, Hale EC, Newman NS, Schibler K, Carlo WA, Kennedy KA, Poindexter BB, Finer NN, Ehrenkranz RA, Duara S, Sánchez PJ, O’Shea TM, Goldberg RN, Van Meurs KP, Faix RG, Phelps DL, Frantz ID, Watterberg KL, Saha S, Das A, Higgins RD: Neonatal outcomes of extremely preterm infants from the NICHD Neonatal Research Network. Pediatrics. 2010, 126: 443-456. 10.1542/peds.2009-2959.View ArticlePubMedPubMed CentralGoogle Scholar
- Adler A, Gottesman G, Dolfin T, Arnon S, Regev R, Bauer S, Litmanovitz I: Bacillus species sepsis in the neonatal intensive care unit. J Infect. 2005, 51: 390-395. 10.1016/j.jinf.2004.12.006.View ArticlePubMedGoogle Scholar
- Singh N, Campbell J, Short B: Control of vancomycin-resistant enterococci in the neonatal intensive care unit. Infect Control. 2005, 26: 646-649.Google Scholar
- Touati A, Achour W, Cherif A, Hmida HB, Afif FB, Jabnoun S, Khrouf N, Hassen AB: Outbreak of Acinetobacter baumannii in a neonatal intensive care unit: antimicrobial susceptibility and genotyping analysis. Ann Epidemiol. 2009, 19: 372-378. 10.1016/j.annepidem.2009.03.010.View ArticlePubMedGoogle Scholar
- Bonora M, Ligozzi M, De Fatima M: Vancomycin-resistant Enterococcus faecium isolates causing hospital outbreaks in northern Italy belong to the multilocus sequence typing C1 lineage. Microb Drug. 2004, 10: 114-123. 10.1089/1076629041310046.View ArticleGoogle Scholar
- Buffet-Bataillon S, Rabier V, Bétrémieux P, Beuchée A, Bauer M, Pladys P, Le Gall E, Cormier M, Jolivet-Gougeon A: Outbreak of Serratia marcescens in a neonatal intensive care unit: contaminated unmedicated liquid soap and risk factors. J Hosp Infect. 2009, 72: 17-22. 10.1016/j.jhin.2009.01.010.View ArticlePubMedGoogle Scholar
- Van den Berg RW, Claahsen HL, Niessen M, Muytjens HL, Liem K, Voss A: Enterobacter cloacae outbreak in the NICU related to disinfected thermometers. J Hosp Infect. 2000, 45: 29-34. 10.1053/jhin.1999.0657.View ArticlePubMedGoogle Scholar
- Naesens R, Jeurissen A, Vandeputte C, Cossey V, Schuermans A: Washing toys in a neonatal intensive care unit decreases bacterial load of potential pathogens. J Hosp Infect. 2009, 71: 197-198. 10.1016/j.jhin.2008.10.018.View ArticlePubMedGoogle Scholar
- Miller CS, Handley KM, Wrighton KC, Frischkorn KR, Thomas BC, Banfield JF: Short-read assembly of full-length 16S amplicons reveals bacterial diversity in subsurface sediments. PLoS One. 2013, 8: e56018-10.1371/journal.pone.0056018.View ArticlePubMedPubMed CentralGoogle Scholar
- Stackebrandt E, Goodfellow M: 16S/23S rRNA sequencing. Nucleic Acid Techniques in Bacterial Systematics. 1991, Chichester, UK: John Wiley & Son Ltd, 115-175.Google Scholar
- Miller CS, Baker BJ, Thomas BC, Singer SW, Banfield JF: EMIRGE: reconstruction of full-length ribosomal genes from microbial community short read sequencing data. Genome Biol. 2011, 12: R44-10.1186/gb-2011-12-5-r44.View ArticlePubMedPubMed CentralGoogle Scholar
- Pruesse E, Quast C, Knittel K, Fuchs BM, Ludwig W, Peplies J, Glöckner FO: SILVA: a comprehensive online resource for quality checked and aligned ribosomal RNA sequence data compatible with ARB. Nucleic Acids Res. 2007, 35: 7188-7196. 10.1093/nar/gkm864.View ArticlePubMedPubMed CentralGoogle Scholar
- Edgar RC: Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010, 26: 2460-2461. 10.1093/bioinformatics/btq461.View ArticlePubMedGoogle Scholar
- Najoshi G: Sickle - A windowed adaptive trimming tool for FASTQ files using quality.https://github.com/najoshi/sickle,
- Wright ES, Yilmaz LS, Noguera DR: DECIPHER, a search-based approach to chimera identification for 16S rRNA sequences. Appl Environ Microbiol. 2012, 78: 717-725. 10.1128/AEM.06516-11.View ArticlePubMedPubMed CentralGoogle Scholar
- Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R: UCHIME improves sensitivity and speed of chimera detection. Bioinformatics. 2011, 27: 2194-2200. 10.1093/bioinformatics/btr381.View ArticlePubMedPubMed CentralGoogle Scholar
- McDonald D, Price MN, Goodrich J, Nawrocki EP, DeSantis TZ, Probst A, Andersen GL, Knight R, Hugenholtz P: An improved Greengenes taxonomy with explicit ranks for ecological and evolutionary analyses of bacteria and archaea. ISME J. 2012, 6: 610-618. 10.1038/ismej.2011.139.View ArticlePubMedPubMed CentralGoogle Scholar
- Caporaso J, Kuczynski J, Stombaugh J: QIIME allows analysis of high-throughput community sequencing data. Nature. 2010, 7: 335-336.Google Scholar
- DeSantis TZ, Hugenholtz P, Larsen N, Rojas M, Brodie EL, Keller K, Huber T, Dalevi D, Hu P, Andersen GL: Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl Environ Microbiol. 2006, 72: 5069-5072. 10.1128/AEM.03006-05.View ArticlePubMedPubMed CentralGoogle Scholar
- Caporaso JG, Bittinger K, Bushman FD, DeSantis TZ, Andersen GL, Knight R: PyNAST: a flexible tool for aligning sequences to a template alignment. Bioinformatics. 2010, 26: 266-267. 10.1093/bioinformatics/btp636.View ArticlePubMedPubMed CentralGoogle Scholar
- Price MN, Dehal PS, Arkin AP: FastTree 2–approximately maximum-likelihood trees for large alignments. PLoS One. 2010, 5: e9490-10.1371/journal.pone.0009490.View ArticlePubMedPubMed CentralGoogle Scholar
- Hamady M, Lozupone C, Knight R: Fast UniFrac: facilitating high-throughput phylogenetic analyses of microbial communities including analysis of pyrosequencing and PhyloChip data. ISME J. 2010, 4: 17-27. 10.1038/ismej.2009.97.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang Q, Garrity GM, Tiedje JM, Cole JR: Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol. 2007, 73: 5261-5267. 10.1128/AEM.00062-07.View ArticlePubMedPubMed CentralGoogle Scholar
- Peng Y, Leung HCM, Yiu SM, Chin FYL: IDBA-UD: A de novo assembler for single-cell and metagenomic sequencing data with highly uneven depth. Bioinformatics. 2012, 28: 1420-1428. 10.1093/bioinformatics/bts174.View ArticlePubMedGoogle Scholar
- Zerbino DR, Birney E: Velvet: algorithms for de novo short read assembly using de Bruijn graphs. Genome Res. 2008, 18: 821-829. 10.1101/gr.074492.107.View ArticlePubMedPubMed CentralGoogle Scholar
- Sharon I, Morowitz MJ, Thomas BC, Costello EK, Relman DA, Banfield JF: Time series community genomics analysis reveals rapid shifts in bacterial species, strains, and phage during infant gut colonization. Genome Res. 2013, 23: 111-120. 10.1101/gr.142315.112.View ArticlePubMedPubMed CentralGoogle Scholar
- Ultsch A, Fabian M: ESOM-Maps: tools for clustering, visualization, and classification with Emergent SOM.http://www.cbs.dtu.dk/courses/27618.chemo/ESOM.pdf,
- Hyatt D, Chen G-L, Locascio PF, Land ML, Larimer FW, Hauser LJ: Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinformatics. 2010, 11: 119-10.1186/1471-2105-11-119.View ArticlePubMedPubMed CentralGoogle Scholar
- Sorek R, Zhu Y, Creevey CJ, Francino MP, Bork P, Rubin EM: Genome-wide experimental determination of barriers to horizontal gene transfer. Science. 2007, 318: 1449-1452. 10.1126/science.1147112.View ArticlePubMedGoogle Scholar
- Wu M, Eisen JA: A simple, fast, and accurate method of phylogenomic inference. Genome Biol. 2008, 9: R151-10.1186/gb-2008-9-10-r151.View ArticlePubMedPubMed CentralGoogle Scholar
- Edgar RC: MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC Bioinformatics. 2004, 5: 113-10.1186/1471-2105-5-113.View ArticlePubMedPubMed CentralGoogle Scholar
- Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32: 1792-1797. 10.1093/nar/gkh340.View ArticlePubMedPubMed CentralGoogle Scholar
- Hug LA, Castelle CJ, Wrighton KC, Thomas BC, Sharon I, Frischkorn KR, Williams KH, Tringe SG, Banfield JF: Community genomic analyses constrain the distribution of metabolic traits across the Chloroflexi phylum and indicate roles in sediment carbon cycling. Microbiome. 2013, 1: 22-10.1186/2049-2618-1-22.View ArticlePubMedPubMed CentralGoogle Scholar
- Sharon I, Morowitz MJ, Thomas BC, Costello EK, Relman DA, Banfield JF: Time series community genomics analysis reveals rapid shifts in bacterial species, strains, and phage during infant gut colonization. Genome Res. 2013, 23: 111-120. 10.1101/gr.142315.112. Supplemental Material. http://genome.cshlp.org/content/suppl/2012/10/15/gr.142315.112.DC2/Supplemental_Material.pdfView ArticlePubMedPubMed CentralGoogle Scholar
- Lagesen K, Hallin P, Rødland EA, Staerfeldt H-H, Rognes T, Ussery DW: RNAmmer: consistent and rapid annotation of ribosomal RNA genes. Nucleic Acids Res. 2007, 35: 3100-3108. 10.1093/nar/gkm160.View ArticlePubMedPubMed CentralGoogle Scholar
- Knights D, Kuczynski J, Charlson E: Bayesian community-wide culture-independent microbial source tracking. Nat Methods. 2011, 8: 6-10. 10.1038/nchembio.740.View ArticleGoogle Scholar
- Van Veen HW: Structural biology: Last of the multidrug transporters. Nature. 2010, 467: 926-927. 10.1038/467926a.View ArticlePubMedGoogle Scholar
- Fernández L, Hancock REW: Adaptive and mutational resistance: role of porins and efflux pumps in drug resistance. Clin Microbiol Rev. 2012, 25: 661-681. 10.1128/CMR.00043-12.View ArticlePubMedPubMed CentralGoogle Scholar
- Buffet-Bataillon S, Tattevin P, Bonnaure-Mallet M, Jolivet-Gougeon A: Emergence of resistance to antibacterial agents: the role of quaternary ammonium compounds - a critical review. Int J Antimicrob Agents. 2012, 39: 381-389. 10.1016/j.ijantimicag.2012.01.011.View ArticlePubMedGoogle Scholar
- Hoffman LR, D’Argenio DA, MacCoss MJ, Zhang Z, Jones RA, Miller SI: Aminoglycoside antibiotics induce bacterial biofilm formation. Nature. 2005, 436: 1171-1175. 10.1038/nature03912.View ArticlePubMedGoogle Scholar
- Mackelprang R, Waldrop MP, DeAngelis KM, David MM, Chavarria KL, Blazewicz SJ, Rubin EM, Jansson JK: Metagenomic analysis of a permafrost microbial community reveals a rapid response to thaw. Nature. 2011, 480: 368-371. 10.1038/nature10576.View ArticlePubMedGoogle Scholar
- Dick GJ, Anantharaman K, Baker BJ, Li M, Reed DC, Sheik CS: The microbiology of deep-sea hydrothermal vent plumes: ecological and biogeographic linkages to seafloor and water column habitats. Front Microbiol. 2013, 4: 124-View ArticlePubMedPubMed CentralGoogle Scholar
- Fierer N, Lauber CL, Zhou N, McDonald D, Costello EK, Knight R: Forensic identification using skin bacterial communities. Proc Natl Acad Sci U S A. 2010, 107: 6477-6481. 10.1073/pnas.1000162107.View ArticlePubMedPubMed CentralGoogle Scholar
- Ong SH, Kukkillaya VU, Wilm A, Lay C, Ho EXP, Low L, Hibberd ML, Nagarajan N: Species identification and profiling of complex microbial communities using shotgun illumina sequencing of 16S rRNA amplicon sequences. PLoS One. 2013, 8: e60811-10.1371/journal.pone.0060811.View ArticlePubMedPubMed CentralGoogle Scholar
- Arias CA, Murray BE: The rise of the Enterococcus: beyond vancomycin resistance. Nat Rev Microbiol. 2012, 10: 266-278. 10.1038/nrmicro2761.View ArticlePubMedPubMed CentralGoogle Scholar
- Manson JM, Hancock LE, Gilmore MS: Mechanism of chromosomal transfer of Enterococcus faecalis pathogenicity island, capsule, antimicrobial resistance, and other traits. Proc Natl Acad Sci U S A. 2010, 107: 12269-12274. 10.1073/pnas.1000139107.View ArticlePubMedPubMed CentralGoogle Scholar
- Murray BE: The Life and Times of the Enterococcus. Clin Microbiol Rev. 1990, 3: 46-65.PubMedPubMed CentralGoogle Scholar
- Bradley CR, Fraise AP: Heat and chemical resistance of enterococci. J Hosp Infect. 1996, 34: 191-196. 10.1016/S0195-6701(96)90065-1.View ArticlePubMedGoogle Scholar
- Condell O, Iversen C, Cooney S, Power KA, Walsh C, Burgess C, Fanning S: Efficacy of biocides used in the modern food industry to control salmonella enterica, and links between biocide tolerance and resistance to clinically relevant antimicrobial compounds. Appl Environ Microbiol. 2012, 78: 3087-3097. 10.1128/AEM.07534-11.View ArticlePubMedPubMed CentralGoogle Scholar
- Weiss-Muszkat M, Shakh D, Zhou Y, Pinto R, Belausov E, Chapman MR, Sela S: Biofilm formation by and multicellular behavior of Escherichia coli O55:H7, an atypical enteropathogenic strain. Appl Environ Microbiol. 2010, 76: 1545-1554. 10.1128/AEM.01395-09.View ArticlePubMedPubMed CentralGoogle Scholar
- Romanova NA, Gawande PV, Brovko LY, Griffiths MW: Rapid methods to assess sanitizing efficacy of benzalkonium chloride to Listeria monocytogenes biofilms. J Microbiol Methods. 2007, 71: 231-237. 10.1016/j.mimet.2007.09.002.View ArticlePubMedGoogle Scholar
- Bokulich NA, Mills DA, Underwood M: Surface microbes in the neonatal intensive care unit: changes with routine cleaning and over time. J Clin Microbiol. 2013, 51: 2617-2624. 10.1128/JCM.00898-13.View ArticlePubMedPubMed CentralGoogle Scholar
- Barnhart MM, Chapman MR: Curli biogenesis and function. Annu Rev Microbiol. 2006, 60: 131-147. 10.1146/annurev.micro.60.080805.142106.View ArticlePubMedPubMed CentralGoogle Scholar
- Buffet-Bataillon S, Branger B, Cormier M, Bonnaure-Mallet M, Jolivet-Gougeon A: Effect of higher minimum inhibitory concentrations of quaternary ammonium compounds in clinical E. coli isolates on antibiotic susceptibilities and clinical outcomes. J Hosp Infect. 2011, 79: 141-146. 10.1016/j.jhin.2011.06.008.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.