Mina Bashir1,2
Mahjabeen Ahmed1,3
Thomas Weinmaier4
Doina Ciobanu5
Natalia Ivanova5
Thomas R. Pieber2
Parag A. Vaishampayan1*- 1Biotechnology and Planetary Protection Group, Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA, USA
- 2Division of Endocrinology and Diabetology, Medical University of Graz, Graz, Austria
- 3Department of Biological Sciences, California State Polytechnic University, Pomona, CA, USA
- 4Division of Computational Systems Biology, Department of Microbiology and Ecosystem Science, University of Vienna, Vienna, Austria
- 5Department of Energy, Joint Genome Institute, Walnut Creek, CA, USA
Strict planetary protection practices are implemented during spacecraft assembly to prevent inadvertent transfer of earth microorganisms to other planetary bodies. Therefore, spacecraft are assembled in cleanrooms, which undergo strict cleaning and decontamination procedures to reduce total microbial bioburden. We wanted to evaluate if these practices selectively favor survival and growth of hardy microorganisms, such as pathogens. Three geographically distinct cleanrooms were sampled during the assembly of three NASA spacecraft: The Lockheed Martin Aeronautics' Multiple Testing Facility during DAWN, the Kennedy Space Center's Payload Hazardous Servicing Facility (KSC-PHSF) during Phoenix, and the Jet Propulsion Laboratory's Spacecraft Assembly Facility during Mars Science Laboratory. Sample sets were collected from the KSC-PHSF cleanroom at three time points: before arrival of the Phoenix spacecraft, during the assembly and testing of the Phoenix spacecraft, and after removal of the spacecraft from the KSC-PHSF facility. All samples were subjected to metagenomic shotgun sequencing on an Illumina HiSeq 2500 platform. Strict decontamination procedures had a greater impact on microbial communities than sampling location Samples collected during spacecraft assembly were dominated by Acinetobacter spp. We found pathogens and potential virulence factors, which determine pathogenicity in all the samples tested during this study. Though the relative abundance of pathogens was lowest during the Phoenix assembly, potential virulence factors were higher during assembly compared to before and after assembly, indicating a survival advantage. Decreased phylogenetic and pathogenic diversity indicates that decontamination and preventative measures were effective against the majority of microorganisms and well implemented, however, pathogen abundance still increased over time. Four potential pathogens, Acinetobacter baumannii, Acinetobacter lwoffii, Escherichia coli and Legionella pneumophila, and their corresponding virulence factors were present in all cleanroom samples. This is the first functional metagenomics study describing presence of pathogens and their corresponding virulence factors in cleanroom environments. The results of this study should be considered for microbial monitoring of enclosed environments such as schools, homes, hospitals and more isolated habitation such the International Space Station and future manned missions to Mars.
Introduction
Detection of signs of life on other planets is of particular interest for many of NASA's planetary missions. In order not to mistake earthborn microorganisms for unknown potential extraterrestrial life, planetary missions are subject to internationally accepted standards of planetary protection, established by the Committee of Space Research (COSPAR) (National Research Council, 2006). Jet Propulsion Laboratory's Planetary Protection Group has undertaken huge efforts (NASA Policy Directive (NPD) 8020.7G, 1999) to avoid inadvertent contamination of other planets with earthborn organisms, and to minimize the bioburden on spacecraft. Spores are of particular interest, given their high resistance to multiple sterilization techniques, including radiation (Venkateswaran et al., 2003; La Duc et al., 2007; Vaishampayan et al., 2012).
All spacecraft parts undergo extensive cleaning and sterilization steps, such as exposure to dry heat, vaporized hydrogen peroxide, radiation and alcohol on surfaces. Additional protocols to reduce the influx of particulate matter include daily vacuuming and mopping of floors, HEPA air filtration, regular replacement of tacky mats at all entry points, and strict gowning procedures. These precautions are routinely taken but with high frequency and stringency during the spacecraft assembly. All personnel that enter the cleanroom are required to put on cleanroom garments. This includes a full body suit, hair and beard nets, facemasks, additional head covering, gloves, shoe covers, and cleanroom boots. These are necessary measures since humans are the major source of contamination in these environments (La Duc et al., 2004; Probst et al., 2013). To monitor contamination levels, cleanrooms are regularly sampled for biological activity, particularly when spacecraft parts are being assembled and cleaned (La Duc et al., 2007; Vaishampayan et al., 2010a).
Multiple sterilization methods are chosen, because there is no known method that can eradicate all microbes, which is still compatible with spacecraft components. Only very resistant microorganisms, such as spores, pathogens, and extremophiles, can overcome these strict decontamination procedures (Ghosh et al., 2009; Derecho et al., 2014). Some microorganisms are even able to survive the harsh conditions of interstellar travel. Researchers placed spore-forming bacteria, isolated from cleanroom environment, outside the International Space Station for 18 months along with exposure to simulated Mars-like conditions, including atmospheric pressure and selective UV-radiation and some of them were still able to survive (Vaishampayan et al., 2012).
Our goal was to elucidate whether decontamination measures lead to selection of hardy microorganisms, including pathogens, in the cleanrooms and therefore posing a potential threat to human health. Pathogens might thrive in these environments perhaps due to their selective phenotypic characteristics, metabolic capabilities and reduced competition for scarce nutrients and niches. We were particularly interested in human pathogens, given that humans are the main source of contamination in cleanrooms (La Duc et al., 2004; Probst et al., 2013), and also because they are exposed to these constantly-evolving microbes. Most studies aiming at determining the microbiome of cleanrooms (La Duc et al., 2009; Sandle, 2011; Vaishampayan et al., 2013; Mahnert et al., 2015; Moissl-Eichinger et al., 2015), other indoor environments (Adams et al., 2015) or even the International Space Station (Checinska et al., 2015) have used 16S rRNA amplicon sequencing. 16S rRNA amplicon sequencing is often used to screen for potential pathogens (Case et al., 2007; Stadlbauer et al., 2015; Bashir et al., 2016). However, the lack of discriminability does not allow differentiating between potential and true pathogens. Previous functional metagenomic studies investigated pathogens in other indoor environments (Tringe et al., 2008; Afshinnekoo et al., 2015), but this is the first study, which focuses on the detection of pathogens as well as virulence factors in cleanrooms.
Three geographically distinct cleanrooms were sampled during the assembly of three NASA spacecraft: Phoenix in Cape Canaveral, Florida, DAWN in Fort Worth, Texas, and Mars Science Laboratory (Curiosity) in Pasadena, California. Sample sets from Phoenix mission were collected from the cleanroom at three time points: before arrival of the spacecraft, during the assembly and testing of the Phoenix spacecraft, and after removal of the spacecraft from the facility. All samples were subjected to whole metagenome shotgun sequencing on an Illumina HiSeq 2500 platform. We screened for pathogens and virulence factors, which determine pathogenicity. Clinically relevant pathogens were identified by searching taxonomic classification and potential virulence factors were identified by comparing reads to the Microbial virulence database.
Materials and Methods
Sample Collection and Processing
Multiple samples were collected from the floor of the Kennedy Space Center's Payload Hazardous Servicing Facility (KSC-PHSF), where the Phoenix spacecraft was assembled. Sample sets were collected from the KSC-PHSF surfaces at three time points: before arrival of the Phoenix spacecraft (10 samples; PHX-B), during the assembly and testing of the Phoenix spacecraft (8 samples; PHX-D), and after removal of the spacecraft from the KSC-PHSF facility (10 samples; PHX-A). 10 samples from the Lockheed Martin Aeronautics' Multiple Testing Facility (LMA-MTF) floor were collected during the DAWN spacecraft assembly. Samples were collected from the Ground Support Equipment (GSE) at Jet Propulsion Laboratory's spacecraft assembly facility (JPL-SAF) during the Mars Science Laboratory (2 samples; MSL) spacecraft assembly. These three cleanroom facilities were certified at ISO 8 (3,520,000 particles >0.5 μm m−3) level and maintained according to the standard cleaning practices. Each sample was collected from 1 m2 of the cleanroom floor or GSE by a wet surface sampling technique using Biological Sampling Kits (BiSKits, QuickSilver Analytics, Abingdon, Md.) and polyester wipes, respectively. Samples from each sampling event were concentrated using Amicon Ultra-15 centrifugal filter tube (Millipore, Jaffrey, NH, Ultracel-50 membrane) as described earlier (La Duc et al., 2009). DNA was extracted from each concentrated sample using bead beating and an automated DNA extraction instrument (Autolyser A-2 DNA, Axcyte Genomics, Menlo Park, CA) and pooled equimolar, as described earlier (Vaishampayan et al., 2010b). DNA samples were archived at −80°C until further use. Negative controls such as field control (sampling devise control), reagent control (during DNA extraction) at each step were collected. None of the negative controls had a sufficient DNA concentration for library preparation and were thus not included in further downstream analysis.
Metagenomic Sequencing
Sample processing was performed in a sodium hypochlorite (bleach) treated laminar flow hood in an ultra-clean environment. Operators were using single-use lab-coats, bleached gloves, hairnets, and booties. Due to low DNA concentrations, samples were subject to multiple displacement amplification (MDA) (Dean et al., 2002). Each sample was divided into 1 ml aliquots, which were amplified via MDA using Repli-g single-cell whole genome amplification kit (Qiagen part #150345). All plastic ware and water were ultraviolet (UV) treated in a Stratalinker 2400 UV Crosslinker (Stratagene, La Jolla, CA) with 254-nm UV for 30–90 min on ice (Woyke et al., 2011). This represents a UV dose range of 5.7–17.1 J/cm2, calculated by measuring the distance from inside the tubes to the light bulb (4 cm). Buffer and enzyme come pre-cleaned and don't require UV-radiation. MDA reaction was prepared following manufacturer protocol for single cells, scaling reaction volume down to 15 μl final volume and addition of Syto13 dye for real-time monitoring. MDA reaction was stopped when sample amplification reached saturation.
Amplified fractions of each sample were combined, and this pooled DNA product (100 μl) was sheared using a Covaris E210 instrument (Covaris, Woburn, MA) set to: 10% duty cycle, intensity 5, and 200 cycles per burst for 1 min. The concentration and fragment size of each sheared product was determined using Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA) in accordance with the manufacturer's recommended conditions. The sheared DNA was end-repaired, A-tailed, and ligated to Illumina adaptors according to standard Illumina PE protocols (Illumina, San Diego, CA). The concentration of the resulting Illumina-indexed libraries was again determined using Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA). Libraries were pooled and normalized to a final concentration of 400 mM each, and the primary bands corresponding to the sizes were gel-purified and dissolved in 30 μl TE. One flow-cell was generated from a pooled library, which was subsequently subjected to sequencing in an Illumina HiSeq2500 instrument (2 × 250 bp), in accordance with manufacturer-provided protocols. The raw sequence data are available within IMG/M (https://img.jgi.doe.gov/cgi-bin/mer/main.cgi) and NCBI's short read archive under the accession number SRP077843.
Sequence Data Analysis
We started with a total of 15,001,132 paired reads for PHX-B, 14,654,014 for PHX-D, and 22,355,430 for PHX-A before quality filtering and pairing. For MSL and DAWN we had 57,892,216 and 2,899,364 reads, respectively (Table 1).
Table 1. Data statistics: number of reads per sample starting with raw reads coming from the sequencer until final taxonomic and functional classification.
FastQC v0.10.1 (Andrews, 2010) was used to determine the base quality throughout the 250 bp HiSeq-generated paired-end reads. PEAR v0.9 (Zhang et al., 2014; default parameters) was used to merge paired reads. Unmerged forward and reverse reads were retained. Merged and unmerged reads were processed using prinseq-lite v0.20.3 (Schmieder and Edwards, 2011) with the following parameters: “-min_len 100 -trim_qual_right 20 -trim_qual_left 20 -trim_left 8.” Adapter sequences and overrepresented sequences were identified with FastQC and removed using Cutadapt v1.1 (Martin, 2011). PhiX174 and a JGI-standard collection of potential contaminant genomes (Supplementary Table 1) were removed by mapping trimmed high-quality reads using BBMap short read aligner v31.18 (Bushnell, 2014) to the respective genomes. Any reads matching any of these contaminant genomes were removed from the dataset.
To generate the human DNA sequence free dataset, all remaining high-quality reads were mapped with BBMap short read aligner against the human genome GRCh38 (including mitochondrial DNA). All positive matches were removed from the dataset.
Both, datasets including and excluding human DNA sequences were compared to NCBI non-redundant database using DIAMOND BLASTX v0.7.1 (Buchfink et al., 2014) with default parameters. Results were imported to MEGAN v5.10.5 (Huson et al., 2007; minimal bit score of 80%; “minscore 80”) for taxonomic binning, functional assignments to KEGG functions, and generation of rarefaction curves (phylogenetic diversity on genus level). After removal of unassigned and unclassified reads, taxonomy and KEGG pathways (Kanehisa and Goto, 2000; Kanehisa et al., 2014) were visualized using Krona Tools v2.4 (Ondov et al., 2011). Taxonomic and metabolic diversity calculations were done in QIIME 1.9.1 (Caporaso et al., 2011) with all samples subsampled to the smallest sample size observed.
Potential virulence factors were identified by comparing contaminant- and human-DNA-sequence-free reads to the Microbial Virulence Database MvirDB (Zhou et al., 2007) using DIAMOND BLASTX (Buchfink et al., 2014) with a 80% sequence similarity cut-off and maximum of target sequences of one. Sequences which passed these criteria were extracted and compared to NCBI non-redundant database using DIAMOND BLASTX (Buchfink et al., 2014) with a maximum of target sequences of one for virulence factor validation (Data Sheet 1 in Supplementary Material). Classified sequences were searched for clinically relevant pathogens (http://www.bode-science-center.com/center/relevant-pathogens-from-a-z.html accessed on Dec 1 2015; Supplementary Table 2).
Results
Phylogenetic Diversity of Cleanroom Samples
Alpha rarefaction curves indicate, besides sufficient sampling efforts, that diversity is drastically lower during the actual spacecraft assemblies (PHX-D, DAWN, and MSL) compared to before or after. This confirms that the very strict gowning, cleaning and sterilization procedures were well executed, and highly effective as previously described (Ghosh et al., 2010). MSL had the highest sampling depth but lowest bioburden (Figure 1A) as GSE undergo stringent cleaning procedures and are exposed to less handling and human contact compared to the cleanroom floors. Interestingly, microbial community profiles during active spacecraft assembly (PHX-D, DAWN, MSL) were more similar to each other than to samples from one location (Figure 1B). Moraxellaceae was the dominating family, making 83, 73, and 62% of all classified sequences for PHX-D, DAWN, and MSL, respectively. The majority of all Moraxellaceae, 94% to 100%, were Acinetobacter spp. (Figure 2 and Presentation S2 in Supplementary Material), making it the most dominating taxa during spacecraft assembly.
Figure 1. Diversity is lower during assembly. (A) Rarefaction curves of samples taken during assembly (PHX-D, DAWN, MSL) show less genera at the same sample size compared to samples taken before (PHX-B) and after (PHX-A) assembly. (B) Principal coordinates analysis on genus taxonomic level based on a Bray-Curtis dissimilarity matrix. Samples taken during spacecraft assembly show a similar community profile although the sampling locations were hundreds of miles apart.
Figure 2. Moraxellaceae dominate cleanroom during spacecraft assembly. Relative abundance of taxa at family level. Human, viral and archaeal impact was minimal. All taxa with a collective abundance of equal or less than 2% from all samples were combined in “Other”. See Presentation S1 in Supplementary Material for more details and lower taxonomic levels.
In general, bacteria were the most dominant kingdom present in all tested cleanrooms with 63 to >99% of all classified sequences. Archaea and viruses on the other hand accounted for less than 0.1% relative abundance combined (Figure 2 and Presentation S2 in Supplementary Material). Surprisingly, the amount of potentially human DNA was minimal. Only 0.04–2% of all sequences were classified as primates: PHX-B 2%, PHX-D 0.08%, PHX-A 0.2, DAWN 0.05, and MSL 0.2% (Presentation S2 in Supplementary Material).
In PHX-B eukaryotes made 36% of all classified sequences. Most of these sequences (22% of total; 60% of eukaryotes) belong to the class of arthropods, such as insects and arachnids. In all other samples less than 0.1% of all classified sequences were arthropods. Probably arthropod sequences originated from free DNA associated with dust particles, given that no living spiders or insects are present in any cleanrooms. In MSL all eukaryotic sequences were assigned to craniate. Fungi were also not prominent in our cleanrooms. The fungal abundance ranged from 0.0008% (MSL) to 1% (PHX-B) (Presentation S2 in Supplementary Material).
Metabolic Diversity During Spacecraft Assembly
Functional assignment resulted in 13,360 KEGG orthologous (KO) for PHX-B, 24,916 KOs for PHX-D, 298,350 KOs for PHX-A, 557 for DAWN and 664,699 for MSL (Table 1 and Data Sheet 1 in Supplementary Material). Figure 3 indicates that the majority of the functional classification was assigned to metabolism (PHX-B 67%, PHX-D 67%, PHX-A 90%, DAWN 75%, MSL 29%; Figure 3 and Data Sheet 1 in Supplementary Material). Although the percentage of sequences assigned to metabolism did not differ much across samples (Figure 3 and Data Sheet 1 in Supplementary Material), we saw a higher metabolic diversity during assembly compared to before or after spacecraft assembly samples (Table 1). We also found that the metabolism of pantothenate and coenzyme A is higher during assembly (PHX-D 4%, DAWN 2%) compared to PHX-B or PHX-A, 0.2% respectively (Data Sheet 1 in Supplementary Material). Nevertheless, no function associated with pantothenate and coenzyme A was found in MSL during assembly. Fifty-two percent of all functional classification from MSL was assigned to Holliday junction DNA helicase RuvB (Genetic Information Processing; Replication and Repair). In PHX-B, 9% were assigned to Genetic Information Processing, 12%, in PHX-D 3%, in PHX-A 11% in DAWN, and 52% in MSL (Figure 3).
Figure 3. Metabolic genes cover majority of functional classifications. KEGG Pathway analysis of PHX-B, PHX-D, PHX-A, DAWN and MSL. During assembly of PHX-D and DAWN, we see that a bigger fraction of all classified sequences have been assigned to metabolism. See Presentation S2 in Supplementary Material for more details.
Potential Pathogens in Cleanroom Samples
After taxonomic classification we selectively screened the classified binned sequences of all our samples for clinically relevant pathogens (Supplementary Table 2). In total we found 48 different human pathogens in all cleanrooms, responsible for various diseases, from gastrointestinal, to the nervous system. Twenty nine different pathogens were detected in PHX-B, 18 in PHX-D, 33 in PHX-A, 10 in DAWN and 11 in MSL (Table 2). Strikingly, four pathogens, namely Acinetobacter baumannii, Acinetobacter lwoffii, Escherichia coli and Legionella pneumophila, were detected in all cleanrooms, even though they are geographically separate. Besides these four pathogens, present in all samples, we found pathogens that were exclusive to KSC-PHSF, during all three time points, namely, Bacillus cereus, Burkholderia pseudomallei, Enterobacter cloacae, Enterococcus faecalis, Listeria monocytogenes, Pseudomonas aeruginosa, Staphylococcus aureus and Staphylococcus epidermidis. In case of PHX-D, 83% of all classified reads were identified as Acinetobacter spp. (Presentation S2 in Supplementary Material), which is listed as a clinically relevant pathogen (Supplementary Table 2). Additionally, we found a decreased pathogen diversity during the actual spacecraft assembly in KSC-PHSF, while pathogen abundance almost triples over time (PHX-B 1.52, PHX-D 2.34, and PHX-A 4.26%; Table 2).
Table 2. Pathogen diversity is lowest during assembly: pathogens found in the different cleanroom samples.
Pathogens and Corresponding Virulence Factors in Cleanrooms
Virulence factors are features, which distinguish pathogens from commensals or symbionts (Das et al., 2011). We found that the fraction of sequences identified as potential virulence factors increased over time in case of Phoenix (Table 3), although overall diversity was lower during assembly (Figure 1A). DAWN had approximately half the virulence factor fraction compared to Phoenix, but MSL, which was sampled from GSEs, had approximately 20 times less potential virulence factors compared to Phoenix.
Table 3. Accumulation of virulence factors over time: total number of virulence factors and hits normalized to hits per million reads found in cleanrooms.
To evaluate if we could find pathogens and their corresponding virulence factors, we identified potential virulence factors that are specifically associated with the pathogens found in our samples. We found 14 different potential virulence factors, which correspond to the classified pathogens in PHX-B, 48 for PHX-D and 41 for PHX-A. Nine different virulence factors were found to correspond with the classified pathogens in DAWN, and 6 were found for the MSL mission.
We were particularly interested in detecting potential virulence factors of the four pathogens, A. baumannii, A. lwoffii, E. coli, and L. pneumophila, which were found in all geographically separated cleanrooms. We found a Acinetobacter sp. specific aminoglycoside 6′-N-acetyltransferase lv and tetA, which are kanamycin B and tetracycline resistance genes, respectively, in PHX-D. Moreover, we found adeABC, which is an A. baumannii specific multidrug efflux pump and beta-lactamase TEM-1 (Supplementary Table 3).
We found that the abundance of potential virulence factors with associated pathogens increased over time (Table 4). Although, virulence factors diversity did not change over time, we observed a trend toward increased pathogens with associated virulence factors (Table 4, pathogenic diversity). Again, MSL had the smallest pathogenic diversity.
Table 4. Virulence factors with their corresponding pathogens.
Discussion
In this study, we demonstrated for the first time the presence of pathogens and their corresponding virulence factors in spacecraft assembly cleanrooms. Our approach allowed us not only to prove the presence of pathogens in the spacecraft assembly cleanrooms, but also their associated potential virulence factors. Most studies investigating the cleanroom microbiome have only used 16S rRNA amplicon sequencing (Vaishampayan et al., 2013; Mahnert et al., 2015). For example, the archived samples from KSC-PHSF during the Phoenix mission used in this study have previously been described using a cultivation based (Ghosh et al., 2010) and cultivation independent technique (Vaishampayan et al., 2010a). On one hand, cultivation based techniques offer a very limited insight into the wide spectrum of microbial diversity, given that most microorganisms are not cultivable, while 16S rRNA amplicon sequencing on the other hand shows a more broad picture, but does not allow a reliable phylogenetic classification below genus level and does not provide any information regarding virulence factors and potential pathogenicity.
We observed that cleanroom samples are dominated by bacteria as reported previously (Weinmaier et al., 2015). Contrary to previous studies, which found substantially more human, archaeal and viral sequences in cleanrooms (Moissl-Eichinger, 2011; Weinmaier et al., 2015), we found significantly less of each taxon in all cleanrooms tested during this study. These previous studies have sampled uncontrolled gowning area and ISO-8 cleanrooms, where no active spacecraft assembly was ongoing. Moreover, each cleanroom is unique, because of factors such as geographical location (Moissl et al., 2007), assembly activities, different decontamination procedures and most importantly, different workers, which are the main source of contamination.
We saw an increased metabolic diversity in samples collected from cleanrooms during spacecraft assembly. Cleanrooms are sometimes referred to as extreme environments; not only due to strict decontamination procedures, but also due to the lack of nutrients, water and cofactors (La Duc et al., 2007; Ghosh et al., 2010). Since there are very few resources to rely on in an area that is maintained to be uninhabitable, any microbes able to survive here would have to synthesize all necessary factors themselves. Sterilization procedures and gowning requirements are even stricter during assembly, making it even harder for microorganisms to survive. Strict gowning protocols and increased frequency of cleaning leads to decrease in introduction of human associated microbes, despite high human activities in the cleanroom during assembly. This might also explain lower phylogenetic diversity during Phoenix spacecraft assembly compared to before or after assembly. In addition to the decreased phylogenetic diversity, pathogenic diversity was also lower during spacecraft assembly, however, we observed an increase in pathogen abundance over time. This suggests that strict decontamination procedure favor the growth of pathogens. Nevertheless, studies with bigger sample sizes need to confirm our descriptive findings. A considerable amount, in case of MSL more than 50%, of all reads was assigned to genetic information processing. This highlights the importance of genetic information processing, including DNA repair in such a harsh environment. Surprisingly, microbial profiles during assembly were very similar. Although, DAWN and MSL samples were collected from geographically distinct locations, they were more similar to PHX-D than PHX-B or PHX-A. This suggests that decontamination procedures have a bigger effect on the cleanroom microbiome than location. Taken together, these results show that decontamination and gowning measures were not only sufficient, but also well executed.
Most virulence factors are organized in so-called pathogenicity-islands (Schmidt and Hensel, 2004). Commensals can turn into pathogens by acquiring pathogenicity-island through phages, or horizontal gene transfer. For example, wild type Vibrio cholerae are not able to cause deadly diarrhea. Only upon infection by the CTX prophage they acquire a pathogenicity island coding for virulence factors, such as the cholera toxin and pili, needed for recognition host and disease induction (Das et al., 2011). Therefore, virulence factor detection is the only reliable method to identify pathogens.
Acinetobacter baumannii, Acinetobacter lwoffii, Escherichia coli and Legionella pneumophila were found in all samples, although samples were collected from three geographically distinct sites. These prevalent pathogens have to be very resistant to overcome all the cleaning and decontamination procedures. Acinetobacter spp., such as A. baumannii and A. lwoffii are non-fastidious and can rely on a single energy source from different substrates (Rathinavelu et al., 2003). They are resistant to radiation (Firstenberg-Eden et al., 1980b) and several disinfectants and can survive in a wide range of temperatures (Firstenberg-Eden et al., 1980a) and even in low pH. These features might explain why Acinetobacter spp. were the most dominating species during spacecraft assembly in this study. Acinetobacter spp. have also been reported in high abundance in cleanrooms in previous studies (Vaishampayan et al., 2010a; Mahnert et al., 2015). Acinetobacter baumannii has been isolated from water and soil (Yeom et al., 2013), but it has also been found in other hostile environments such as intensive care units. Although A. baumannii is not pathogenic to healthy individuals, it can be an opportunistic pathogen in immunocompromised patients. Hence, it is one of the ESKAPE pathogens (Boucher et al., 2009), which are multidrug-resistant bacteria, responsible for the majority of nosocomial infections (Rice, 2008).
We found A. baumannii specific beta-Lactamase TEM-1, AdeABC and another cation/multidrug efflux pump, which might be responsible for A. baumannii's resistance against all decontamination measures. AdeABC alone is responsible for resistance to aminoglycosides, tetracyclines, erythromycin, chloramphenicol, trimethoprim, fluoroquinolones, some beta-lactams, and also recently tigecycline since they have been described as substrates for this multidrug efflux pump (Wieczorek et al., 2008). We did not detect A. Iwoffii associated potential virulence factors in our data set. MvirDB has only three A. Iwoffii (formerly known as Acinetobacter calcoaceticus) associated virulence factors (two beta-lactamases and a chloramphenicol acetyl transferase). Nevertheless, the presence of this opportunistic pathogen in all our sample collection from locations separated by hundreds of miles, its resistant features, and our finding that Acinetobacter spp. were dominating in all three locations during assembly, suggests that A. lwoffii and A. baumannii are actually viable in the spacecraft cleanroom environment. L. pneumophila, another pathogen present in all three distinct locations, is the causative agent of the Legionnaires' disease (Shevchuk et al., 2011), with symptoms such as fever, chills, and coughing. We found Legionella secretion pathway protein E (LspE), which is part of a type II secretion system required for its full virulence and environmental persistence (Hales and Shuman, 1999). In addition, other L. pneumophila associated virulence factors, such as the catalase-peroxidase KatB and superoxide dismutase were present, potentially explaining why this species is resistant to hydrogen peroxide treatment, one of the decontamination procedures. The last potential pathogen we found in all cleanrooms was E. coli. Although, we have found several virulence factors such as transposases and antimicrobial resistance genes, we cannot confirm whether or not this specific E. coli is a pathogen, given that more and more antimicrobial resistance genes are being found in commensal E. coli (Kaesbohrer et al., 2012; Tadesse et al., 2012; Wasyl et al., 2013). While we think that the four pathogens found in all geographically separated cleanrooms are alive, given their resistant features, we are not able to tell if the classified taxa and functions derive from intact living cells or if they are originating from dead cells. In an ongoing study we're including propidium monoazide staining, enabling us to differentiate between sequences coming from intact live and dead microorganisms.
Interestingly, potential virulence factor abundance increased over time, despite lower phylogenetic diversity during assembly. We have concluded that virulence factors may provide a survival advantage in this very hostile environment (Rathinavelu et al., 2003). Multidrug efflux pumps might be pumping out harmful compounds before they are able to execute their deadly effect (Yoon et al., 2013). This virulence factor accumulation seems to be species dependent, as we also see an increase in pathogenic abundance over time. We also found pathogens not belonging to the bacterial kingdom; such as Candida parapsilosis, a fungus, which plays an important role in wound and tissue sepsis of immunocompromised patients and makes up to 15% of all Candida infections.
One limitation of this study is the low biomass in cleanroom samples, due to the repeated strict cleaning and decontamination practices. MDA was necessary to acquire DNA concentrations sufficient for library preparation. MDA can introduce bias, by favoring some DNA fragments over others (Direito et al., 2014). Therefore, some microorganisms might not have been detected in our approach, while others might be overrepresented. Although, we performed stringent quality filtering of our reads, it's impossible to get rid of all errors and biases. Homology based approaches such as BLASTx against specialized databases such as MvirDB are biased, because a sequence with an 80% sequence similarity might have a better hit to a reference which is not in the database. However, we searched all positive MvirDB hits against NCBI non-redundant database, and the majority was classified as virulence factors (see Presentation S1 in Supplementary Material). Moreover, the circumstantial evidences of the presence of virulence factors associated with human pathogens in cleanroom samples should be confirmed by implementation of selective cultivation based approaches and viability-based molecular assays in future missions.
Humans spend most of their lives indoors (Höppe and Martinac, 1998). Recent studies have speculated that human microbiome is the major contributor to the overall indoor microbiome (Lax et al., 2014). Stringent cleaning and maintenance practices in highly controlled indoor environments such as cleanrooms, hospitals and intensive care units may lead to a relative increase of human pathogens in these environments. This may have serious impact on health of the inhabitants. Monitoring pathogens and virulence factors in these indoor environments may prevent diseases such as nosocomial infections and sustain human health.
The results of this study could be used to develop fast and cost-efficient tests (Craw et al., 2015) to detect the presence of specific pathogens or their virulence factors in enclosed environments such as public transport, pharmaceutical cleanrooms, hospitals, and intensive care units. This study has broadened our understanding of the role of pathogens in such highly controlled environments and should be considered for microbial monitoring of the ISS during sustained presence of humans in space and future manned missions to Mars.
Author Contributions
Designed project: PV. Performed wet or computational experiments: MB, PV, DC and NI. Analyzed data: MB, NI and TW. Drafting the manuscript: MB and MA. Generated figures and tables: MB and MA. Wrote and critically reviewed the manuscript: all authors.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
Part of the research described in this study was carried out at the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration. This research was funded by Planetary Protection Research program in NNH11ZDA001N, ROSES 2011 awarded to PV and NI. The authors are grateful to Drs. Catharine Conley and Ying Lin for valuable discussion and input. We would like to thank Dr. Kasthuri Venkateswaran (JPL) for making archived DNA available. MB is thankful for the financial support of the FWF Austrian Science Fund (W1241) and the Medical University of Graz through the PhD Program Molecular Fundamentals of Inflammation (DK-MOLIN), as well as the Bank Austria Visiting Scientist Program of the Medical University of Graz.
Supplementary Material
The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fmicb.2016.01321
References
Adams, R. I., Bateman, A. C., Bik, H. M., and Meadow, J. F. (2015). Microbiota of the indoor environment: a meta-analysis. Microbiome 3, 49. doi: 10.1186/s40168-015-0108-3
Afshinnekoo, E., Meydan, C., Levy, S., and Mason, C. E. (2015). Geospatial resolution of human and bacterial diversity with city-scale metagenomics. Cell Syst. 1, 1–15. doi: 10.1016/j.cels.2015.01.001
Andrews, S. (2010). FastQC: a Quality Control Tool for High Throughput Sequence Data. Available online at: http://www.citeulike.org/user/nailest/article/11583827
Bashir, M., Prietl, B., Tauschmann, M., Mautner, S. I., Kump, P. K., Treiber, G., et al. (2016). Effects of high doses of vitamin D3 on mucosa-associated gut microbiome vary between regions of the human gastrointestinal tract. Eur. J. Nutr. 55, 1479–1489. doi: 10.1007/s00394-015-0966-2
Boucher, H. W., Talbot, G. H., Bradley, J. S., Edwards, J. E., Gilbert, D., Rice, L. B., et al. (2009). Bad Bugs, no drugs: no ESKAPE! an update from the Infectious Diseases Society of America. Clin. Infect. Dis. 48, 1–12. doi: 10.1086/595011
Buchfink, B., Xie, C., and Huson, D. H. (2014). Fast and sensitive protein alignment using DIAMOND. Nat. Methods 12, 59–60. doi: 10.1038/nmeth.3176
Bushnell, B. (2014). BBMap Short Read Aligner, and Other Bioinformatics Tools. Available online at: http://www.osti.gov/scitech/biblio/1241166
Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2011). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. doi: 10.1038/nmeth.f.303
Case, R. J., Boucher, Y., Dahllöf, I., Holmström, C., Doolittle, W. F., and Kjelleberg, S. (2007). Use of 16S rRNA and rpoB genes as molecular markers for microbial ecology studies. Appl. Environ. Microbiol. 73, 278–288. doi: 10.1128/AEM.01177-06
Checinska, A., Probst, A. J., Vaishampayan, P., White, J. R., Kumar, D., Stepanov, V. G., et al. (2015). Microbiomes of the dust particles collected from the International Space Station and Spacecraft Assembly Facilities. Microbiome 3, 50. doi: 10.1186/s40168-015-0116-3
Craw, P., Mackay, R. E., Naveenathayalan, A., Hudson, C., Branavan, M., Sadiq, S. T., et al. (2015). A simple, low-cost platform for real-time isothermal nucleic acid amplification. Sensors (Basel). 15, 23418–23430. doi: 10.3390/s150923418
Das, B., Bischerour, J., and Barre, F.-X. (2011). Molecular mechanism of acquisition of the cholera toxin genes. Indian J. Med. Res. 133, 195–200.
Dean, F. B., Hosono, S., Fang, L., Wu, X., Faruqi, A. F., Bray-Ward, P., et al. (2002). Comprehensive human genome amplification using multiple displacement amplification. Proc. Natl. Acad. Sci. U.S.A. 99, 5261–5266. doi: 10.1073/pnas.082089499
Derecho, I., McCoy, K. B., Vaishampayan, P., Venkateswaran, K., and Mogul, R. (2014). Characterization of Hydrogen Peroxide-Resistant Acinetobacter Species Isolated during the Mars Phoenix Spacecraft Assembly. Astrobiology 14, 837–847. doi: 10.1089/ast.2014.1193
Direito, S. O. L., Zaura, E., Little, M., Ehrenfreund, P., and Röling, W. F. M. (2014). Systematic evaluation of bias in microbial community profiles induced by whole genome amplification. Environ. Microbiol. 16, 643–657. doi: 10.1111/1462-2920.12365
Firstenberg-Eden, R., Rowley, D. B., and Shattuck, E. (1980a). Thermal inactivation and injury of Moraxella-Acinetobacter cells in ground beef. Appl. Environ. Microbiol. 39, 159–164.
Firstenberg-Eden, R., Rowley, D. B., and Shattuck, G. E. (1980b). Factors affecting inactivation of Moraxella-Acinetobacter cells in an irradiation process. Appl. Environ. Microbiol. 40, 480–485.
Ghosh, S., Osman, S., Vaishampayan, P., and Venkateswaran, K. (2010). Recurrent isolation of extremotolerant bacteria from the clean room where Phoenix spacecraft components were assembled. Astrobiology 10, 325–335. doi: 10.1089/ast.2009.0396
Ghosh, S., Zhang, P., Li, Y. Q., and Setlow, P. (2009). Superdormant spores of Bacillus species have elevated wet-heat resistance and temperature requirements for heat activation. J. Bacteriol. 191, 5584–5591. doi: 10.1128/JB.00736-09
Hales, L. M., and Shuman, H. A. (1999). Legionella pneumophila contains a type II general secretion pathway required for growth in amoebae as well as for secretion of the Msp protease. Infect. Immun. 67, 3662–3666.
Höppe, P., and Martinac, I. (1998). Indoor climate and air quality. Review of current and future topics in the field of ISB study group 10. Int. J. Biometeorol. 42, 1–7.
Huson, D. H., Auch, A. F., Qi, J., and Schuster, S. C. (2007). MEGAN analysis of metagenomic data. Genome Res. 17, 377–386. doi: 10.1101/gr.5969107
Kaesbohrer, A., Schroeter, A., Tenhagen, B. A., Alt, K., Guerra, B., and Appel, B. (2012). Emerging antimicrobial resistance in commensal Escherichia coli with public health relevance. Zoonoses Public Health 59, 158–165. doi: 10.1111/j.1863-2378.2011.01451.x
Kanehisa, M., and Goto, S. (2000). KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 28, 27–30. doi: 10.1093/nar/28.1.27
Kanehisa, M., Goto, S., Sato, Y., Kawashima, M., Furumichi, M., and Tanabe, M. (2014). Data, information, knowledge and principle: back to metabolism in KEGG. Nucleic Acids Res. 42, D199–D205. doi: 10.1093/nar/gkt1076
La Duc, M., Kern, R., Venkateswaran, K., and La Duc, M. T. (2004). Microbial monitoring of spacecraft and associated environments. Microb. Ecol. 47, 150–158. doi: 10.1007/s00248-003-1012-0
La Duc, M. T., Dekas, A., Osman, S., Moissl, C., Newcombe, D., and Venkateswaran, K. (2007). Isolation and characterization of bacteria capable of tolerating the extreme conditions of clean room environments. Appl. Environ. Microbiol. 73, 2600–2611. doi: 10.1128/AEM.03007-06
La Duc, M. T., Osman, S., Vaishampayan, P., Piceno, Y., Andersen, G., Spry, J. A., et al. (2009). Comprehensive census of bacteria in clean rooms by using DNA microarray and cloning methods. Appl. Environ. Microbiol. 75, 6559–6567. doi: 10.1128/AEM.01073-09
Lax, S., Smith, D. P., Hampton-Marcell, J., Owens, S. M., Handley, K. M., Scott, N. M., et al. (2014). Longitudinal analysis of microbial interaction between humans and the indoor environment. Science 345, 1048–1052. doi: 10.1126/science.1254529
Mahnert, A., Vaishampayan, P., Probst, A. J., Auerbach, A., Moissl-Eichinger, C., Venkateswaran, K., et al. (2015). Cleanroom maintenance significantly reduces abundance but not diversity of indoor microbiomes. PLoS ONE 10:e0134848. doi: 10.1371/journal.pone.0134848
Martin, M. (2011). Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet. J. 17, 10–12. doi: 10.14806/ej.17.1.200
Moissl, C., Osman, S., La Duc, M. T., Dekas, A., Brodie, E., DeSantis, T., et al. (2007). Molecular bacterial community analysis of clean rooms where spacecraft are assembled. FEMS Microbiol. Ecol. 61, 509–521. doi: 10.1111/j.1574-6941.2007.00360.x
Moissl-Eichinger, C. (2011). Archaea in artificial environments: their presence in global spacecraft clean rooms and impact on planetary protection. ISME J. 5, 209–219. doi: 10.1038/ismej.2010.124
Moissl-Eichinger, C., Auerbach, A. K., Probst, A. J., Mahnert, A., Tom, L., Piceno, Y., et al. (2015). Quo vadis? Microbial profiling revealed strong effects of cleanroom maintenance and routes of contamination in indoor environments. Sci. Rep. 5:9156. doi: 10.1038/srep09156
National Research Council (2006). Preventing the Forward Contamination of Mars. Washington, DC: National Academies Press.
Ondov, B. D., Bergman, N. H., and Phillippy, A. M. (2011). Interactive metagenomic visualization in a Web browser. BMC Bioinformatics 12:385. doi: 10.1186/1471-2105-12-385
Probst, A. J., Auerbach, A. K., and Moissl-Eichinger, C. (2013). Archaea on Human Skin. PLoS ONE 8:e65388. doi: 10.1371/journal.pone.0065388
Rathinavelu, S., Zavros, Y., and Merchant, J. L. (2003). Acinetobacter lwoffii infection and gastritis. Microbes Infect. 5, 651–657. doi: 10.1016/S1286-4579(03)00099-6
Rice, L. B. (2008). Federal funding for the study of antimicrobial resistance in nosocomial pathogens : no ESKAPE. J. Infect. Dis. 197, 1079–1081. doi: 10.1086/533452
Sandle, T. (2011). A review of cleanroom Microflora: types, trends, and patterns. PDA J. Pharm. Sci. Technol. 65, 392–403. doi: 10.5731/pdajpst.2011.00765
Schmidt, H., and Hensel, M. (2004). pathogenicity islands in bacterial pathogenesis pathogenicity islands in bacterial pathogenesis. Clin. Microbiol. Rev. 17, 14–56. doi: 10.1128/CMR.17.1.14-56.2004
Schmieder, R., and Edwards, R. (2011). Quality control and preprocessing of metagenomic datasets. Bioinformatics 27, 863–864. doi: 10.1093/bioinformatics/ btr026
Shevchuk, O., Jäger, J., and Steinert, M. (2011). Virulence properties of the legionella pneumophila cell envelope. Front. Microbiol. 2:74. doi: 10.3389/fmicb.2011.00074
Stadlbauer, V., Leber, B., Lemesch, S., Trajanoski, S., Bashir, M., Horvath, A., et al. (2015). Lactobacillus casei Shirota supplementation does not restore Gut Microbiota composition and Gut Barrier in Metabolic syndrome: a randomized pilot study. PLoS ONE 10:e0141399. doi: 10.1371/journal.pone.0141399
Tadesse, D. A., Zhao, S., Tong, E., Ayers, S., Singh, A., Bartholomew, M. J., et al. (2012). Antimicrobial drug resistance in Escherichia coli from humans and food animals, United States, 1950-2002. Emerg. Infect. Dis. 18, 741–749. doi: 10.3201/eid1805.111153
Tringe, S. G., Zhang, T., Liu, X., Yu, Y., Lee, W. H., Yap, J., et al. (2008). The airborne metagenome in an indoor urban environment. PLoS ONE 3:e1862. doi: 10.1371/journal.pone.0001862
Vaishampayan, P., Osman, S., Andersen, G., and Venkateswaran, K. (2010a). High-density 16S microarray and clone library-based microbial community composition of the Phoenix spacecraft assembly clean room. Astrobiology 10, 499–508. doi: 10.1089/ast.2009.0443
Vaishampayan, P., Probst, A. J., La Duc, M. T., Bargoma, E., Benardini, J. N., Andersen, G. L., et al. (2013). New perspectives on viable microbial communities in low-biomass cleanroom environments. ISME J. 7, 312–324. doi: 10.1038/ismej.2012.114
Vaishampayan, P., Probst, A., Krishnamurthi, S., Ghosh, S., Osman, S., McDowall, A., et al. (2010b). Bacillus horneckiae sp. nov., isolated from a spacecraft-assembly clean room. Int. J. Syst. Evol. Microbiol. 60, 1031–1037. doi: 10.1099/ijs.0.008979-0
Vaishampayan, P. A., Rabbow, E., Horneck, G., and Venkateswaran, K. J. (2012). Spores for a prolonged period of time in real space conditions. Astrobiology 12, 487–497. doi: 10.1089/ast.2011.0738
Venkateswaran, K., Kempf, M., Chen, F., Satomi, M., Nicholson, W., and Kern, R. (2003). Bacillus nealsonii sp. nov., isolated from a spacecraft-assembly facility, whose spores are gamma-radiation resistant. Int. J. Syst. Evol. Microbiol. 53, 165–172. doi: 10.1099/ijs.0.02311-0
Wasyl, D., Hoszowski, A., Zając, M., and Szulowski, K. (2013). Antimicrobial resistance in commensal Escherichia coli isolated from animals at slaughter. Front. Microbiol. 4:221. doi: 10.3389/fmicb.2013.00221
Weinmaier, T., Probst, A. J., La Duc, M. T., Ciobanu, D., Cheng, J.-F., Ivanova, N., et al. (2015). A viability-linked metagenomic analysis of cleanroom environments: eukarya, prokaryotes, and viruses. Microbiome 3, 62. doi: 10.1186/s40168-015-0129-y
Wieczorek, P., Sacha, P., Hauschild, T., Zórawski, M., Krawczyk, M., and Tryniszewska, E. (2008). Multidrug resistant Acinetobacter baumannii–the role of AdeABC (RND family) efflux pump in resistance to antibiotics. Folia Histochem. Cytobiol. 46, 257–267. doi: 10.2478/v10042-008-0056-x
Woyke, T., Sczyrba, A., Lee, J., Rinke, C., Tighe, D., Clingenpeel, S., et al. (2011). Decontamination of MDA reagents for single cell whole genome amplification. PLoS ONE 6:e26161. doi: 10.1371/journal.pone.0026161
Yeom, J., Shin, J.-H., Yang, J.-Y., Kim, J., and Hwang, G.-S. (2013). (1)H NMR-based metabolite profiling of planktonic and biofilm cells in Acinetobacter baumannii 1656-2. PLoS ONE 8:e57730. doi: 10.1371/journal.pone.0057730
Yoon, E.-J., Courvalin, P., and Grillot-Courvalin, C. (2013). RND-type efflux pumps in multidrug-resistant clinical isolates of Acinetobacter baumannii: major role for AdeABC overexpression and AdeRS mutations. Antimicrob. Agents Chemother. 57, 2989–2995. doi: 10.1128/AAC.02556-12
Zhang, J., Kobert, K., Flouri, T., and Stamatakis, A. (2014). PEAR: a fast and accurate Illumina Paired-End reAd mergeR. Bioinformatics 30, 614–620. doi: 10.1093/bioinformatics/btt593
Keywords: cleanroom, pathogens, indoor environments, microbiome, spacecraft, virulence factors, Acinetobacter, functional metagenomics
Citation: Bashir M, Ahmed M, Weinmaier T, Ciobanu D, Ivanova N, Pieber TR and Vaishampayan PA (2016) Functional Metagenomics of Spacecraft Assembly Cleanrooms: Presence of Virulence Factors Associated with Human Pathogens. Front. Microbiol. 7:1321. doi: 10.3389/fmicb.2016.01321
Received: 09 March 2016; Accepted: 10 August 2016;
Published: 09 September 2016.
Edited by:
Martin Grube, University of Graz, AustriaReviewed by:
Tomislav Cernava, Austrian Centre of Industrial Biotechnology, AustriaTim Sandle, University of Manchester, UK
Copyright © 2016 Bashir, Ahmed, Weinmaier, Ciobanu, Ivanova, Pieber and Vaishampayan. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Parag A. Vaishampayan, vaishamp@jpl.nasa.gov