Front. Microbiol. Frontiers in Microbiology Front. Microbiol. 1664-302X Frontiers Media S.A. 10.3389/fmicb.2021.608478 Microbiology Original Research Characterization of Spacesuit Associated Microbial Communities and Their Implications for NASA Missions Danko David 1 2 Malli Mohan Ganesh Babu 3 Sierra Maria A. 1 4 Rucker Michelle 5 Singh Nitin K. 1 Regberg Aaron B. 6 Bell Mary S. 7 O’Hara Niamh B. 2 Ounit Rachid 8 Mason Christopher E. 2 4 9 10 * Venkateswaran Kasthuri 3 * 1Tri-Institutional Computational Biology & Medicine Program, Weill Cornell Medicine of Cornell University, Manhattan, NY, United States 2The HRH Prince Alwaleed Bin Talal Bin Abdulaziz Alsaud Institute for Computational Biomedicine, Weill Cornell Medicine, New York, NY, United States 3Biotechnology and Planetary Protection Group, Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA, United States 4Department of Physiology and Biophysics, Weill Cornell Medicine, New York, NY, United States 5Exploration Mission Planning Office, Johnson Space Center, Houston, TX, United States 6Astromaterials Research and Exploration Science Division, Johnson Space Center, Houston, TX, United States 7Jacobs@NASA/Johnson Space Center, Houston, TX, United States 8Department of Computer Science and Engineering, University of California, Riverside, Riverside, CA, United States 9The WorldQuant Initiative for Quantitative Prediction, Weill Cornell Medicine, New York, NY, United States 10The Feil Family Brain and Mind Research Institute, Weill Cornell Medicine, New York, NY, United States

Edited by: Andreas Teske, University of North Carolina at Chapel Hill, United States

Reviewed by: Tim Sandle, The University of Manchester, United Kingdom; Rakesh Mogul, California State Polytechnic University, Pomona, United States

*Correspondence: Christopher E. Mason, chm2042@med.cornell.edu Kasthuri Venkateswaran, kjvenkat@jpl.nasa.gov

These authors have contributed equally to this work and share first authorship

This article was submitted to Extreme Microbiology, a section of the journal Frontiers in Microbiology

29 07 2021 2021 12 608478 20 09 2020 16 06 2021 Copyright © 2021 Danko, Malli Mohan, Sierra, Rucker, Singh, Regberg, Bell, O’Hara, Ounit, Mason and Venkateswaran. 2021 Danko, Malli Mohan, Sierra, Rucker, Singh, Regberg, Bell, O’Hara, Ounit, Mason and Venkateswaran

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) and the copyright owner(s) 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.

Background

Crewed National Aeronautics and Space Administration (NASA) missions to other solar system bodies are currently being planned. One high-profile scientific focus during such expeditions would be life detection, specifically the discovery of past or present microbial life, if they exist. However, both humans and associated objects typically carry a high microbial burden. Thus, it is essential to distinguish between microbes brought with the expedition and those present on the exploring planets. Modern spacesuits are unique, customized spacecraft which provide protection, mobility and life support to crew during spacewalks, yet they vent, and the mobility of microbes through spacesuits has not been studied.

Results

To evaluate the microbial colonization of spacesuits, NASA used an Extravehicular Activity swab kit to examine viable microbial populations of 48 samples from spacesuits using both traditional microbiological methods and molecular sequencing methods. The cultivable microbial population ranged from below the detection limit to 9 × 102 colony forming units per 25 cm2 of sample and also significantly varied by the location. The cultivable microbial diversity was dominated by members of Bacillus, Arthrobacter, and Ascomycota. However, 16S rRNA-based viable bacterial burden ranged from 105 to 106 copies per 25 cm2 of sample. Shotgun metagenome sequencing revealed the presence of a diverse microbial population on the spacesuit surfaces, including Curtobacterium and Methylobacterium from across all sets of spacesuits in high abundance. Among bacterial species identified, higher abundance of Cutibacterium acnes, Methylobacterium oryzae, and M. phyllosphaerae reads were documented.

Conclusion

The results of this study provide evidence that identical microbial strains may live on the wrist joint, inner gauntlet, and outer gauntlet of spacesuits. This raises the possibility, but does not confirm that microbial contaminants on the outside of the suits could contaminate planetary science operations unless additional measures are taken. Overall, these data provide the first estimate of microbial distribution associated with spacesuit surfaces, which will help future mission planners develop effective planetary protection strategies.

spacesuit microbial diversity ISS metagemonic metagenome assembled genomes (MAGs) microbial ecology

香京julia种子在线播放

    1. <form id=HxFbUHhlv><nobr id=HxFbUHhlv></nobr></form>
      <address id=HxFbUHhlv><nobr id=HxFbUHhlv><nobr id=HxFbUHhlv></nobr></nobr></address>

      Introduction

      Several spacefaring nations and private corporations are planning to send humans and spacecraft to other planets such as Mars, to search for evidence of habitats that could support life (NRC, 2014). Planetary Protection research efforts at National Aeronautics and Space Administration (NASA) seek to develop technologies to minimize any terrestrial microbial contamination to ensure the safety and health of astronauts, while also preserving scientific integrity of exoplanetary samples (NASA, 2019a). Planetary Protection aims involve the study and prevention of forward and back contamination, meaning the interchange of microbes and organic materials from Earth to other solar system bodies and vice versa (Debus and Arnould, 2008).

      When astronauts will be sent to search for life on other planets, it will be necessary to understand what microorganisms they may bring with them. It is estimated that 85% of all microbial isolates recovered from spacecraft and supported facilities are microorganisms associated with the human microbiota (Nicholson et al., 2009). Accordingly, a team at the Johnson Space Center (JSC) at NASA has developed a prototype Extravehicular Activity (EVA) swab kit that is suitable for handling by the astronauts in spacesuits to collect microbial samples aseptically, aiming to profile microorganisms associated with spacesuits (Rucker et al., 2018). In this communication, a microbial characterization associated with wrist joints of flight Extravehicular Mobility Unit (EMU), Modified Advanced Crew Escape System and Orion Crew Survival System (MACES/OCSS) spacesuits was carried out to evaluate the form, fit and function of the EVA swab tool; that functional testing provided an opportunity to characterize the typical microbial contamination on spacesuits.

      To explore and work in space, crew members must take their environment with them because there is no atmospheric pressure and no oxygen to sustain life. Inside the human crew vehicle, the atmosphere can be controlled so that special clothing is not necessary, but when outside exploring in space, astronauts need protection (Schwartz et al., 2002). Since various materials including fabrics and clothing are known to harbor specific microbiomes (Breuker et al., 2003; Cappitelli and Sorlini, 2008; Cataño et al., 2012; Callewaert et al., 2014; Sterndorff et al., 2020), it is of the highest interest to the NASA scientific community to explore the microbiome of the spacesuit (National Academies of Sciences, Engineering and Medicine, 2018). This study is not designed to understand the indigenous microbiome of the spacesuit when manufactured; instead spacesuit microbiome was measured when crew wear them after nominal handling and use to see how microorganisms might persist on the suits. Thorough characterization of spacesuit microbiome will enable the design of appropriate spacesuit architecture to minimize human commensal microorganism, which cannot be sterilized, from leaking into the external environment thus compromising life detection missions. Currently, all NASA spacesuits are designed to be flexible and which could lead to leakage. However, leak paths are not well-characterized, and it remains unclear what fraction of leakage occurs through mechanisms that would transport microbes. Characterization of spacesuits will also allow NASA to better understand cleaning process effectiveness for the spacesuits.

      Since 2006, the field of genomics has been revolutionized by the development of next-generation sequencing technologies, enabling the comprehensive understanding of the microbial ecology of built environments such as offices (Chase et al., 2016), hospitals (Westwood et al., 2014), and transportation system environments (Hsu et al., 2016; Danko et al., 2021a) where humans spend a significant fraction of their time. Subsequently, molecular microbial community analyses were implemented to monitor the International Space Station (ISS) (Singh et al., 2018; Checinska Sielaff et al., 2019) and spacecraft assembly cleanrooms (Danko et al., 2021b) but this is the first report measuring spacesuit microbiome. While these technologies for microbial detection have elucidated the prevalence of microbial species, it was not until recently significant efforts have been pointed at developing sampling methods that enable sample collection in microgravity or a vacuum, that are simple to handle by crew members donned with large gloves, and that could preserve samples appropriately before performing subsequent molecular methods (Sandle, 2011; Rucker et al., 2018).

      Since bulky EVA suits can restrict movement and limit visibility through the helmet visor, the primary objective was aimed to evaluate the interface between a fully suited test subject handling the EVA swab tool by the crew. Fully suited testing is important for identifying tool design issues prior to flight. At exploration destinations, such as Mars, suited crew may be required to periodically sample their suits as part of an environmental monitoring protocol. In addition, a benefit of this test was an opportunity to characterize the microorganisms found on or near selected suit pressure joints under vacuum and when the spacesuits were positively pressured, enabling NASA to assess exploration mission operations and hardware design to mitigate microbial leakage.

      In this study EVA swab tools were used to collect several samples from variety of spacesuits (n = 7 sets; 48 samples) in a JSC training session. Spacesuit samples were treated with (allowing measurement of viable/intact cells) or without propidium monoazide (PMA, dead and alive cells) (48 samples each of PMA and no PMA; total n = 98 samples), a DNA intercalating dye before utilizing molecular technologies (Vaishampayan et al., 2013). The viable microbial burden targeting 16S rRNA gene (for bacteria/archaea) and internal transcribed region (ITS; for fungi) were estimated using quantitative polymerase chain reaction (qPCR) assay and shotgun metagenome sequencing (Singh et al., 2018). Furthermore, culturable microbial burden associated with spacesuits was measured using the traditional culture-based colony counts. This study will provide NASA with the ability to evaluate the spectrum of microbial diversity associated with spacesuits.

      Materials and Methods EVA Swab Material Selection

      Validation of the macrofoam swabs (EVA swab tool material) to collect microorganisms from various material types was not a part of this study. However, a comprehensive study was performed previously to understand the suitable swab materials (cotton, polyester, and macrofoam) in the efficient removal of the microorganisms from the aluminum and titanium surfaces (Kwan et al., 2011). Briefly, a model microbial community comprised 11 distinct species of bacterial, archaeal, and fungal lineages, was used to examine the effects of variables in sampling matrices, target cell density/molecule concentration, and cryogenic storage on the overall efficacy of the sampling regimen. The biomolecules and cells/spores recovered from each collection device were assessed by cultivable and microscopic enumeration, and quantitative and species-specific PCR assays. rRNA gene-based quantitative PCR analysis showed that cotton swabs were superior to nylon-flocked swabs and macrofoam swabs significantly outperformed polyester wipes. Furthermore, macrofoam swab materials were found to withstand extreme temperature fluctuations of the space conditions including varying pressure, and vacuum (Rucker et al., 2018).

      EVA Swab Sample Kit Preparation and Sample Collection

      Three different kinds of spacesuits were sampled (Figure 1). Briefly, the EMU suits are currently used for EVA on ISS, but are not designed for use in planetary missions. We sampled stainless steel wrist joints and cloth gauntlets covering the joints on these suits. The outer fabric of the EMU is made of Ortho-Fabric, which is a blend of Gortex (ePTFE), Kevlar (a para-aramid type fiber related to nylon) and Nomex (a meta-aramid type fiber) (Newman et al., 2000). The MACES and OCSS suits designed for internal cabin use, such as inside Orion during launch and reentry through Earth’s atmosphere, use similar wrist joint as the EMU but without a gauntlet to cover it. The outer layer of the MACES and OCSS suits is comprised of orange Nomex (Watson, 2014; NASA, 2019b). NASA has conducted a series of ground tests intended to evaluate the EVA swab kit’s form, fit, and function under mission operations scenarios, in preparation for eventual sample collection from outside the ISS (Rucker et al., 2018). For samples collected from the EMU, EVA swabbing was an add-on to a routine suit familiarization test that all flight crew are required to perform. Familiarization involves suit fit and functional checks, followed by a 4-h prebreathe protocol (to mitigate potential for decompression sickness) before exposure to vacuum in the Space Station Airlock Test Article chamber. Spacesuit samples were collected during the prebreathe protocol, when the crew member was breathing pure oxygen at a suit internal pressure 4.3 psi higher than ambient external pressure, but not yet at external vacuum pressure. Although standard laboratory swabs could have been used under these conditions, this test provided an opportunity for suited crew to practice self-swabbing with the flight-like EVA swab kit, which will be necessary in future studies where samples will be taken under external vacuum conditions. A second series of tests was conducted with the MACES and OCSS suits. In these tests, four test subjects sampled their own suits (two MACES suits and two OCSS suits) inside the 11-foot vacuum chamber while the chamber was at vacuum (0.01 torr). The internal pressure inside the suits was 4.3 psi. Samples collected during these tests were exposed to a maximum of 4 h of vacuum.

      Photographs showing collection of various types of samples from different spacesuit types. (A–I) Images of sample collection and associated equipment, intended for use in flight. (J) An image of an Extravehicular Mobility Unit (EMU) suit. (K) An image of an Advanced Crew Escape Suit (ACES). Total number of samples collected are tabulated. N/A: Not available.

      Sample kit cleaning, sterilization, and assembly were performed at JSC according to a purpose-developed protocol. Each sample canister (assembled with filter and ball plungers) and swab end effector assembly was placed into separate autoclave bags. Bagged components were placed into Steris LV 250 laboratory steam sterilizer and sterilized using a gravity cycle of 45 min at 121°C (250°F) and 103.4 kPa (15 psi). Note that neither the sample caddy itself nor the tool handle which was never in contact with the swab head were autoclaved. Bagged components were allowed up to 1 h of cool-down time at approximately 22°C (71.6°F) for safe handling. Following autoclaving, bagged components were transferred to a Labconco Horizontal Clean Bench (Model # 36100000, ISO Class 5). With the commercial swab inside its sterile packaging, the swab stem was cut to optimal length [approximately 6.0 cm (2.4 in)] using sterilized wire cutters, making sure the swab head remained inside its packaging until the final assembly step. The cut end of the swab was then inserted into the end effector slot, and set screws were tightened to hold the swab in place. Sterile packaging was removed from the swab head immediately before inserting each swab assembly into its sterile container. Each container/swab assembly was then mounted into the tool caddy, wiped clean with isopropanol and placed into bonded storage until the test.

      During swab assembly, technicians wore sterile gloves, and both the gloves and assembly tools (Allen wrench, scissors, and forceps) were sprayed with 70% ethanol surface disinfectant. All parts were handled either with sterile forceps or the autoclave bags, with no contact between the gloves and tool areas that must remain sterile. After assembly, the EVA sample kits were transported to the test site packed inside hard-sided storage cases. Once at the test site, the analog crew were briefed on tool usage and were given an opportunity to practice with a spare handle and sample caddy assembly. Over a period of 7 months between December 2016 and June 2017, 176 spacesuits, environmental control, and floor samples were collected during eight sampling time periods at JSC. Figure 1 shows sample collection from various parts of the spacesuits, EVA sampling kits, and number of samples associated with various spacesuits. The specific location for each sampling event of these 48 samples, surface area, and collection dates are given in Table 1 and detailed metadata about spacesuit types used, fabrics composition, microbial burden, cultivable diversity, are given in the Supplementary Tables 1, 2.

      Controls

      Among 48 spacesuit samples including five controls samples were further analyzed for various microbiological characteristics using traditional and shotgun metagenomic sequence analyses. Environmental controls were the swabs that were removed from the canister during testing but not touched to any surface. Negative controls were swabs that were not opened at all during testing. Among these 48 samples, 36 were from EMU and 12 were from MACES spacesuits.

      Sample Processing

      After sample collection, sample processing took place in an ISO 7 (10K class) cleanroom at JPL. Under ISO 5 certified biosafety cabinet, each swab was aseptically severed with a sterile cutter and transferred to a 50 mL Falcon tube containing 15 mL of sterile phosphate-buffered saline (PBS; pH 7.4). The tube with the swab was shaken for 2 min followed by concentration with a Concentrating Pipette (Innova Prep, Drexel, MO, United States) using 0.22 μm Hollow Fiber Polysulfone tips (Cat #: CC08022) and PBS as elution fluid. Each sample was concentrated to 5 mL. A 100 μL concentrated aliquots were plated on various agar plates to estimate cultivable population using traditional plate count methods (described below). One mL of the diluted solution (200 μL plus 1.8 mL PBS) was used to conduct an ATP assay (Kikkoman Corp., Noda, Japan) to rapidly measure total and viable microbial population (Venkateswaran et al., 2003), enabling appropriate serial dilutions. Furthermore, 3 mL of each concentrated sample was split into two 1.5 mL- aliquots and one aliquot was treated with PMA to assess viability (Vaishampayan et al., 2013), while the second aliquot was handled similarly but without the addition of PMA. Briefly, 18.25 μL of 2 mM PMA was added to one half of the 3-mL sample (final concentration 25 μM) followed by 5 min incubation at room temperature in the dark and 15 min exposure to the activation system (PMA LED device, Biotium, Hayward, CA, United States). Each sample was then split into two 0.75 mL aliquots. One aliquot was transferred to bead beating tubes containing Lysing Matrix E (MP Biomedicals, Santa Ana, CA, United States), followed by bead beating for 60 s using the vortex sample holder (MO Bio, Carlsbad, CA, United States). The bead-beaten aliquot and the aliquot without bead beating were combined for their corresponding PMA-treated and non-treated samples. DNA extraction was accomplished with the Maxwell 16 automated system (Promega, Madison, WI, United States), in accordance with manufacturer instructions. A Maxwell control (MC) without any sample added in its cartridge was run concurrently with each flight sample set to account for microbial contamination associated with reagents (kitome) used in the automated DNA extraction. The extracted DNA was eluted in 50 μL of water and stored at −20°C and processed with the rest of the samples later.

      Estimation and Identification of Cultivable Microbial Population

      The 100 μl of each concentrated sample were plated on Reasoner’s 2A agar (R2A for environmental microbes), Potato Dextrose Agar with chloramphenicol (100 μg/mL; PDA for fungi), and blood agar (BA for human commensals; Hardy Diagnostics, Santa Maria, CA, United States) in duplicate. R2A and PDA plates were incubated at 25°C for 7 days and BA plates at 37°C for 2 days at which time colony forming units (CFU) were counted. All colonies were picked from each plate and from each suit sampling location. The isolates were then archived in semisolid R2A or PDA slants (agar media diluted 1:10) and stored at room temperature. Once a culture was confirmed to be pure, two cryobead stocks (Copan Diagnostics, Murrieta, CA, United States) were prepared for each isolate and stored at –80°C. A loopful of purified microbial culture was directly subjected to PCR, and the targeted fragment was amplified (colony PCR), or DNA was extracted with the UltraClean DNA kit (MO Bio, Carlsbad, CA, United States) or Maxwell 16 instrument. The extracted DNA was used for PCR to amplify the 1.5 kb 16S rRNA gene to identify bacterial strains. The following primers were used for the 16S rRNA gene amplification to estimate bacterial population. The forward primer, 27F (5′-AGA GTT TGA TCC TGG CTC AG-3′) and the reverse primer, 1492R (5′-GGT TAC CTT GTT ACG ACT T-3′) (Lane, 1991; Turner et al., 1999). The PCR conditions were as follows: denaturation at 95°C for 5 min, followed by 35 cycles consisting of denaturation at 95°C for 50 s, annealing at 55°C for 50 s, and extension at 72°C for 1 min 30 s and finalized by extension at 72°C for 10 min. For fungal population estimation, the forward primer ITS 1F (5′-TTG GTC ATT TAG AGG AAG TAA-3′) (Lai et al., 2007) and reverse primer Tw13 (5′-GGT CCG TGT TTC AAG ACG-3′) (Taylor and Bruns, 1999) were used to obtain ∼1.2 kb ITS product. The PCR conditions were as follows: Initial denaturation at 95°C for 3 min followed by 25 cycles of 95°C for 50 s, annealing at 58°C for 30 s, and extension at 72°C for 2 min, followed by a final extension at 72°C for 10 min. The amplicons were inspected by gel electrophoresis in 1% agarose gel. When bands for products were visible, amplification products were treated with Antarctic phosphatase and exonuclease to remove 5′- and 3′- phosphates from unused dNTPs before sequencing. The sequencing was performed (Rockville, MD, United States) using 27F and 1492R primers for Bacteria, and ITS1F and Tw13 primers for fungi. The sequences were assembled using SeqMan Pro from DNAStar Lasergene Package (DNASTAR Inc., Madison, WI, United States). The bacterial sequences were searched against EzTaxon-e database (Kim et al., 2012) and the fungal sequences against the UNITE database (Koljalg et al., 2013). The identification was based on the closest percentage similarity (>97%) to previously identified microbial type strains.

      qPCR Assay

      Following the DNA extraction, quantitative polymerase chain reaction (qPCR), targeting the partial 16S rRNA gene (bacteria) or partial ITS region (fungi), was performed with SmartCycler (Cepheid, CA, United States) to quantify the microbial burden as previously established (Checinska Sielaff et al., 2019). Each 25-μL reaction consisted of 12.5 μL of 2X iQ SYBR Green Supermix (BioRad, Hercules, CA, United States), 1 μL each of forward and reverse oligonucleotide primers (10 μM each), and 1 μL of template DNA (PMA treated and non-treated samples). Each sample was run in triplicate; the average and standard deviation were calculated based on these results. Purified DNA from a model microbial community (Kwan et al., 2011) served as the positive control and DNase/RNase free molecular-grade distilled water (Promega, Madison, WI, United States) was used as the negative control in each run. The number of gene copies was determined from the standard curve as described previously with a modification where synthetic fragments of B. pumilus (1.4 kb 16S rRNA gene) or Aureobasidium pullulans (1-kb ITS region) were used instead of genomic DNA (Checinska et al., 2015). The qPCR efficiency was ∼98%. The negative control values were not deducted since the values were at ∼100 copies per 1 or 10 μL and not scalable (yielded the same results despite using 1 μL and 10 μL of DNA templates).

      Illumina Based DNA Sequencing and Analysis

      The initial DNA yield and metagenome library quantitation of all 96 samples tested (48 samples PMA treated and 48 samples PMA untreated) were measured by Qbit (Thermo Fisher Scientific Inc., United States). DNA libraries for all 96 samples were prepared for shotgun metagenome sequencing using the Nextera DNA Library Preparation Kit from Illumina. The quality and fragment size of each library were assessed on a Bioanalyzer 2100 (Agilent). Separate adapters were added to the DNA from each library, normalized to 2 nM, pooled, denatured, and diluted to 1.8 pM according to the standard recommendations by Illumina. The HiSeq4000 platform (Illumina) was used for sequencing, resulting in 100-bp paired-end reads.

      Bioinformatics Analysis

      Bioinformatic analyses were performed on Weill Cornell Medicine’s Athena compute cluster, a typical high-performance grid compute (Slurm) system. The secondary analysis was performed on Linux and MacOS systems. Unless otherwise noted programs were run with default settings.

      Data Quality Control and Filtering

      Sequence data were processed with AdapterRemoval (v2.17) to remove low-quality reads and reads with ambiguous bases. Subsequently, reads were aligned to the human genome (hg38, including alternate contigs) using Bowtie2 (v2.3.0, fast preset). Read pairs where both ends mapped to the human genome were separated from read pairs where neither mate mapped. Read pairs where only one mate mapped were discarded. Hereafter, we refer to these read sets as human reads and non-human reads. We did not process human reads beyond counting the total fraction of DNA from our samples which mapped to the human genome.

      Taxonomic Profiling and Analysis

      Taxonomic profiles were generated by processing non-human reads KrakenUniq (v0.3.2) with a database based on all draft and reference genomes in RefSeq Microbial (bacteria, fungi, virus, and archaea) ca. March 2017. KrakenUniq uses a k-mer based approach to identify reads. Reads are broken into k-mers of 31 bases. Each k-mer is mapped to a database that lists the lowest common ancestor of all genomes which contained the k-mer. Each read is classified by identifying the best supported path in the taxonomic tree of markers. Finally, the taxonomic makeup of a sample is given by concatenating annotations for reads without further processing. KrakenUniq counts the number of unique marker k-mers assigned to each taxa, and we filtered taxa with fewer than 512 unique markers. Differential abundance estimation (where applicable) using the ALDEx2 R package was performed (Fernandes et al., 2013). Briefly, ALDEx2 transforms read count matrices using a centered log ratio transformation that models samples as Dirichlet-Multinomial distributions over taxa then compares taxonomic abundances across groups. If two groups are given, comparison is done with a Wilcoxon rank sum test, more than two groups are tested via a generalized linear model. All p-values are multiple hypotheses corrected using Benjamini-Hochberg. We considered a taxon to have differential abundance in a given condition if its corrected p-value was less than or equal to p = 0.05.

      Dimensionality reduction of taxonomic profiles was performed with Uniform Manifold Approximation and Projection UMAP (McInnes et al., 2018) based on a matrix of Jensen-Shannon Divergences (JSD) between samples. Analysis of intersample diversity (beta-diversity) was achieved using the same matrix of JSD. Intrasample diversity (alpha-diversity) was measured by taking Shannon’s Entropy of the total sum normalized taxonomic profile of each sample. Rarefaction analysis of taxa was performed by grouping samples by location and setting and selecting 16 uniform random groups for each value. A curve of best fit was found by fitting a logarithmic model to the series.

      Profile of Eukaryotic species were generated using CLARK-S (v1.2.5) (Ounit and Lonardi, 2016) using sequences classified with high confidence (i.e., confidence score > 0.75, and gamma score > 0.10) as defined in the CLARK manual. Identification of taxa was further restricted to species with relative abundance at least 0.01% of the total sequences.

      Samples were compared to eight representative samples of human body sites selected from the Human Microbiome Project (HMP) (Turnbaugh et al., 2007) for each of five body sites: oral, skin, airways, gastrointestinal, and urogenital. Using MetaPhlAn2 (v2.2) (Truong et al., 2015), we generated taxonomic profiles for HMP samples and our samples and compared profiles using Cosine Similarity.

      Functional Profiling and Analysis

      HUMANn2 (Franzosa et al., 2018) was used to generate functional metabolic profiles of the genes in our samples. Non-human reads were aligned to Uniref90 (ca. March 2017) using the DIAMOND aligner (v0.8.6) (Buchfink et al., 2015). Subsequently, alignments were processed using HUMANn2 (v0.11.1) to produce profiles of pathway abundance. Pathways were tested for differential expression using the Wilcoxon rank sum corrected by Benjamini Hochberg. Dimensionality reduction of pathways was performed using PCoA over normalized pathway abundances.

      Profiling Antimicrobial Resistance Genes

      Profiles of antimicrobial resistance (AMR) genes using MegaRes (v1.0.1) (Lakin et al., 2017) were created. To generate profiles from MegaRes, we mapped non-human reads to the database using Bowtie2 (v2.3.0, very-sensitive presets) (Langmead and Salzberg, 2012). Subsequently, alignments were analyzed using ResistomeAnalyzer (commit 15a52dd) (Dean, 2018) and normalized by total reads per sample and gene length to give Reads per kilo base per million mapped reads (RPKMs). MegaRes includes an ontology grouping resistance genes into, gene classes, AMR mechanisms, and gene groups.

      Identification of Genomes and Strains

      We assembled contigs from all PMA treated samples using MegaHIT (v1.1.3) (Li et al., 2015) then clustered the resulting contigs into draft genomes using MetaBAT2 (Kang et al., 2019). Draft genomes were quality controlled and assigned a rough taxonomic rank using CheckM (Parks et al., 2015). Genomes with less than 50% completeness or more than 20% contamination were discarded. We aligned all genomes to one another to using Nucleotide MUMmer (Delcher et al., 2003) and processed the results to generate an Average Nucleotide Identity (ANI) between all pairs of draft genomes. We discarded all alignments that covered less than half the average lengths of the genomes being aligned. We further discarded alignments with less than 99% ANI so that we would only be left with pairs of nearly identical genomes. We grouped these alignments into connected components and analyzed the sites where each component was found.

      Results Microbial Abundance

      A total of 48 samples (36 EMU and 12 MACES) were collected from six different surfaces of the spacesuits or environmental controls. Sampling surfaces include: left wrist joint (12 samples), left inner glove gauntlet (5 samples), left outer glove gauntlet (5 samples), right wrist joint (11 samples), right inner glove gauntlet (4 samples), right outer glove gauntlet (6 samples). All controls were analyzed for all microbiological and molecular biological examinations (5 samples, Table 1). All these 48 samples were categorized into sets (n = 7 sets) based on the suit types or sample collection dates (Table 1). In addition, metadata such as locations, type of suits, materials of spacesuits, and date of collection are given in Supplementary Table 1.

      Characteristics of various spacesuites sampled during this study and associated metadata.

      Set # Sampled locations Metagenome sample ID Sampling date Suit types sampled Material type sampled Sample pressure PMA or no PMA
      SET-1 Control Swab – Not removed from canister JC-044 12/28/2016 EMU N/A 760 No PMA
      SET-1 Control Swab – Not removed from canister JC-092 12/28/2016 EMU N/A 760 PMA
      SET-1 Outside, left wrist gauntlet JC-018 12/28/2016 EMU Beta cloth 760 No PMA
      SET-1 Outside, left wrist gauntlet JC-066 12/28/2016 EMU Beta cloth 760 PMA
      SET-1 Inside, left wrist gauntlet JC-013 12/28/2016 EMU Beta cloth 760 No PMA
      SET-1 Inside, left wrist gauntlet JC-061 12/28/2016 EMU Beta cloth 760 PMA
      SET-1 Left wrist joint groove JC-001 12/28/2016 EMU Stainless Steel 760 No PMA
      SET-1 Left wrist joint groove JC-049 12/28/2016 EMU Stainless Steel 760 PMA
      SET-1 Outside, right wrist gauntlet JC-038 12/28/2016 EMU Beta cloth 760 No PMA
      SET-1 Outside, right wrist gauntlet JC-086 12/28/2016 EMU Beta cloth 760 PMA
      SET-1 Inside, right wrist gauntlet JC-034 12/28/2016 EMU Beta cloth 760 No PMA
      SET-1 Inside, right wrist gauntlet JC-082 12/28/2016 EMU Beta cloth 760 PMA
      SET-1 Right wrist joint groove JC-023 12/28/2016 EMU Stainless Steel 760 No PMA
      SET-1 Right wrist joint groove JC-071 12/28/2016 EMU Stainless Steel 760 PMA
      SET-2 Control Swab – Not removed from canister JC-045 12/14/2016 EMU N/A 760 No PMA
      SET-2 Control Swab – Not removed from canister JC-093 12/14/2016 EMU N/A 760 PMA
      SET-2 Outside, left wrist gauntlet JC-019 12/14/2016 EMU Beta cloth 760 No PMA
      SET-2 Outside, left wrist gauntlet JC-067 12/14/2016 EMU Beta cloth 760 PMA
      SET-2 Left wrist joint groove JC-002 12/14/2016 EMU Stainless Steel 760 No PMA
      SET-2 Left wrist joint groove JC-050 12/14/2016 EMU Stainless Steel 760 PMA
      SET-2 Outside, Right wrist gauntlet JC-039 12/14/2016 EMU Beta cloth 760 No PMA
      SET-2 Outside, Right wrist gauntlet JC-087 12/14/2016 EMU Beta cloth 760 PMA
      SET-2 Right wrist joint groove JC-024 12/14/2016 EMU Stainless Steel 760 No PMA
      SET-2 Right wrist joint groove JC-072 12/14/2016 EMU Stainless Steel 760 PMA
      SET-4 Long Term Control assembled 2/6 tested 3/16 JC-047 3/15/2017 EMU N/A 760 No PMA
      SET-4 Long Term Control assembled 2/6 tested 3/16 JC-095 3/15/2017 EMU N/A 760 PMA
      SET-3 Exterior, palm-side left wrist gauntlet JC-020 2/6/2017 EMU Beta cloth 760 No PMA
      SET-3 Exterior, palm-side left wrist gauntlet JC-068 2/6/2017 EMU Beta cloth 760 PMA
      SET-3 Left wrist joint groove JC-003 2/6/2017 EMU Stainless Steel 760 No PMA
      SET-3 Left wrist joint groove JC-051 2/6/2017 EMU Stainless Steel 760 PMA
      SET-3 Interior, left wrist gauntlet JC-014 2/6/2017 EMU Beta cloth 760 No PMA
      SET-3 Interior, left wrist gauntlet JC-062 2/6/2017 EMU Beta cloth 760 PMA
      SET-3 Right wrist joint groove JC-025 2/6/2017 EMU Stainless Steel 760 No PMA
      SET-3 Right wrist joint groove JC-073 2/6/2017 EMU Stainless Steel 760 PMA
      SET-3 Control Swab – Not removed from canister JC-046 2/6/2017 EMU N/A 760 No PMA
      SET-3 Control Swab – Not removed from canister JC-094 2/6/2017 EMU N/A 760 PMA
      SET-4 Left wrist outer gauntlet JC-021 3/15/2017 EMU Beta cloth 760 No PMA
      SET-4 Left wrist outer gauntlet JC-069 3/15/2017 EMU Beta cloth 760 PMA
      SET-4 Left wrist inner gauntlet JC-015 3/15/2017 EMU Beta cloth 760 No PMA
      SET-4 Left wrist inner gauntlet JC-063 3/15/2017 EMU Beta cloth 760 PMA
      SET-4 Left glove/lower arm groove JC-004 3/15/2017 EMU Stainless Steel 760 No PMA
      SET-4 Left glove/lower arm groove JC-052 3/15/2017 EMU Stainless Steel 760 PMA
      SET-4 Right wrist outer gauntlet JC-040 3/15/2017 EMU Beta cloth 760 No PMA
      SET-4 Right wrist outer gauntlet JC-088 3/15/2017 EMU Beta cloth 760 PMA
      SET-4 Right wrist inner gauntlet JC-035 3/15/2017 EMU Beta cloth 760 No PMA
      SET-4 Right wrist inner gauntlet JC-083 3/15/2017 EMU Beta cloth 760 PMA
      SET-4 Right glove/lower arm groove JC-026 3/15/2017 EMU Stainless Steel 760 No PMA
      SET-4 Right glove/lower arm groove JC-074 3/15/2017 EMU Stainless Steel 760 PMA
      SET-5 Control Swab – Not removed from canister JC-048 3/16/2017 MACES/OCCS N/A 0.01 No PMA
      SET-5 Control Swab – Not removed from canister JC-096 3/16/2017 MACES/OCCS N/A 0.01 PMA
      SET-5 Left wrist JC-005 5/16/2017 OCCS Stainless Steel 0.01 No PMA
      SET-5 Left wrist JC-053 5/16/2017 OCCS Stainless Steel 0.01 PMA
      SET-5 Right wrist JC-027 5/16/2017 OCCS Stainless Steel 0.01 No PMA
      SET-5 Right wrist JC-075 5/16/2017 OCCS Stainless Steel 0.01 PMA
      SET-5 Left wrist JC-006 5/30/2017 MACES Stainless Steel 0.01 No PMA
      SET-5 Left wrist JC-054 5/30/2017 MACES Stainless Steel 0.01 PMA
      SET-5 Right wrist JC-028 5/30/2017 MACES Stainless Steel 0.01 No PMA
      SET-5 Right wrist JC-076 5/30/2017 MACES Stainless Steel 0.01 PMA
      SET-6 Left wrist outer gauntlet JC-022 5/30/2017 EMU Beta cloth 760 No PMA
      SET-6 Left wrist outer gauntlet JC-070 5/30/2017 EMU Beta cloth 760 PMA
      SET-6 Left wrist inner gauntlet JC-016 5/30/2017 EMU Beta cloth 760 No PMA
      SET-6 Left wrist inner gauntlet JC-064 5/30/2017 EMU Beta cloth 760 PMA
      SET-6 Left glove/lower arm groove JC-007 5/30/2017 EMU Stainless Steel 760 No PMA
      SET-6 Left glove/lower arm groove JC-055 5/30/2017 EMU Stainless Steel 760 PMA
      SET-6 Right wrist outer gauntlet JC-041 5/30/2017 EMU Beta cloth 760 No PMA
      SET-6 Right wrist outer gauntlet JC-089 5/30/2017 EMU Beta cloth 760 PMA
      SET-6 Right wrist inner gauntlet JC-036 5/30/2017 EMU Beta cloth 760 No PMA
      SET-6 Right wrist inner gauntlet JC-084 5/30/2017 EMU Beta cloth 760 PMA
      SET-6 Right glove/lower arm groove JC-029 5/30/2017 EMU Stainless Steel 760 No PMA
      SET-6 Right glove/lower arm groove JC-077 5/30/2017 EMU Stainless Steel 760 PMA
      SET-6 Right wrist outer gauntlet JC-042 6/14/2017 EMU Beta cloth 760 No PMA
      SET-6 Right wrist outer gauntlet JC-090 6/14/2017 EMU Beta cloth 760 PMA
      SET-6 Left Wrist Crew 3 JC-008 6/8/2017 MACES Stainless Steel 0.01 No PMA
      SET-6 Left Wrist Crew 3 JC-056 6/8/2017 MACES Stainless Steel 0.01 PMA
      SET-6 Right Wrist Crew 3 JC-030 6/8/2017 MACES Stainless Steel 0.01 No PMA
      SET-6 Right Wrist Crew 3 JC-078 6/8/2017 MACES Stainless Steel 0.01 PMA
      SET-6 L Wrist Crew 4 JC-009 6/8/2017 MACES Stainless Steel 0.01 No PMA
      SET-6 L Wrist Crew 4 JC-057 6/8/2017 MACES Stainless Steel 0.01 PMA
      SET-6 Left wrist crew 3 JC-010 6/12/2017 OCCS Stainless Steel 0.01 No PMA
      SET-6 Left wrist crew 3 JC-058 6/12/2017 OCCS Stainless Steel 0.01 PMA
      SET-6 Right wrist crew 3 JC-031 6/12/2017 OCCS Stainless Steel 0.01 No PMA
      SET-6 Right wrist crew 3 JC-079 6/12/2017 OCCS Stainless Steel 0.01 PMA
      SET-6 Left wrist crew 4 JC-011 6/12/2017 MACES Stainless Steel 0.01 No PMA
      SET-6 Left wrist crew 4 JC-059 6/12/2017 MACES Stainless Steel 0.01 PMA
      SET-6 Right wrist crew 4 JC-032 6/12/2017 MACES Stainless Steel 0.01 No PMA
      SET-6 Right wrist crew 4 JC-080 6/12/2017 MACES Stainless Steel 0.01 PMA
      SET-7 Right wrist outer gauntlet JC-043 6/26/2017 EMU Beta cloth 760 No PMA
      SET-7 Right wrist outer gauntlet JC-091 6/26/2017 EMU Beta cloth 760 PMA
      SET-7 Right wrist inner gauntlet JC-037 6/26/2017 EMU Beta cloth 760 No PMA
      SET-7 Right wrist inner gauntlet JC-085 6/26/2017 EMU Beta cloth 760 PMA
      SET-7 Right glove/lower arm groove JC-033 6/26/2017 EMU Stainless Steel 760 No PMA
      SET-7 Right glove/lower arm groove JC-081 6/26/2017 EMU Stainless Steel 760 PMA
      SET-7 Left wrist inner gauntlet JC-017 6/26/2017 EMU Beta cloth 760 No PMA
      SET-7 Left wrist inner gauntlet JC-065 6/26/2017 EMU Beta cloth 760 PMA
      SET-7 Left glove/lower arm groove JC-012 6/26/2017 EMU Stainless Steel 760 No PMA
      SET-7 Left glove/lower arm groove JC-060 6/26/2017 EMU Stainless Steel 760 PMA
      N/A, Not applicable; PMA, Propidium monoazide treated to measure viable microorganisms and no-PMA are the samples constitute both dead and live microbes; EMU, Extravehicular Mobility Unit; MACES, Modified Advanced Crew Escape System; OCCS, Onboard Complex Control System.

      Our samples contained viable bacterial populations which were estimated by culture-dependent and independent analyses and are summarized in Supplementary Table 1 and Figure 2A. We cultured various microorganisms from our samples on three different types of media: blood agar, R2A, and PDA. The number of cultivable bacterial counts on R2A plates ranged from no growth to 9.0 × 102 CFU per 25 cm2 (Figure 2A). The bacterial counts on blood agar were ranged from no growth to 3.5 × 102 CFU per 25 cm2. No bacterial colonies were observed in any of the controls during this study. The phylogenetic affiliation of the bacterial strains isolated in this study was shown in Supplementary Figure 1A. Among 24 bacterial strains isolated and identified, the microorganisms belonged to the members of the phyla Firmicutes (13 strains), Actinobacteria (10 strains) and Proteobacteria (1 strain). Bacillus species represented the highest number of isolates, followed by Arthrobacter species. Comparatively, fungal isolates were not abundant and only six strains belonging to six different species were isolated. The ITS-based sequence analyses identified them as Epicorum nigrum, Alternaria sp., Penicillium fagi, Aureobasidium pullulans, Naganishia adeliensis, and Neonectria sp. The results of ATP-assay were not shown but ATP contents were used to estimate the microbial burden which further helped to determine appropriate serial dilutions.

      (A) Number of CFU found when plating samples on three different media. R2A, Reasoner’s 2A agar (for environmental bacteria); PDA, potato dextrose agar with chloramphenicol (100 μg/mL; for fungi); BA, blood agar (for human commensals). (B) Number of copies of the 16S rRNA gene found on different suit component: inner gauntlet (IG), outer gauntlet (OG), and wrist joint (WJ). (C) Number of copies of the ITS gene found on different components. (D) Total number of species detected by metagenomics on different components. All refers to microbes which were found in all components from a given set while Any refers to species on at least one component in a set. (E) Rarefaction curves showing the diversity of different numbers of components pooled together.

      A qPCR assay of the 16S and ITS genes was performed to measure the absolute microbial population of both viable (PMA treated) and total (PMA-untreated) microorganisms. This assay did not show a statistically significant difference in the microbial load among various locations sampled on spacesuits tested nor in various sets categorized. Viable bacterial load (PMA treated samples) was estimated at approximately 105 to 106 16S rRNA copies per 25 cm2, Figure 2B. Viable bacterial population was an order of magnitude less abundant than total bacterial burden that include both dead and live microorganisms (Supplementary Table 1). Viable fungal population was measured at approximately 102 to 104 ITS copies per 25 cm2, Figure 2C. No significant difference was observed between EMU and MACES suits in either cultivable and culture-independent microbial burden assays.

      Molecular Microbial Diversity

      Spacesuits are modular, each set refers to a single assembled set of components operated and sampled on a given day. The 48 samples including five controls were either treated with PMA or left untreated, resulting in an analysis of 96 samples. Among the 96 samples subjected for shotgun library preparation, all samples yielded enough DNA fragments except four PMA-treated samples and one non-PMA treated sample, hence 91 samples were subsequently assayed for shotgun metagenome sequencing. The PMA treated samples that did not yield any shotgun metagenome reads were SET-2 outside, left wrist gauntlet; SET-3 interior, left wrist gauntlet; SET-3 right wrist joint groove; and SET-7 left glove/lower arm groove.

      In total, 319M reads were generated from all 91 samples. Human (∼38.2%) and animal (∼30%) associated reads were removed from the analyses. The PMA (49.8M) and non-PMA (54.7M) reads were ∼30% of the total reads. Approximately 104.5M reads associated with microorganisms were generated after high quality trimming from PMA (44 samples) and non-PMA treated (47 samples) samples. Dimensionality reduction comparing microbial taxonomic profiles between PMA treated and untreated samples showed an average shift based on PMA treatment (Supplementary Figure 2) suggesting that some types of microbes may be present on spacesuits as non-viable detritus. PMA treated samples were the focus of this study as they represent the intact/viable cells and information about PMA untreated samples were presented in supplementary datasets (Supplementary Tables 2, 3). The PMA-based analyses revealed that there were no microbial diversity differences among the EMU and MACES suits.

      For all PMA treated samples, at domain level, the majority of the reads were assigned to bacteria (98.6%), followed by eukaryotes (0.9%), then archaea (0.24%), and viral signatures were 0.17%. For samples not treated with PMA, these reads were assigned to bacteria (98.6%), followed by eukaryotes (0.9%), archaea (0.5%), and viruses (0.1%). The proportional abundance of bacteria and fungi were similar in both PMA treated and non-PMA treated samples. When the relative abundance of all metagenomics reads was summed, ∼80% of the reads were attributed to the species whose reads were >100K.

      None of the control samples yielded microbes that could be cultured in the media employed during this study which confirms that the EVA tool kit prepared for this study was sterile. But when all samples were considered for molecular analyses, ∼5% of the total metagenomics reads associated with bacteria, fungi, and viruses were present in control samples (n = 5). Among 993 microbial species observed in all spacesuits including control during this study (Supplementary Table 2), 13 bacterial taxa of control samples exhibited >10K reads and they were identified as Bacillus pumilus, Cutibacterium acnes, Janthinobacterium species (n = 3), Micrococcus luteus, Negativicoccus massiliensis, Pseudomonas species (n = 5), and Ralstonia insidiosa. Among them, C. acnes, Janthinobacterium species, Pseudomonas species, and R. insidiosa members were present in all five control samples. The bacterial species associated with controls that exhibited >100K reads were N. massiliensis (512K reads), C. acnes (448K reads), and Pseudomonas sp. NC02 (347K reads). Hence, few contaminant species were found as “kitomes” during this study and our finding is based on identifying microbial species/strains that are not in controls.

      When various sets of spacesuits were compared, some differences were observed. Set #7 samples consist of members of the genera Methylobacterium and Curtobacterium whereas Pseudomonas species were prevalent in samples collected from set #5. Among 350 bacterial genera constituting 660 bacterial taxa identified, sequences of the members of the genera Curtobacterium and Methylobacterium were retrieved across all sets of spacesuits in high abundance. The compositional analysis showed a higher abundance of Curtobacterium, Methylobacterium, Negativicoccus, and Pseudomonas that exhibited more than two million reads. Among bacterial species identified (60 species > 100K reads; 239 species > 10K reads), higher abundance (>2M reads) of Curtobacterium acnes (8.9M reads), Methylobacterium oryzae (4.4M reads), and M. phyllosphaerae (4.2M reads) sequences were observed. Low fungal, archaeal, and viral reads were retrieved during this study and their sequence abundances and taxa characteristics are presented in Supplementary Tables 2–4.

      Molecular Microbial Diversity Indices

      The total number of microbial species (species richness) found on each type of component (Inner Glove Gauntlet, Outer Glove Gauntlet, and Wrist Joint) was similar and typically between 100 and 200 (Figure 2C). A subset of these species could be found on all components in a set (typically 50–100 species found in all three components of either the left or right side of the suit) establishing a shared community. The inner and outer suit gauntlet had higher richness than the wrist joint (p < 2–16, one-way ANOVA).

      To establish the total number of microbial species in the entire study (Figure 2D), a rarefaction analysis was performed (Figure 2E). Suits were considered as a whole and separately by component. A total of 660 microbial species were observed across all samples but a curve fit to the subsamples did not flatten which suggests that more microbial diversity would be seen with more samples collected. However, an analogous curve fit to subsets of species that occurred in all part in set did flatten, suggesting there may be a core community of 100–200 organisms common to spacesuits. Individual component types necessarily had more species than were found in all parts in set but fewer than were found in any part of a set.

      To address the study design of collecting multiple samples from the same suit, we conducted a nested analysis using a regression Generalized Linear Mixed Model, and found that alpha diversity (Shannon Index) varied significantly across spacesuits for the PMA untreated group (F5,35 = 4.84, P = 0.002) but did not vary significantly for the PMA treated group. This may be due to the higher power demands of nested models and the limited number of samples collected.

      Taxonomic Analysis of Spacesuits

      Microbial taxa were categorized based on a number of different conditions (1) differential abundance between PMA treated samples and untreated samples (determined by ALDEx2), (2) high prevalence taxa found in 31 out of 32 PMA treated EMU samples (excluding controls), (3) increased abundance in MACES suits compared to EMU suits for PMA treated samples (determined by ALDEx2), (4) differential abundance between suit components (wrist, inner, and outer gauntlets) in EMU samples treated with PMA (determined by ALDEx2). Differential abundance was defined as a Benjamini–Hochberg corrected q-value of 0.05 or less based on ALDEx2. Among PMA treated EMU samples one species Corynebacterium kroppenstedtii was identified as being significantly (q = 0.031) less abundant in wrist joint samples compared to other microbial species. Ninety-nine species were identified as differentially abundant in samples treated with PMA and untreated samples.

      UMAP plot on the taxonomic profiles of samples (Supplementary Figure 2) and PMA treated samples only are depicted in Figure 3A. As expected, PMA treated sample clearly separated from untreated samples. This shows that there is a distinct likely viable set of microbes present on the sampled spacesuits. Within PMA treated samples, generally samples from the same suit clustered together with Sets 1 and 7 as notably tight clusters. Set 7 (all EMU suit sampled on June 26, 2017) was a definite outlier in relative abundance matching the pattern observed for alpha diversity.

      (A) UMAP of taxonomic profiles from PMA treated samples from EMU suits. Color indicates the set a sample came from. (B,C) Distance between different types of suit components [inner gauntlet (IG), outer gauntlet (OG), wrist joint (WJ), and controls (cntrl)] and samples from different sets (B) or the same set (C). (D) The similarity of taxa from different components to representative samples from the Human Microbiome Project.

      Beta Diversity and Sample Differentiation

      The distance between taxonomic profiles of PMA treated samples of EMU suites was compared using JSD analysis. Dimensionality reduction of these distances using UMAP showed limited clustering by suit (Figure 3A). EMU suits were subdivided into eight triplets of samples that contained precisely one wrist, inner gauntlet, and outer gauntlet from the same suit. These triplets were in physical proximity to one another when sampled. We then compared two distributions: the distribution of distances between components in the same set and the distribution of distances between components in different sets (Figures 3B,C). The average JSD between components in the same set was 0.355 compared to 0.542 between components in different sets. A two-sided Welch’s t-test showed that these distributions did not share the same mean with p-value less than 2.0–16.

      We also compared these distributions to distance distributions for control samples. The mean JSD between suit components and control samples collected at the same time was 0.365 while the mean distance between control samples and suit components from other sets was 0.470. The distance between components from different sets was larger than the difference between controls and components from other sets based on a two-sided Welch’s t-test with p of 1.01–4. Analogously, the distances between components from the same set were less than controls with a p of 2.17–7.

      Taxonomic profiles of PMA treated samples during this study were compared to exemplar samples from the HMP. Spacesuit samples were found to be most similar to HMP skin and airway samples, suggesting that spacesuit microbiomes could originate from human skin or airway communities (Figure 3D). Notably the similarity to human body sites was not found to significantly vary by suit component or by which suit was being tested (one-way ANOVA), suggesting all components of all suits are exposed to human skin and airways microbiomes.

      Antimicrobial Resistance Genes

      Sequences of EMU PMA samples were mapped to known AMR genes and performed a rarefaction analysis of potential AMR genes (Figure 4A). Suits were considered as a whole and separately by component. Left and right gauntlet samples from the same component of the same suit were grouped together. A total observed richness of 40 AMR genes was noticed, but a curve fit to subsamples did not flatten, suggesting more diversity of AMR genes could be found. Samples from the outer gauntlet had more abundance than samples taken from the wrist, which in turn showed more AMR genes from the inner gauntlet. We grouped identified AMR genes by resistance class according to the MegaRes ontology. Five samples contained genes from the macrolide, lincosamide and streptogramin (MLS) class, 3 from the elfamycins, and just one sample contained a resistance gene from the beta-lactams (Figure 4B).

      (A) Rarefaction plot of the number of AMR genes in different suit subsets. (B) Plot of the number of samples with AMR genes from a given resistance class. MLS stands for macrolide-lincosamide-streptogramin. (C) Relative abundance of major taxa in parts per million (PPM) across PMA treated samples from EMU suits.

      Identified Microbial Species

      Among the viable microbial species, Pseudomonas species were abundant in spacesuits studied. A core microbiome (occurring in 90% of samples or more) of 40 species in EMU samples treated with PMA was determined with several species in Staphylococcus, Streptococcus, Pseudomonas, and Burkholderiales. The distribution of abundances for microbial species with the highest median relative abundances was identified. Cutibacterium acnes was the most abundant taxa followed by several Pseudomonas, Janthinobacterium, and Ralstonia insidiosa (Figure 4C).

      Fungal species (identified using CLARK-S, see methods) were also prevalent with 13 species identified in two or more samples. These include Malassezia restricta (found in all samples) which was associated with the skin microbiome of astronauts after their missions on the ISS by Sugita et al. (2016) and a number of other human commensal species. The full list of fungal species identified is given in Supplementary Table 3.

      Strain Specific Metagenome Assembled Genomes

      We built Metagenome Assembled Genomes (MAGs) from assemblies of all PMA treated samples, including controls. We identified MAGs that were found in more than one sample (99% ANI, see methods). These MAGs corresponded to two groups. One group of draft genomes was found in seven samples and was roughly categorized as a Propionibacterium species, the other group was found in five samples and was categorized as a member of Rhizobiales (Figures 5A,B). Both genome groups were fully connected, each draft genome from each sample had 99% ANI to each other sample in the group.

      Both samples where identical taxa were assembled. Names are structured as Left (L) or Right (R), the suit component, and the set number. Edges indicate identical assemblies based on 99.5% ANI. (A) Samples where a Propionibacterium taxa was found. (B) Samples where a Rhizobiales taxa was found.

      Both genomes were found in multiple samples from the same spacesuit. The Propionibacterium group was found in three samples from Set 6: the right wrist, and the inner and outer left gauntlets (Figure 5A). The Rhizobiales group was found in five samples from Set 7: the right wrist, inner and outer gauntlets and the left wrist and inner gauntlet (Figure 5B). Neither genome was found in any control sample. Since the samples where the genomes were found were treated with PMA, we further conclude that these microorganisms were likely viable.

      The incidence of multiple examples of the same genome at wrist, inner, and outer gauntlets of two suits is consistent with the possibility that viable bacteria might have migrated from the inside to the outside of a spacesuit, but could also have been deposited in all three locations during the suit donning process. For crew health, surfaces inside the space suits were cleaned using stericide wipes; however exterior surfaces (including the surfaces that were sampled), were not cleaned. The stericide wipes chemical components were: 1.5% n-alkyl dimethyl benzyl ammonium chloride (60% C14, 30% C16, 5% C12, 5% C18) and 1.5% n-alkyl ethylbenzyl ammonium chloride (50% C12, 30% C14, 17% C16, 5% C18) as maintenance procedure. Since neither MAG was found in any control it is concluded that presence of these genomes is not due to contamination but might be due to migration from one location to another location of the spacesuit.

      Discussion

      In this study, we established that viable microbes (and their MAGs) are present on the wrist assemblages of spacesuits, that certain microbial strains can survive on all three components of the wrist assembly without being found in corresponding controls, and that these microbes broadly resemble those of human commensal skin microbiomes.

      Microorganisms associated with fabrics have not historically been studied in detail (Buschle-Diller et al., 1994; Cappitelli and Sorlini, 2008; Linacre et al., 2010; Daly et al., 2012; Lee et al., 2016), whereas microbiome of human (Turnbaugh et al., 2007; Nelson et al., 2010; Jensen, 2013; Shafquat et al., 2014) and built environments (Danko et al., 2021a), including closed habitat of ISS environment (Singh et al., 2018; Checinska Sielaff et al., 2019), has received much attention (Westwood et al., 2014; Kettleson et al., 2015; Chase et al., 2016; Lax et al., 2017). One of the objectives of this study was determining if a human within a spacesuit could act as a source for the unintentional microbial contamination and pass on microbial signatures out of the spacesuits. Future in-depth sampling and testing of various spacesuit components are necessary to conclude the transmission route. A detailed microbiome analyses of spacesuit before crew occupation should be carried out and such samples are not available for this study.

      Microorganisms might escape through the clothing fibers via adherence, growth, and damage to the fibers. When synthetic fabrics were evaluated for microbial composition, micrococci were predominantly isolated both via culture and molecular methods (Callewaert et al., 2014). Prior studies have also revealed abundance of staphylococci on both cotton and synthetic fabrics, but corynebacteria were not enriched on any textile types, indicating that clothing fiber composition might promote differential growth of distinct microbes. Similarly, the spacesuits examined during the study revealed isolation of actinobacteria and staphylococci, but members of corynebacteria were not isolated using traditional methods (Supplementary Figure 1A). However, the culture-independent analyses showed presence of corynebacteria (Figure 4C).

      Characterization of microorganisms degrading the synthetic polymers such as polyvinyl chloride (PVC), polyurethane, nylon, and acrylics and their mode of action have been reviewed (Cappitelli and Sorlini, 2008). As shown in this study black fungi were isolated (Supplementary Figure 1B), however colonization of PVCs by these fungi and their degrading capabilities of plasticizers should be assessed before their concluding potential polymer degradation (Roberts and Davidson, 1986; Webb et al., 2000). The sequences associated with Candida albicans, a common skin microbiota, were not retrieved during this study, but retrieval of sequence from taxonomically related Candida dubliniensis was found from the majority of suit samples and not from controls.

      The composition of the fabrics and particles associated with them might determine the interaction of the microbes and fiber, but such phenomenon is not always uniform and large discrepancies exist. Fabric materials that are plant-based (e.g., cotton) might provide nutrients for microorganisms to degrade and also cotton fabrics were reported to adsorb sweat components thus promotes growth of microorganisms (Szostak-Kotowa, 2004). Several enzymes were reported to degrade fibers like cellulose and notably fungi secreting cellulolytic enzymes (Buschle-Diller et al., 1994). Even though synthetic fibers last longer than natural fabrics like cotton, they collect moisture between the fibers and become nutrients source for microorganisms (Szostak-Kotowa, 2004). However, during this study, no differences in microbial burden was noticed when EMU and MACES/OCCS suits were compared.

      When all PMA-treated samples were pooled into various locations of the spacesuit such as wrist joint (n = 23), inner gauntlet (n = 10), and outer gauntlet (n = 10), opportunistic pathogens were found to be associated with wrist joint rather than the inner or outer gauntlet samples (Table 2). Notably, C. kroppenstedtii, an opportunistic pathogen, has a high relative abundance in wrist joints (∼10K reads) compared to inner gauntlets (∼1K reads) and outer gauntlets (630 reads). Rothia dentocariosa, an oral cavity microbe, was found more in outer gauntlet. Similarly, the microbiome of outer and inner gauntlets appeared to harbor microorganisms associated with soil (Pseudomonas stutzeri) as well as the radiation-resistant Methylobacterium radiotolerans. The possibility of microbes surviving harsh conditions associated with space by migrating on or through spacesuits should be explored with controlled experiment before drawing the movement of microbes from one location to another location.

      Potential transmission of viable microorganisms among various locations of space suits.

      Taxa Number of PMA reads* retrieved from:
      Wrist joint (n=23) Inner gauntlet (n=10) Outer gauntlet (n=10)
      Acinetobacter baumannii 1,364 1,088
      Enterobacter cloacae 402 736
      Corynebacterium kroppenstedtii 9,987 1,102 630
      Rothia dentocariosa 1,132 804 6,972
      Kocuria rhizophila 1,022 7,258 534
      Pseudomonas stutzeri 1,917 763 214,233
      Methylobacterium radiotolerans 102,052 214,816 110,274
      Methylobacterium extorquens 12,163 26,274 12,994
      Methylobacterium nodulans 1,245 2,710 1,335
      Methylobacterium populi 2,274 4,842 2,413
      Methylobacterium sp 4-46 2,221 4,709 2,306
      Streptococcus thermophilus 2,200 1,027 625
      *The shotgun metagenome reads generated after PMA treatment andretrieved from various locations were pooled to generate potential transmission of the viable microorganisms from one location to another.

      Detecting microbes that were reported to be potentially harmful to astronaut health might be a concern. Members of Methylobacterium (12 species) dominated microbial communities on set #7 suits in this study, were reported to be opportunistic pathogens and might cause infections to immunocompromised patients (Kovaleva et al., 2014). Moreover, since astronauts are in close contact with suits while in use and shared suits present a hypothetical route for pathogen transmission, these measures can help inform potential risk. Though our work was limited to the exterior of suits, this study documented that spacesuit could harbor viable microbes. None of the microbes discovered are thought to present a health risk, but some belong to clades that contain potential pathogens. All of them represent organisms that may be relevant for NASA planetary protection (NASA, 2019a), since they may persist on the suit.

      Conclusion

      The microbial characterization of spacesuits examined during this study established the following scenarios. (1) Viable microbes are present on the exterior of spacesuits. (2) Certain microbes exist on suit joints and gauntlets. (3) The microbiomes on suit surfaces resemble those of human skin and airways. More sophisticated approaches can help to conclusively establish whether microbes do migrate from the interior to the exterior of pressurized suits in space.

      Additional work to better determine microbial origin and evaluate microbial contamination mitigation techniques is warranted. This report is a first step in establishing a catalog of microbial sequences known to occur on spacesuits and equipment. Gene specific marker, single nucleotide polymorphism (gene property), single nucleotide variation, and with deep coverage (×100) should be performed to track the source of microbial passage between the exterior and interior of currently existing spacesuits.

      Data Availability Statement

      The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/genbank/, PRJNA545796.

      Author Contributions

      MR developed EVA tool, collected samples, coordinated and designed the study with input from KV. MR, MB, and AR collected samples at JSC. MB and AR wrote the cleaning, sterilization, and assembly protocol for the EVA Swab Kit. GM contributed to sample processing, traditional microbiology assays, extracted DNA, assayed molecular microbial burden, generated corresponding figures, tables, and text associated with these analyses. CM group generated shotgun metagenome sequencing. DD designed, coordinated, and carried out computational analyses. NO’H and RO assisted with strain level analyses. NS conducted metagenome-based phylogenetic and functional analyses and interpreted the results. KV, DD, and MS drafted the manuscript and responsible for data analysis and interpretation. CM edited the manuscript. All authors read and approved the final manuscript.

      Conflict of Interest

      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.

      Publisher’s Note

      All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

      Funding. This research was supported by the JPL Advanced Concept Development fund awarded to KV that funded a fellowship for GM. DD was supported by the Tri-Institutional Training Program in Computational Biology and Medicine (CBM) funded by the NIH grant 1T32GM083937. We thank XSEDE and Philip Blood for their support of this project. AR and MB were supported by the Astromaterials Acquisition and Curation Office at NASA’s Johnson Space Center. EVA Swab Kit development and spacesuit testing was funded via a Johnson Space Center Independent Research and Development award for MR.

      Part of the research described in this publication was carried out at the Jet Propulsion Laboratory (JPL), California Institute of Technology, under a contract with NASA. Government sponsorship acknowledged.

      Supplementary Material

      The Supplementary Material for this article can be found online at: /articles/10.3389/fmicb.2021.608478/full#supplementary-material

      Phylogenetic tree of cultivated strains (bacteria [A] and fungi [B]) isolated from spacesuit surfaces.

      UMAP of taxonomic profiles from all samples colored by whether the samples were treated with PMA or not.

      Metadata for samples including what suit they were taken from and when.

      Bacterial abundances (in number of reads) from metagenomics for all samples.

      Eukaryote abundances (in number of reads) from metagenomics for all samples.

      Percentage of reads assigned to different major groups.

      References Breuker M. McNamara C. Young L. Perry T. Young A. Mitchell R. (2003). Fungal growth on synthetic cloth from Apollo spacesuits. Ann. Microbiol. 53 4754. Buchfink B. Xie C. Huson D. H. (2015). Fast and sensitive protein alignment using DIAMOND. Nat. Methods 12 5960. 10.1038/nmeth.3176 25402007 Buschle-Diller G. Zeronian S. H. Pan N. Yoon M. Y. (1994). Enzymatic hydrolysis of cotton, linen, ramie, and viscose rayon fabrics. Text. Res. J. 64 270279. 10.1177/004051759406400504 Callewaert C. De Maeseneire E. Kerckhof F.-M. Verliefde A. Van de Wiele T. Boon N. (2014). Microbial odor profile of polyester and cotton clothes after a fitness session. Appl. Environ. Microbiol. 80 66116619. 10.1128/AEM.01422-14 25128346 Cappitelli F. Sorlini C. (2008). Microorganisms attack synthetic polymers in items representing our cultural heritage. Appl. Environ. Microbiol. 74 564569. 10.1128/aem.01768-07 18065627 Cataño J. C. Echeverri L. M. Szela C. (2012). Bacterial contamination of clothes and environmental items in a third-level hospital in Colombia. Interdiscipl. Perspect. Infect. Dis. 2012:507640. 10.1155/2012/507640 22536231 Chase J. Fouquier J. Zare M. Sonderegger D. L. Knight R. Kelley S. T. (2016). Geography and location are the primary drivers of office microbiome composition. mSystems 1:e0022-16. 10.1128/mSystems.00022-16 27822521 Checinska A. Probst A. J. Vaishampayan P. White J. R. Kumar D. Stepanov V. G. (2015). Microbiomes of the dust particles collected from the international space station and spacecraft assembly facilities. Microbiome 3:50. 10.1186/s40168-015-0116-3 26502721 Checinska Sielaff A. Urbaniak C. Mohan G. B. M. Stepanov V. G. Tran Q. Wood J. M. (2019). Characterization of the total and viable bacterial and fungal communities associated with the International Space Station surfaces. Microbiome 7:50. 10.1186/s40168-019-0666-x 30955503 Daly D. J. Murphy C. McDermott S. D. (2012). The transfer of touch DNA from hands to glass, fabric and wood. Forens. Sci. Int. Genet. 6 4146. 10.1016/j.fsigen.2010.12.016 21330229 Danko D. Bezdan D. Afshin E. E. Ahsanuddin S. Bhattacharya C. Butler D. J. (2021a). A global metagenomic map of urban microbiomes and antimicrobial resistance. Cell 184 33763393e17. 10.1016/j.cell.2021.05.002 34043940 Danko D. C. Sierra M. A. Benardini J. N. Guan L. Wood J. M. Singh N. (2021b). A comprehensive metagenomics framework to characterize organisms relevant for planetary protection. Microbiome 9 115. Dean C. (2018). Resistome Analyzer. Available online at: https://github.com/cdeanj/resistomeanalyzer Debus A. Arnould J. (2008). Planetary protection issues related to human missions to Mars. Adv. Space Res. 42 11201127. 10.1016/j.asr.2007.10.005 Delcher A. L. Salzberg S. L. Phillippy A. M. (2003). Using MUMmer to identify similar regions in large sequence sets. Curr. Protoc. Bioinform. 10:Unit10.3. Fernandes A. D. Macklaim J. M. Linn T. G. Reid G. Gloor G. B. (2013). ANOVA-like differential expression (ALDEx) analysis for mixed population RNA-Seq. PLoS One 8:e67019. 10.1371/journal.pone.0067019 23843979 Franzosa E. A. McIver L. J. Rahnavard G. Thompson L. R. Schirmer M. Weingart G. (2018). Species-level functional profiling of metagenomes and metatranscriptomes. Nat. Methods 15 962968. 10.1038/s41592-018-0176-y 30377376 Hsu T. Joice R. Vallarino J. Abu-Ali G. Hartmann E. M. Shafquat A. (2016). Urban transit system microbial communities differ by surface type and interaction with humans and the environment. mSystems 1:e0018-16. 10.1128/mSystems.00018-16 27822528 Jensen G. L. (2013). The Human Microbiome, Diet, and Health: Workshop Summary. Washington, DC: National Academies Press. Kang D. D. Li F. Kirton E. Thomas A. Egan R. An H. (2019). MetaBAT 2: an adaptive binning algorithm for robust and efficient genome reconstruction from metagenome assemblies. PeerJ 7:e7359. 10.7717/peerj.7359 31388474 Kettleson E. M. Adhikari A. Vesper S. Coombs K. Indugula R. Reponen T. (2015). Key determinants of the fungal and bacterial microbiomes in homes. Environ. Res. 138 130135. 10.1016/j.envres.2015.02.003 25707017 Kim O.-S. Cho Y.-J. Lee K. Yoon S.-H. Kim M. Na H. (2012). Introducing EzTaxon-e: a prokaryotic 16S rRNA gene sequence database with phylotypes that represent uncultured species. Intern. J. Syst. Evol. Microbiol. 62 716721. 10.1099/ijs.0.038075-0 22140171 Koljalg U. Nilsson R. H. Abarenkov K. Tedersoo L. Taylor A. F. Bahram M. (2013). Towards a unified paradigm for sequence-based identification of fungi. Mol. Ecol. 22 52715277. 10.1111/mec.12481 24112409 Kovaleva J. Degener J. E. van der Mei H. C. (2014). Methylobacterium and its role in health care-associated infection. J. Clin. Microbiol. 52 13171321. 10.1128/jcm.03561-13 24430456 Kwan K. Cooper M. La Duc M. T. Vaishampayan P. Stam C. Benardini J. N. (2011). Evaluation of procedures for the collection, processing, and analysis of biomolecules from low-biomass surfaces. Appl. Environ. Microbiol. 77 29432953. 10.1128/aem.02978-10 21398492 Lai X. Cao L. Tan H. Fang S. Huang Y. Zhou S. (2007). Fungal communities from methane hydrate-bearing deep-sea marine sediments in South China See. ISME J. 1 756762. 10.1038/ismej.2007.51 18059498 Lakin S. M. Dean C. Noyes N. R. Dettenwanger A. Ross A. S. Doster E. (2017). MEGARes: an antimicrobial resistance database for high throughput sequencing. Nucleic Acids Res. 45 D574D580. Lane D. J. (1991). “Nucleic acid techniques in bacterial systematics,” in Nucleic Acid Techniques in Bacterial Systematics, Vol. 1 eds Stackebrandt E. Goodfellow M. (New York, NY: Wiley), 115175. Langmead B. Salzberg S. L. (2012). Fast gapped-read alignment with Bowtie 2. Nat. Methods 9:357. 10.1038/nmeth.1923 22388286 Lax S. Sangwan N. Smith D. Larsen P. Handley K. M. Richardson M. (2017). Bacterial colonization and succession in a newly opened hospital. Sci. Transl. Med. 9:eaah6500. 10.1126/scitranslmed.aah6500 28539477 Lee S.-Y. Woo S.-K. Lee S.-M. Eom Y.-B. (2016). Forensic analysis using microbial community between skin bacteria and fabrics. Toxicol. Environ. Health Sci. 8 263270. 10.1007/s13530-016-0284-y Li D. Liu C.-M. Luo R. Sadakane K. Lam T.-W. (2015). MEGAHIT: an ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph. Bioinformatics 31 16741676. 10.1093/bioinformatics/btv033 25609793 Linacre A. Pekarek V. Swaran Y. C. Tobe S. S. (2010). Generation of DNA profiles from fabrics without DNA extraction. Forens. Sci. Int. Genet. 4 137141. 10.1016/j.fsigen.2009.07.006 20129473 McInnes L. Healy J. Melville J. (2018). Umap: Uniform manifold approximation and projection for dimension reduction. arXiv [Preprint]. Available online at: https://arxiv.org/abs/1802.03426 NASA (2019a). NASA Policy Instruction-8020.7G: NASA Policy on Planetary Protection Requirements for Human Extraterrestrial Missions. Washington, DC: NASA. NASA (2019b). Orion Suit Equipped to Expect the Unexpected on Artemis Missions. Washington, DC: NASA. National Academies of Sciences, Engineering and Medicine (2018). A Midterm Assessment of Implementation of the Decadal Survey on Life and Physical Sciences Research at NASA. Washington, DC: The National Academies Press. Nelson K. E. Weinstock G. M. Highlander S. K. Worley K. C. Creasy H. H. Wortman J. R. (2010). A catalog of reference genomes from the human microbiome. Science 328 994999. Newman D. Schmidt P. Rahn D. (2000). Modeling the Extravehicular Mobility Unit (EMU) Space Suit: Physiological Implications for Extravehicular Activity (EVA) (0148-7191). Available online at: http://web.mit.edu/aeroastro/www/people/dnewman/pdfs/DJN_ICES2000-3.26.pdf Nicholson W. L. Schuerger A. C. Race M. S. (2009). Migrating microbes and planetary protection. Trends Microbiol. 17 389392. 10.1016/j.tim.2009.07.001 19726193 NRC (2014). Pathways to Exploration: Rationales and Approaches for a U.S. Program of Human Space Exploration. Washington, DC: The National Academies Press. Ounit R. Lonardi S. (2016). Higher classification sensitivity of short metagenomic reads with CLARK-S. Bioinformatics 32 38233825. 10.1093/bioinformatics/btw542 27540266 Parks D. H. Imelfort M. Skennerton C. T. Hugenholtz P. Tyson G. W. (2015). CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res. 25 10431055. 10.1101/gr.186072.114 25977477 Roberts W. T. Davidson P. M. (1986). Growth characteristics of selected fungi on polyvinyl chloride film. Appl. Environ. Microbiol. 51 673676. 10.1128/aem.51.4.673-676.1986 3707118 Rucker M. A. Hood D. Walker M. Venkateswaran K. J. Schuerger A. C. (2018). EVA swab tool to support planetary protection and astrobiology evaluations. Paper Presented at the 2018 IEEE Aerospace Conference, New York, NY. Sandle T. (2011). A review of cleanroom microflora: types, trends, and patterns. PDA J. Pharm. Sci. Technol. 65 392403. 10.5731/pdajpst.2011.00765 22293526 Schwartz S. J. Hoffman J. A. Hodgson E. Ronzani P. A. (2002). Is there space for wearables?. Paper Presented at the Proceedings of the Sixth International Symposium on Wearable Computers, London. Shafquat A. Joice R. Simmons S. L. Huttenhower C. (2014). Functional and phylogenetic assembly of microbial communities in the human microbiome. Trends Microbiol. 22 261266. 10.1016/j.tim.2014.01.011 24618403 Singh N. K. Wood J. M. Karouia F. Venkateswaran K. (2018). Succession and persistence of microbial communities and antimicrobial resistance genes associated with International Space Station environmental surfaces. Microbiome 6:214. 10.1186/s40168-018-0609-y 30514368 Sterndorff E. B. Russel J. Jakobsen J. Mortensen M. S. Gori K. Herschend J. (2020). The T-shirt microbiome is distinct between individuals and shaped by washing and fabric type. Environ. Res. 185:109449. 10.1016/j.envres.2020.109449 32278157 Sugita T. Yamazaki T. Makimura K. Cho O. Yamada S. Ohshima H. (2016). Comprehensive analysis of the skin fungal microbiota of astronauts during a half-year stay at the International Space Station. Sabouraudia 54 232239. 10.1093/mmy/myv121 26773135 Szostak-Kotowa J. (2004). Biodeterioration of textiles. Intern. Biodeteriorat. Biodegrad. 53 165170. 10.1016/S0964-8305(03)00090-8 Taylor D. L. Bruns T. D. (1999). Community structure of ectomycorrhizal fungi in a Pinus muricata forest: minimal overlap between the mature forest and resistant propagule communities. Mol. Ecol. 8 18371850. 10.1046/j.1365-294x.1999.00773.x 10620228 Truong D. T. Franzosa E. A. Tickle T. L. Scholz M. Weingart G. Pasolli E. (2015). MetaPhlAn2 for enhanced metagenomic taxonomic profiling. Nat. Methods 12 902903. 10.1038/nmeth.3589 26418763 Turnbaugh P. J. Ley R. E. Hamady M. Fraser-Liggett C. M. Knight R. Gordon J. I. (2007). The human microbiome project. Nature 449 804810. 10.1038/nature06244 17943116 Turner S. Pryer K. M. Miao V. P. Palmer J. D. (1999). Investigating deep phylogenetic relationships among cyanobacteria and plastids by small subunit rRNA sequence analysis. J. Eukaryot. Microbiol. 46 327338. 10.1111/j.1550-7408.1999.tb04612.x 10461381 Vaishampayan P. Probst A. J. La Duc M. T. Bargoma E. Benardini J. N. Andersen G. L. (2013). New perspectives on viable microbial communities in low-biomass cleanroom environments. ISME J. 7 312324. 10.1038/ismej.2012.114 23051695 Venkateswaran K. Hattori N. La Duc M. T. Kern R. (2003). ATP as a biomarker of viable microorganisms in clean-room facilities. J. Microbiol. Methods 52 367377. 10.1016/s0167-7012(02)00192-6 Watson R. D. (2014). Modified Advanced Crew Escape Suit Intravehicular Activity Suit for Extravehicular Activity Mobility Evaluations. Available online at: http://hdl.handle.net/2346/59685 Webb J. S. Nixon M. Eastwood I. M. Greenhalgh M. Robson G. D. Handley P. S. (2000). Fungal colonization and biodeterioration of plasticized polyvinyl chloride. Appl. Environ. Microbiol. 66 31943200. 10.1128/aem.66.8.3194-3200.2000 10919769 Westwood J. Burnett M. Spratt D. Ball M. Wilson D. J. Wellsteed S. (2014). The hospital microbiome project: meeting report for the UK science and innovation network UK-USA workshop ‘beating the superbugs: hospital microbiome studies for tackling antimicrobial resistance’, October 14th 2013. Stand. Genom. Sci. 9:12. 10.1186/1944-3277-9-12
      ‘Oh, my dear Thomas, you haven’t heard the terrible news then?’ she said. ‘I thought you would be sure to have seen it placarded somewhere. Alice went straight to her room, and I haven’t seen her since, though I repeatedly knocked at the door, which she has locked on the inside, and I’m sure it’s most unnatural of her not to let her own mother comfort her. It all happened in a moment: I have always said those great motor-cars shouldn’t be allowed to career about the streets, especially when they are all paved with cobbles as they are at Easton Haven, which are{331} so slippery when it’s wet. He slipped, and it went over him in a moment.’ My thanks were few and awkward, for there still hung to the missive a basting thread, and it was as warm as a nestling bird. I bent low--everybody was emotional in those days--kissed the fragrant thing, thrust it into my bosom, and blushed worse than Camille. "What, the Corner House victim? Is that really a fact?" "My dear child, I don't look upon it in that light at all. The child gave our picturesque friend a certain distinction--'My husband is dead, and this is my only child,' and all that sort of thing. It pays in society." leave them on the steps of a foundling asylum in order to insure [See larger version] Interoffice guff says you're planning definite moves on your own, J. O., and against some opposition. Is the Colonel so poor or so grasping—or what? Albert could not speak, for he felt as if his brains and teeth were rattling about inside his head. The rest of[Pg 188] the family hunched together by the door, the boys gaping idiotically, the girls in tears. "Now you're married." The host was called in, and unlocked a drawer in which they were deposited. The galleyman, with visible reluctance, arrayed himself in the garments, and he was observed to shudder more than once during the investiture of the dead man's apparel. HoME香京julia种子在线播放 ENTER NUMBET 0016lianme.com.cn
      jlgopi.com.cn
      jdzrctc.com.cn
      www.idcqek.com.cn
      imlark.com.cn
      www.kschain.com.cn
      kychain.com.cn
      eekpls.com.cn
      mjdcks.com.cn
      www.wypiano.com.cn
      处女被大鸡巴操 强奸乱伦小说图片 俄罗斯美女爱爱图 调教强奸学生 亚洲女的穴 夜来香图片大全 美女性强奸电影 手机版色中阁 男性人体艺术素描图 16p成人 欧美性爱360 电影区 亚洲电影 欧美电影 经典三级 偷拍自拍 动漫电影 乱伦电影 变态另类 全部电 类似狠狠鲁的网站 黑吊操白逼图片 韩国黄片种子下载 操逼逼逼逼逼 人妻 小说 p 偷拍10幼女自慰 极品淫水很多 黄色做i爱 日本女人人体电影快播看 大福国小 我爱肏屄美女 mmcrwcom 欧美多人性交图片 肥臀乱伦老头舔阴帝 d09a4343000019c5 西欧人体艺术b xxoo激情短片 未成年人的 插泰国人夭图片 第770弾み1 24p 日本美女性 交动态 eee色播 yantasythunder 操无毛少女屄 亚洲图片你懂的女人 鸡巴插姨娘 特级黄 色大片播 左耳影音先锋 冢本友希全集 日本人体艺术绿色 我爱被舔逼 内射 幼 美阴图 喷水妹子高潮迭起 和后妈 操逼 美女吞鸡巴 鸭个自慰 中国女裸名单 操逼肥臀出水换妻 色站裸体义术 中国行上的漏毛美女叫什么 亚洲妹性交图 欧美美女人裸体人艺照 成人色妹妹直播 WWW_JXCT_COM r日本女人性淫乱 大胆人艺体艺图片 女同接吻av 碰碰哥免费自拍打炮 艳舞写真duppid1 88电影街拍视频 日本自拍做爱qvod 实拍美女性爱组图 少女高清av 浙江真实乱伦迅雷 台湾luanlunxiaoshuo 洛克王国宠物排行榜 皇瑟电影yy频道大全 红孩儿连连看 阴毛摄影 大胆美女写真人体艺术摄影 和风骚三个媳妇在家做爱 性爱办公室高清 18p2p木耳 大波撸影音 大鸡巴插嫩穴小说 一剧不超两个黑人 阿姨诱惑我快播 幼香阁千叶县小学生 少女妇女被狗强奸 曰人体妹妹 十二岁性感幼女 超级乱伦qvod 97爱蜜桃ccc336 日本淫妇阴液 av海量资源999 凤凰影视成仁 辰溪四中艳照门照片 先锋模特裸体展示影片 成人片免费看 自拍百度云 肥白老妇女 女爱人体图片 妈妈一女穴 星野美夏 日本少女dachidu 妹子私处人体图片 yinmindahuitang 舔无毛逼影片快播 田莹疑的裸体照片 三级电影影音先锋02222 妻子被外国老头操 观月雏乃泥鳅 韩国成人偷拍自拍图片 强奸5一9岁幼女小说 汤姆影院av图片 妹妹人艺体图 美女大驱 和女友做爱图片自拍p 绫川まどか在线先锋 那么嫩的逼很少见了 小女孩做爱 处女好逼连连看图图 性感美女在家做爱 近距离抽插骚逼逼 黑屌肏金毛屄 日韩av美少女 看喝尿尿小姐日逼色色色网图片 欧美肛交新视频 美女吃逼逼 av30线上免费 伊人在线三级经典 新视觉影院t6090影院 最新淫色电影网址 天龙影院远古手机版 搞老太影院 插进美女的大屁股里 私人影院加盟费用 www258dd 求一部电影里面有一个二猛哥 深肛交 日本萌妹子人体艺术写真图片 插入屄眼 美女的木奶 中文字幕黄色网址影视先锋 九号女神裸 和骚人妻偷情 和潘晓婷做爱 国模大尺度蜜桃 欧美大逼50p 西西人体成人 李宗瑞继母做爱原图物处理 nianhuawang 男鸡巴的视屏 � 97免费色伦电影 好色网成人 大姨子先锋 淫荡巨乳美女教师妈妈 性nuexiaoshuo WWW36YYYCOM 长春继续给力进屋就操小女儿套干破内射对白淫荡 农夫激情社区 日韩无码bt 欧美美女手掰嫩穴图片 日本援交偷拍自拍 入侵者日本在线播放 亚洲白虎偷拍自拍 常州高见泽日屄 寂寞少妇自卫视频 人体露逼图片 多毛外国老太 变态乱轮手机在线 淫荡妈妈和儿子操逼 伦理片大奶少女 看片神器最新登入地址sqvheqi345com账号群 麻美学姐无头 圣诞老人射小妞和强奸小妞动话片 亚洲AV女老师 先锋影音欧美成人资源 33344iucoom zV天堂电影网 宾馆美女打炮视频 色五月丁香五月magnet 嫂子淫乱小说 张歆艺的老公 吃奶男人视频在线播放 欧美色图男女乱伦 avtt2014ccvom 性插色欲香影院 青青草撸死你青青草 99热久久第一时间 激情套图卡通动漫 幼女裸聊做爱口交 日本女人被强奸乱伦 草榴社区快播 2kkk正在播放兽骑 啊不要人家小穴都湿了 www猎奇影视 A片www245vvcomwwwchnrwhmhzcn 搜索宜春院av wwwsee78co 逼奶鸡巴插 好吊日AV在线视频19gancom 熟女伦乱图片小说 日本免费av无码片在线开苞 鲁大妈撸到爆 裸聊官网 德国熟女xxx 新不夜城论坛首页手机 女虐男网址 男女做爱视频华为网盘 激情午夜天亚洲色图 内裤哥mangent 吉沢明歩制服丝袜WWWHHH710COM 屌逼在线试看 人体艺体阿娇艳照 推荐一个可以免费看片的网站如果被QQ拦截请复制链接在其它浏览器打开xxxyyy5comintr2a2cb551573a2b2e 欧美360精品粉红鲍鱼 教师调教第一页 聚美屋精品图 中韩淫乱群交 俄罗斯撸撸片 把鸡巴插进小姨子的阴道 干干AV成人网 aolasoohpnbcn www84ytom 高清大量潮喷www27dyycom 宝贝开心成人 freefronvideos人母 嫩穴成人网gggg29com 逼着舅妈给我口交肛交彩漫画 欧美色色aV88wwwgangguanscom 老太太操逼自拍视频 777亚洲手机在线播放 有没有夫妻3p小说 色列漫画淫女 午间色站导航 欧美成人处女色大图 童颜巨乳亚洲综合 桃色性欲草 色眯眯射逼 无码中文字幕塞外青楼这是一个 狂日美女老师人妻 爱碰网官网 亚洲图片雅蠛蝶 快播35怎么搜片 2000XXXX电影 新谷露性家庭影院 深深候dvd播放 幼齿用英语怎么说 不雅伦理无需播放器 国外淫荡图片 国外网站幼幼嫩网址 成年人就去色色视频快播 我鲁日日鲁老老老我爱 caoshaonvbi 人体艺术avav 性感性色导航 韩国黄色哥来嫖网站 成人网站美逼 淫荡熟妇自拍 欧美色惰图片 北京空姐透明照 狼堡免费av视频 www776eom 亚洲无码av欧美天堂网男人天堂 欧美激情爆操 a片kk266co 色尼姑成人极速在线视频 国语家庭系列 蒋雯雯 越南伦理 色CC伦理影院手机版 99jbbcom 大鸡巴舅妈 国产偷拍自拍淫荡对话视频 少妇春梦射精 开心激动网 自拍偷牌成人 色桃隐 撸狗网性交视频 淫荡的三位老师 伦理电影wwwqiuxia6commqiuxia6com 怡春院分站 丝袜超短裙露脸迅雷下载 色制服电影院 97超碰好吊色男人 yy6080理论在线宅男日韩福利大全 大嫂丝袜 500人群交手机在线 5sav 偷拍熟女吧 口述我和妹妹的欲望 50p电脑版 wwwavtttcon 3p3com 伦理无码片在线看 欧美成人电影图片岛国性爱伦理电影 先锋影音AV成人欧美 我爱好色 淫电影网 WWW19MMCOM 玛丽罗斯3d同人动画h在线看 动漫女孩裸体 超级丝袜美腿乱伦 1919gogo欣赏 大色逼淫色 www就是撸 激情文学网好骚 A级黄片免费 xedd5com 国内的b是黑的 快播美国成年人片黄 av高跟丝袜视频 上原保奈美巨乳女教师在线观看 校园春色都市激情fefegancom 偷窥自拍XXOO 搜索看马操美女 人本女优视频 日日吧淫淫 人妻巨乳影院 美国女子性爱学校 大肥屁股重口味 啪啪啪啊啊啊不要 操碰 japanfreevideoshome国产 亚州淫荡老熟女人体 伦奸毛片免费在线看 天天影视se 樱桃做爱视频 亚卅av在线视频 x奸小说下载 亚洲色图图片在线 217av天堂网 东方在线撸撸-百度 幼幼丝袜集 灰姑娘的姐姐 青青草在线视频观看对华 86papa路con 亚洲1AV 综合图片2区亚洲 美国美女大逼电影 010插插av成人网站 www色comwww821kxwcom 播乐子成人网免费视频在线观看 大炮撸在线影院 ,www4KkKcom 野花鲁最近30部 wwwCC213wapwww2233ww2download 三客优最新地址 母亲让儿子爽的无码视频 全国黄色片子 欧美色图美国十次 超碰在线直播 性感妖娆操 亚洲肉感熟女色图 a片A毛片管看视频 8vaa褋芯屑 333kk 川岛和津实视频 在线母子乱伦对白 妹妹肥逼五月 亚洲美女自拍 老婆在我面前小说 韩国空姐堪比情趣内衣 干小姐综合 淫妻色五月 添骚穴 WM62COM 23456影视播放器 成人午夜剧场 尼姑福利网 AV区亚洲AV欧美AV512qucomwwwc5508com 经典欧美骚妇 震动棒露出 日韩丝袜美臀巨乳在线 av无限吧看 就去干少妇 色艺无间正面是哪集 校园春色我和老师做爱 漫画夜色 天海丽白色吊带 黄色淫荡性虐小说 午夜高清播放器 文20岁女性荫道口图片 热国产热无码热有码 2015小明发布看看算你色 百度云播影视 美女肏屄屄乱轮小说 家族舔阴AV影片 邪恶在线av有码 父女之交 关于处女破处的三级片 极品护士91在线 欧美虐待女人视频的网站 享受老太太的丝袜 aaazhibuo 8dfvodcom成人 真实自拍足交 群交男女猛插逼 妓女爱爱动态 lin35com是什么网站 abp159 亚洲色图偷拍自拍乱伦熟女抠逼自慰 朝国三级篇 淫三国幻想 免费的av小电影网站 日本阿v视频免费按摩师 av750c0m 黄色片操一下 巨乳少女车震在线观看 操逼 免费 囗述情感一乱伦岳母和女婿 WWW_FAMITSU_COM 偷拍中国少妇在公车被操视频 花也真衣论理电影 大鸡鸡插p洞 新片欧美十八岁美少 进击的巨人神thunderftp 西方美女15p 深圳哪里易找到老女人玩视频 在线成人有声小说 365rrr 女尿图片 我和淫荡的小姨做爱 � 做爱技术体照 淫妇性爱 大学生私拍b 第四射狠狠射小说 色中色成人av社区 和小姨子乱伦肛交 wwwppp62com 俄罗斯巨乳人体艺术 骚逼阿娇 汤芳人体图片大胆 大胆人体艺术bb私处 性感大胸骚货 哪个网站幼女的片多 日本美女本子把 色 五月天 婷婷 快播 美女 美穴艺术 色百合电影导航 大鸡巴用力 孙悟空操美少女战士 狠狠撸美女手掰穴图片 古代女子与兽类交 沙耶香套图 激情成人网区 暴风影音av播放 动漫女孩怎么插第3个 mmmpp44 黑木麻衣无码ed2k 淫荡学姐少妇 乱伦操少女屄 高中性爱故事 骚妹妹爱爱图网 韩国模特剪长发 大鸡巴把我逼日了 中国张柏芝做爱片中国张柏芝做爱片中国张柏芝做爱片中国张柏芝做爱片中国张柏芝做爱片 大胆女人下体艺术图片 789sss 影音先锋在线国内情侣野外性事自拍普通话对白 群撸图库 闪现君打阿乐 ady 小说 插入表妹嫩穴小说 推荐成人资源 网络播放器 成人台 149大胆人体艺术 大屌图片 骚美女成人av 春暖花开春色性吧 女亭婷五月 我上了同桌的姐姐 恋夜秀场主播自慰视频 yzppp 屄茎 操屄女图 美女鲍鱼大特写 淫乱的日本人妻山口玲子 偷拍射精图 性感美女人体艺木图片 种马小说完本 免费电影院 骑士福利导航导航网站 骚老婆足交 国产性爱一级电影 欧美免费成人花花性都 欧美大肥妞性爱视频 家庭乱伦网站快播 偷拍自拍国产毛片 金发美女也用大吊来开包 缔D杏那 yentiyishu人体艺术ytys WWWUUKKMCOM 女人露奶 � 苍井空露逼 老荡妇高跟丝袜足交 偷偷和女友的朋友做爱迅雷 做爱七十二尺 朱丹人体合成 麻腾由纪妃 帅哥撸播种子图 鸡巴插逼动态图片 羙国十次啦中文 WWW137AVCOM 神斗片欧美版华语 有气质女人人休艺术 由美老师放屁电影 欧美女人肉肏图片 白虎种子快播 国产自拍90后女孩 美女在床上疯狂嫩b 饭岛爱最后之作 幼幼强奸摸奶 色97成人动漫 两性性爱打鸡巴插逼 新视觉影院4080青苹果影院 嗯好爽插死我了 阴口艺术照 李宗瑞电影qvod38 爆操舅母 亚洲色图七七影院 被大鸡巴操菊花 怡红院肿么了 成人极品影院删除 欧美性爱大图色图强奸乱 欧美女子与狗随便性交 苍井空的bt种子无码 熟女乱伦长篇小说 大色虫 兽交幼女影音先锋播放 44aad be0ca93900121f9b 先锋天耗ばさ无码 欧毛毛女三级黄色片图 干女人黑木耳照 日本美女少妇嫩逼人体艺术 sesechangchang 色屄屄网 久久撸app下载 色图色噜 美女鸡巴大奶 好吊日在线视频在线观看 透明丝袜脚偷拍自拍 中山怡红院菜单 wcwwwcom下载 骑嫂子 亚洲大色妣 成人故事365ahnet 丝袜家庭教mp4 幼交肛交 妹妹撸撸大妈 日本毛爽 caoprom超碰在email 关于中国古代偷窥的黄片 第一会所老熟女下载 wwwhuangsecome 狼人干综合新地址HD播放 变态儿子强奸乱伦图 强奸电影名字 2wwwer37com 日本毛片基地一亚洲AVmzddcxcn 暗黑圣经仙桃影院 37tpcocn 持月真由xfplay 好吊日在线视频三级网 我爱背入李丽珍 电影师傅床戏在线观看 96插妹妹sexsex88com 豪放家庭在线播放 桃花宝典极夜著豆瓜网 安卓系统播放神器 美美网丝袜诱惑 人人干全免费视频xulawyercn av无插件一本道 全国色五月 操逼电影小说网 good在线wwwyuyuelvcom www18avmmd 撸波波影视无插件 伊人幼女成人电影 会看射的图片 小明插看看 全裸美女扒开粉嫩b 国人自拍性交网站 萝莉白丝足交本子 七草ちとせ巨乳视频 摇摇晃晃的成人电影 兰桂坊成社人区小说www68kqcom 舔阴论坛 久撸客一撸客色国内外成人激情在线 明星门 欧美大胆嫩肉穴爽大片 www牛逼插 性吧星云 少妇性奴的屁眼 人体艺术大胆mscbaidu1imgcn 最新久久色色成人版 l女同在线 小泽玛利亚高潮图片搜索 女性裸b图 肛交bt种子 最热门有声小说 人间添春色 春色猜谜字 樱井莉亚钢管舞视频 小泽玛利亚直美6p 能用的h网 还能看的h网 bl动漫h网 开心五月激 东京热401 男色女色第四色酒色网 怎么下载黄色小说 黄色小说小栽 和谐图城 乐乐影院 色哥导航 特色导航 依依社区 爱窝窝在线 色狼谷成人 91porn 包要你射电影 色色3A丝袜 丝袜妹妹淫网 爱色导航(荐) 好男人激情影院 坏哥哥 第七色 色久久 人格分裂 急先锋 撸撸射中文网 第一会所综合社区 91影院老师机 东方成人激情 怼莪影院吹潮 老鸭窝伊人无码不卡无码一本道 av女柳晶电影 91天生爱风流作品 深爱激情小说私房婷婷网 擼奶av 567pao 里番3d一家人野外 上原在线电影 水岛津实透明丝袜 1314酒色 网旧网俺也去 0855影院 在线无码私人影院 搜索 国产自拍 神马dy888午夜伦理达达兔 农民工黄晓婷 日韩裸体黑丝御姐 屈臣氏的燕窝面膜怎么样つぼみ晶エリーの早漏チ○ポ强化合宿 老熟女人性视频 影音先锋 三上悠亚ol 妹妹影院福利片 hhhhhhhhsxo 午夜天堂热的国产 强奸剧场 全裸香蕉视频无码 亚欧伦理视频 秋霞为什么给封了 日本在线视频空天使 日韩成人aⅴ在线 日本日屌日屄导航视频 在线福利视频 日本推油无码av magnet 在线免费视频 樱井梨吮东 日本一本道在线无码DVD 日本性感诱惑美女做爱阴道流水视频 日本一级av 汤姆avtom在线视频 台湾佬中文娱乐线20 阿v播播下载 橙色影院 奴隶少女护士cg视频 汤姆在线影院无码 偷拍宾馆 业面紧急生级访问 色和尚有线 厕所偷拍一族 av女l 公交色狼优酷视频 裸体视频AV 人与兽肉肉网 董美香ol 花井美纱链接 magnet 西瓜影音 亚洲 自拍 日韩女优欧美激情偷拍自拍 亚洲成年人免费视频 荷兰免费成人电影 深喉呕吐XXⅩX 操石榴在线视频 天天色成人免费视频 314hu四虎 涩久免费视频在线观看 成人电影迅雷下载 能看见整个奶子的香蕉影院 水菜丽百度影音 gwaz079百度云 噜死你们资源站 主播走光视频合集迅雷下载 thumbzilla jappen 精品Av 古川伊织star598在线 假面女皇vip在线视频播放 国产自拍迷情校园 啪啪啪公寓漫画 日本阿AV 黄色手机电影 欧美在线Av影院 华裔电击女神91在线 亚洲欧美专区 1日本1000部免费视频 开放90后 波多野结衣 东方 影院av 页面升级紧急访问每天正常更新 4438Xchengeren 老炮色 a k福利电影 色欲影视色天天视频 高老庄aV 259LUXU-683 magnet 手机在线电影 国产区 欧美激情人人操网 国产 偷拍 直播 日韩 国内外激情在线视频网给 站长统计一本道人妻 光棍影院被封 紫竹铃取汁 ftp 狂插空姐嫩 xfplay 丈夫面前 穿靴子伪街 XXOO视频在线免费 大香蕉道久在线播放 电棒漏电嗨过头 充气娃能看下毛和洞吗 夫妻牲交 福利云点墦 yukun瑟妃 疯狂交换女友 国产自拍26页 腐女资源 百度云 日本DVD高清无码视频 偷拍,自拍AV伦理电影 A片小视频福利站。 大奶肥婆自拍偷拍图片 交配伊甸园 超碰在线视频自拍偷拍国产 小热巴91大神 rctd 045 类似于A片 超美大奶大学生美女直播被男友操 男友问 你的衣服怎么脱掉的 亚洲女与黑人群交视频一 在线黄涩 木内美保步兵番号 鸡巴插入欧美美女的b舒服 激情在线国产自拍日韩欧美 国语福利小视频在线观看 作爱小视颍 潮喷合集丝袜无码mp4 做爱的无码高清视频 牛牛精品 伊aⅤ在线观看 savk12 哥哥搞在线播放 在线电一本道影 一级谍片 250pp亚洲情艺中心,88 欧美一本道九色在线一 wwwseavbacom色av吧 cos美女在线 欧美17,18ⅹⅹⅹ视频 自拍嫩逼 小电影在线观看网站 筱田优 贼 水电工 5358x视频 日本69式视频有码 b雪福利导航 韩国女主播19tvclub在线 操逼清晰视频 丝袜美女国产视频网址导航 水菜丽颜射房间 台湾妹中文娱乐网 风吟岛视频 口交 伦理 日本熟妇色五十路免费视频 A级片互舔 川村真矢Av在线观看 亚洲日韩av 色和尚国产自拍 sea8 mp4 aV天堂2018手机在线 免费版国产偷拍a在线播放 狠狠 婷婷 丁香 小视频福利在线观看平台 思妍白衣小仙女被邻居强上 萝莉自拍有水 4484新视觉 永久发布页 977成人影视在线观看 小清新影院在线观 小鸟酱后丝后入百度云 旋风魅影四级 香蕉影院小黄片免费看 性爱直播磁力链接 小骚逼第一色影院 性交流的视频 小雪小视频bd 小视频TV禁看视频 迷奸AV在线看 nba直播 任你在干线 汤姆影院在线视频国产 624u在线播放 成人 一级a做爰片就在线看狐狸视频 小香蕉AV视频 www182、com 腿模简小育 学生做爱视频 秘密搜查官 快播 成人福利网午夜 一级黄色夫妻录像片 直接看的gav久久播放器 国产自拍400首页 sm老爹影院 谁知道隔壁老王网址在线 综合网 123西瓜影音 米奇丁香 人人澡人人漠大学生 色久悠 夜色视频你今天寂寞了吗? 菲菲影视城美国 被抄的影院 变态另类 欧美 成人 国产偷拍自拍在线小说 不用下载安装就能看的吃男人鸡巴视频 插屄视频 大贯杏里播放 wwwhhh50 233若菜奈央 伦理片天海翼秘密搜查官 大香蕉在线万色屋视频 那种漫画小说你懂的 祥仔电影合集一区 那里可以看澳门皇冠酒店a片 色自啪 亚洲aV电影天堂 谷露影院ar toupaizaixian sexbj。com 毕业生 zaixian mianfei 朝桐光视频 成人短视频在线直接观看 陈美霖 沈阳音乐学院 导航女 www26yjjcom 1大尺度视频 开平虐女视频 菅野雪松协和影视在线视频 华人play在线视频bbb 鸡吧操屄视频 多啪啪免费视频 悠草影院 金兰策划网 (969) 橘佑金短视频 国内一极刺激自拍片 日本制服番号大全magnet 成人动漫母系 电脑怎么清理内存 黄色福利1000 dy88午夜 偷拍中学生洗澡磁力链接 花椒相机福利美女视频 站长推荐磁力下载 mp4 三洞轮流插视频 玉兔miki热舞视频 夜生活小视频 爆乳人妖小视频 国内网红主播自拍福利迅雷下载 不用app的裸裸体美女操逼视频 变态SM影片在线观看 草溜影院元气吧 - 百度 - 百度 波推全套视频 国产双飞集合ftp 日本在线AV网 笔国毛片 神马影院女主播是我的邻居 影音资源 激情乱伦电影 799pao 亚洲第一色第一影院 av视频大香蕉 老梁故事汇希斯莱杰 水中人体磁力链接 下载 大香蕉黄片免费看 济南谭崔 避开屏蔽的岛a片 草破福利 要看大鸡巴操小骚逼的人的视频 黑丝少妇影音先锋 欧美巨乳熟女磁力链接 美国黄网站色大全 伦蕉在线久播 极品女厕沟 激情五月bd韩国电影 混血美女自摸和男友激情啪啪自拍诱人呻吟福利视频 人人摸人人妻做人人看 44kknn 娸娸原网 伊人欧美 恋夜影院视频列表安卓青青 57k影院 如果电话亭 avi 插爆骚女精品自拍 青青草在线免费视频1769TV 令人惹火的邻家美眉 影音先锋 真人妹子被捅动态图 男人女人做完爱视频15 表姐合租两人共处一室晚上她竟爬上了我的床 性爱教学视频 北条麻妃bd在线播放版 国产老师和师生 magnet wwwcctv1024 女神自慰 ftp 女同性恋做激情视频 欧美大胆露阴视频 欧美无码影视 好女色在线观看 后入肥臀18p 百度影视屏福利 厕所超碰视频 强奸mp magnet 欧美妹aⅴ免费线上看 2016年妞干网视频 5手机在线福利 超在线最视频 800av:cOm magnet 欧美性爱免播放器在线播放 91大款肥汤的性感美乳90后邻家美眉趴着窗台后入啪啪 秋霞日本毛片网站 cheng ren 在线视频 上原亚衣肛门无码解禁影音先锋 美脚家庭教师在线播放 尤酷伦理片 熟女性生活视频在线观看 欧美av在线播放喷潮 194avav 凤凰AV成人 - 百度 kbb9999 AV片AV在线AV无码 爱爱视频高清免费观看 黄色男女操b视频 观看 18AV清纯视频在线播放平台 成人性爱视频久久操 女性真人生殖系统双性人视频 下身插入b射精视频 明星潜规测视频 mp4 免賛a片直播绪 国内 自己 偷拍 在线 国内真实偷拍 手机在线 国产主播户外勾在线 三桥杏奈高清无码迅雷下载 2五福电影院凸凹频频 男主拿鱼打女主,高宝宝 色哥午夜影院 川村まや痴汉 草溜影院费全过程免费 淫小弟影院在线视频 laohantuiche 啪啪啪喷潮XXOO视频 青娱乐成人国产 蓝沢润 一本道 亚洲青涩中文欧美 神马影院线理论 米娅卡莉法的av 在线福利65535 欧美粉色在线 欧美性受群交视频1在线播放 极品喷奶熟妇在线播放 变态另类无码福利影院92 天津小姐被偷拍 磁力下载 台湾三级电髟全部 丝袜美腿偷拍自拍 偷拍女生性行为图 妻子的乱伦 白虎少妇 肏婶骚屄 外国大妈会阴照片 美少女操屄图片 妹妹自慰11p 操老熟女的b 361美女人体 360电影院樱桃 爱色妹妹亚洲色图 性交卖淫姿势高清图片一级 欧美一黑对二白 大色网无毛一线天 射小妹网站 寂寞穴 西西人体模特苍井空 操的大白逼吧 骚穴让我操 拉好友干女朋友3p