Microbes in the neonatal intensive care unit resemble those found in the gut of premature infants
Microbiomevolume 2, Article number: 1 (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.
Fecal samples were collected every third day, starting on the third day of life, for 1 month from two infants. Infants were enrolled in the study based on the criteria that they were <31 weeks’ gestation, <1,250 g at birth, and were housed in the same physical location within the NICU during the first month of life. A summary of health-related metadata including antibiotics exposure is provided in Table 1. Fecal samples were collected using a previously established perineal stimulation procedure and were stored at -80°C within 10 minutes . All samples were collected after signed guardian consent was obtained, as outlined in our protocol to the ethical research board of the University of Pittsburgh (IRB PRO11060238). This consent included sample collection permissions and consent to publish study findings.
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
After sample preparation, 57 and 36 samples amplified successfully and were subsequently analyzed for infant 1 and infant 2, respectively (Table 2). EMIRGE generated approximately 12,000 full-length 16S rRNA sequences and OTUs for each room-infant pair (clustered at the 97% nucleotide identity level). Broadly speaking, species richness decreased from electronics < sinks < surfaces < incubators < hands < tubes, a finding that was corroborated with several alpha diversity indexes (Table 3). Nearly 300 genera were detected in the NICU. To broadly visualize temporal stability of environments across time and space, the phylum level classifications are plotted in Figure 1. Actinobacteria, Firmicutes, and Proteobacteria dominate the sampled environments, with areas most exposed to human skin deposition having the most variation over time. At lower taxonomic levels, similar trends are observed. Based on the 20 most abundant families, frequently touched surfaces are distinct from infrequently touched surfaces (Figure 1). UniFrac distance-based community composition PCoA reveals four discernible ecosystem types (skin associated communities, sinks, tubes, and feces) and confirms clustering of samples prone to skin deposition via touching (Figure 2).
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.
The de novo assemblies reconstructed a majority of the genomes for 4 of the 5 and 8 of the 11 most abundant bacterial colonists from infant 1 and infant 2’s metagenomes, respectively. For infant 1, time-series organism abundance patterns in the sample sets analyzed via ESOM (Figure 3) defined five major genome bins for which between 37% and 99% of the single copy genes were identified, based on standard analyses of the single copy gene inventory (Table 4). For infant 2, time-series organism abundance patterns in the sample sets analyzed via ESOM (Figure 3) defined 11 major genome bins for which between 27 and 99% of the single copy genes were identified (Table 4).
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 distribution of shared OTUs across sampled sites was visualized through a spring-weighted edge-embedded network plot. To limit the noise from infrequently detected microorganism types, we restricted the plot to OTUs occurring in two or more samples from each infant (Figure 4). The spring weight is derived from EMIRGE generated abundances, and the distribution of OTUs in the plot is governed both by frequency of occurrence and abundance. In Figure 4, the circular white nodes (representing OTUs) found in many environment types (more edges) are pulled closer to the middle of the network whereas OTUs shared by only two samples (fewer edges) are positioned closer to the periphery of the network. The top 5% of most frequently occurring OTUs aggregate in a central cluster in the middle of the network. Similar to the PCoA plot, general clustering is observed based on environment type (that is, skin-associated sites cluster together, as do sink samples). When restricting the network for OTUs only found in fecal samples (Figure 4, enlargements), one can visualize the OTU distribution across the sampled NICU environments. Three highly connected OTUs are present in fecal samples, two of which are in the top 5% most frequently occurring OTUs in infant 1’s room samples. Several of the OTUs in infant 2’s fecal samples fall within the top ten most frequently occurring OTUs in the room environment. Interestingly, infant 2’s most abundant gut colonists, Staphylococcus sp. and E. faecalis, are the two most frequently occurring OTUs in the room environment.
The NICU as a reservoir for gut colonists
Figure 5 summarizes the gut colonizing organisms found in room samples at the genera level. Typically, for both infants, electronics had the lowest relative abundance of organisms detected in the gut whereas tubing had the highest. Temporal variation of gut genera was extreme in most environments.
The use of Bayesian microbial source tracking software , with the perspective of room samples as the source and fecal samples as the sink, produced mixed results in terms of finding likely gut reservoirs (Figure 6). In infant 1, tubing, surfaces, and electronics had the highest probabilities as sources, but the bloom of Bacteroides fragilis, from a source not detected by our sampling regime, lowered the probability of sampled source environments for the latter half of the sampling period. Infant 2’s samples showed the opposite pattern in that early gut colonists migrated from an unknown reservoir, whereas later in sampling, incubator, tubing, surfaces, and hands were the most probable reservoir.
Shared gut colonizers
The infant cohort shared only one gut colonizer, E. faecalis, which contained 100% 16S rRNA gene level sequence identity. A higher resolution analysis using a concatenated alignment of 32 highly conserved, single-copy genes show the strains differ by only 2 amino acids across the 4,101 positions. These two E. faecalis strains phylogenetically cluster most closely to each other, but are very closely related to other E. faecalis strains (Figure 7).
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.
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.
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.
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.
Kelley ST, Gilbert JA: Studying the microbiology of the indoor environment. Genome Biol. 2013, 14: 202-10.1186/gb-2013-14-2-202.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Singh N, Campbell J, Short B: Control of vancomycin-resistant enterococci in the neonatal intensive care unit. Infect Control. 2005, 26: 646-649.
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.
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.
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.
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.
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.
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.
Stackebrandt E, Goodfellow M: 16S/23S rRNA sequencing. Nucleic Acid Techniques in Bacterial Systematics. 1991, Chichester, UK: John Wiley & Son Ltd, 115-175.
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.
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.
Edgar RC: Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010, 26: 2460-2461. 10.1093/bioinformatics/btq461.
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.
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.
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.
Caporaso J, Kuczynski J, Stombaugh J: QIIME allows analysis of high-throughput community sequencing data. Nature. 2010, 7: 335-336.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32: 1792-1797. 10.1093/nar/gkh340.
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.
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.pdf
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.
Knights D, Kuczynski J, Charlson E: Bayesian community-wide culture-independent microbial source tracking. Nat Methods. 2011, 8: 6-10. 10.1038/nchembio.740.
Van Veen HW: Structural biology: Last of the multidrug transporters. Nature. 2010, 467: 926-927. 10.1038/467926a.
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.
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.
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.
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.
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-
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.
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.
Arias CA, Murray BE: The rise of the Enterococcus: beyond vancomycin resistance. Nat Rev Microbiol. 2012, 10: 266-278. 10.1038/nrmicro2761.
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.
Murray BE: The Life and Times of the Enterococcus. Clin Microbiol Rev. 1990, 3: 46-65.
Bradley CR, Fraise AP: Heat and chemical resistance of enterococci. J Hosp Infect. 1996, 34: 191-196. 10.1016/S0195-6701(96)90065-1.
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.
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.
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.
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.
Barnhart MM, Chapman MR: Curli biogenesis and function. Annu Rev Microbiol. 2006, 60: 131-147. 10.1146/annurev.micro.60.080805.142106.
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.
Funding was provided through the Alfred P Sloan Foundation and the National Science Foundation’s Graduate Research Fellowship Program.
The authors declare that they have no competing interests.
JFB, MJM, and BB conceived of the project. RB organized cohort recruitment and sample collections. BAF conducted nucleic acid extractions and BB conducted the amplification reactions. IS and BB conducted the metagenomic assemblies and BCT provided annotations. CSM and BB conducted the 16S rRNA gene reconstructions. BB and JFB wrote the final manuscript. All authors have read and approve the manuscript.