SRX7112007 A-NE-NV-1-125-1 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622918 A-NE-NV-1-125-1 A-NE-NV-1-125-1 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112008 A-NE-NV-1-125-2 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622919 A-NE-NV-1-125-2 A-NE-NV-1-125-2 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112009 A-NE-NV-5-62-1 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622920 A-NE-NV-5-62-1 A-NE-NV-5-62-1 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112010 NA-NE-NV-5-170-1 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622921 NA-NE-NV-5-170-1 NA-NE-NV-5-170-1 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112011 NA-NE-NV-6-47-2 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622922 NA-NE-NV-6-47-2 NA-NE-NV-6-47-2 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112012 NA-NE-NV-6-47-1 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622923 NA-NE-NV-6-47-1 NA-NE-NV-6-47-1 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112013 NA-NE-NV-7-66-2 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622924 NA-NE-NV-7-66-2 NA-NE-NV-7-66-2 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112014 NA-NE-NV-7-66-1 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622925 NA-NE-NV-7-66-1 NA-NE-NV-7-66-1 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112015 NA-NE-NV-8-144-2 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622926 NA-NE-NV-8-144-2 NA-NE-NV-8-144-2 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112016 NA-NE-NV-8-144-1 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622927 NA-NE-NV-8-144-1 NA-NE-NV-8-144-1 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112017 NA-NE-NV-9-119-2 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622928 NA-NE-NV-9-119-2 NA-NE-NV-9-119-2 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112018 NA-NE-NV-9-119-1 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622929 NA-NE-NV-9-119-1 NA-NE-NV-9-119-1 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112019 A-NE-NV-5-62-2 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622930 A-NE-NV-5-62-2 A-NE-NV-5-62-2 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112020 A-NE-NV-6-35-1 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622931 A-NE-NV-6-35-1 A-NE-NV-6-35-1 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112021 A-NE-NV-6-35-2 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622932 A-NE-NV-6-35-2 A-NE-NV-6-35-2 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112022 A-NE-NV-7-95-1 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622933 A-NE-NV-7-95-1 A-NE-NV-7-95-1 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112023 A-NE-NV-7-95-2 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622934 A-NE-NV-7-95-2 A-NE-NV-7-95-2 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112024 A-NE-NV-8-59-1 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622935 A-NE-NV-8-59-1 A-NE-NV-8-59-1 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112025 A-NE-NV-9-64-1 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622936 A-NE-NV-9-64-1 A-NE-NV-9-64-1 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112026 A-NE-NV-9-64-2 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622937 A-NE-NV-9-64-2 A-NE-NV-9-64-2 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112027 CP1-1 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622938 CP1-1 CP1-1 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112028 A-NE-NV-10-120-1 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622939 A-NE-NV-10-120-1 A-NE-NV-10-120-1 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112029 CP1-10 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622940 CP1-10 CP1-10 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112030 CP1-2 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622941 CP1-2 CP1-2 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112031 CP1-3 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622942 CP1-3 CP1-3 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112032 CP1-4 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622943 CP1-4 CP1-4 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112033 CP1-5 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622944 CP1-5 CP1-5 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112034 CP1-6 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622945 CP1-6 CP1-6 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112035 CP1-7 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622946 CP1-7 CP1-7 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112036 CP1-8 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622947 CP1-8 CP1-8 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112037 CP1-9 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622948 CP1-9 CP1-9 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112038 CP2-1 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622949 CP2-1 CP2-1 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112039 A-NE-NV-10-120-2 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622950 A-NE-NV-10-120-2 A-NE-NV-10-120-2 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112040 CP2-10 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622951 CP2-10 CP2-10 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112041 CP2-2 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622952 CP2-2 CP2-2 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112042 CP2-3 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622953 CP2-3 CP2-3 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112043 CP2-4 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622954 CP2-4 CP2-4 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112044 CP2-5 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622955 CP2-5 CP2-5 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112045 CP2-6 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622956 CP2-6 CP2-6 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112046 CP2-7 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622957 CP2-7 CP2-7 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112047 CP2-8 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622958 CP2-8 CP2-8 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112048 CP2-9 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622959 CP2-9 CP2-9 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112049 M11A-1 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622960 M11A-1 M11A-1 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112050 A-NE-NV-2-40-1 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622961 A-NE-NV-2-40-1 A-NE-NV-2-40-1 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112051 M11A-2 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622962 M11A-2 M11A-2 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112052 M11A-3 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622963 M11A-3 M11A-3 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112053 M11A-4 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622964 M11A-4 M11A-4 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112054 M11A-5 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622965 M11A-5 M11A-5 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112055 M13NA-1 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622966 M13NA-1 M13NA-1 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112056 M13NA-2 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622967 M13NA-2 M13NA-2 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112057 M13NA-3 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622968 M13NA-3 M13NA-3 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112058 M13NA-4 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622969 M13NA-4 M13NA-4 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112059 M13NA-5 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622970 M13NA-5 M13NA-5 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112060 M14NA-1 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622971 M14NA-1 M14NA-1 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112061 A-NE-NV-2-40-2 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622972 A-NE-NV-2-40-2 A-NE-NV-2-40-2 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112062 M14NA-2 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622973 M14NA-2 M14NA-2 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112063 M14NA-3 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622974 M14NA-3 M14NA-3 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112064 M14NA-4 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622975 M14NA-4 M14NA-4 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112065 M14NA-5 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622976 M14NA-5 M14NA-5 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112066 M24NA-1 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622977 M24NA-1 M24NA-1 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112067 M24NA-2 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622978 M24NA-2 M24NA-2 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112068 M24NA-3 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622979 M24NA-3 M24NA-3 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112069 M24NA-4 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622980 M24NA-4 M24NA-4 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112070 M24NA-5 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622981 M24NA-5 M24NA-5 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112071 M28A-1 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622982 M28A-1 M28A-1 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112072 A-NE-NV-3-3-1 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622983 A-NE-NV-3-3-1 A-NE-NV-3-3-1 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112073 M28A-2 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622984 M28A-2 M28A-2 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112074 M28A-3 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622985 M28A-3 M28A-3 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112075 M28A-4 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622986 M28A-4 M28A-4 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112076 M28A-5 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622987 M28A-5 M28A-5 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112077 M32A-1 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622988 M32A-1 M32A-1 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112078 M32A-2 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622989 M32A-2 M32A-2 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112079 M32A-3 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622990 M32A-3 M32A-3 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112080 M32A-4 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622991 M32A-4 M32A-4 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112081 M32A-5 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622992 M32A-5 M32A-5 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112082 M33A-1 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622993 M33A-1 M33A-1 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112083 A-NE-NV-3-3-2 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622994 A-NE-NV-3-3-2 A-NE-NV-3-3-2 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112084 M33A-2 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622995 M33A-2 M33A-2 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112085 M33A-3 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622996 M33A-3 M33A-3 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112086 M33A-4 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622997 M33A-4 M33A-4 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112087 M33A-5 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622998 M33A-5 M33A-5 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112088 M38A-1 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5622999 M38A-1 M38A-1 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112089 M38A-2 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5623000 M38A-2 M38A-2 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112090 M38A-3 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5623001 M38A-3 M38A-3 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112091 M38A-4 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5623002 M38A-4 M38A-4 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112092 M38A-5 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5623003 M38A-5 M38A-5 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112093 M39NA-1 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5623004 M39NA-1 M39NA-1 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112094 A-NE-NV-4-51-1 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5623005 A-NE-NV-4-51-1 A-NE-NV-4-51-1 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112095 M39NA-2 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5623006 M39NA-2 M39NA-2 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112096 M39NA-3 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5623007 M39NA-3 M39NA-3 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112097 M39NA-4 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5623008 M39NA-4 M39NA-4 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112098 M39NA-5 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5623009 M39NA-5 M39NA-5 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112099 M7NA-1 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5623010 M7NA-1 M7NA-1 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112100 M7NA-2 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5623011 M7NA-2 M7NA-2 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112101 M7NA-3 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5623012 M7NA-3 M7NA-3 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112102 M7NA-4 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5623013 M7NA-4 M7NA-4 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112103 M7NA-5 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5623014 M7NA-5 M7NA-5 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112104 NA-NE-NV-1-102-1 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5623015 NA-NE-NV-1-102-1 NA-NE-NV-1-102-1 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112105 A-NE-NV-4-51-2 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5623016 A-NE-NV-4-51-2 A-NE-NV-4-51-2 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112106 NA-NE-NV-1-102-2 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5623017 NA-NE-NV-1-102-2 NA-NE-NV-1-102-2 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112107 NA-NE-NV-10-90-2 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5623018 NA-NE-NV-10-90-2 NA-NE-NV-10-90-2 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112108 NA-NE-NV-10-90-1 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5623019 NA-NE-NV-10-90-1 NA-NE-NV-10-90-1 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112109 NA-NE-NV-2-53-2 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5623020 NA-NE-NV-2-53-2 NA-NE-NV-2-53-2 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112110 NA-NE-NV-2-53-1 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5623021 NA-NE-NV-2-53-1 NA-NE-NV-2-53-1 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112111 NA-NE-NV-3-37-2 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5623022 NA-NE-NV-3-37-2 NA-NE-NV-3-37-2 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112112 NA-NE-NV-3-37-1 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5623023 NA-NE-NV-3-37-1 NA-NE-NV-3-37-1 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112113 NA-NE-NV-4-171-2 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5623024 NA-NE-NV-4-171-2 NA-NE-NV-4-171-2 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112114 NA-NE-NV-4-171-1 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5623025 NA-NE-NV-4-171-1 NA-NE-NV-4-171-1 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX7112115 NA-NE-NV-5-170-2 SRP128901 PRJNA429454 We isolated DNA from skin swabs using the PowerSoil DNA Isolation Kit (MoBio Labora- tories Inc., Carlsbad, CA) following the protocol described in Herna ndez-Go mez et al. (2017a). In brief, we implemented two sequential PCRs to prepare our 16S rRNA V2 region amplicon sequencing library using the primer pair 27F/338R (Fierer et al., 2008) and following the methodology described in Herna ndez- Go mez et al. (2017b). We shipped the barcoded sample pool on dry ice overnight to the Cornell University Biotechnology Resource Center (Ithaca, NY). We sequenced the sample pool on a MiSeq machine (Illumina, Inc. San Diego, CA) using the reagent kit V2 to produce 250 base pair (bp) paired-end reads. SRS5623026 NA-NE-NV-5-170-2 NA-NE-NV-5-170-2 AMPLICON METAGENOMIC PCR Illumina MiSeq