SRX20293960 57_S57 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620112 S57 57_S57 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20293961 58_S58 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620113 S58 58_S58 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20293962 67_S67 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620114 S67 67_S67 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20293963 161_S161 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620115 S161 161_S161 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20293964 162_S162 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620116 S162 162_S162 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20293965 163_S163 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620117 S163 163_S163 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20293966 164_S164 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620118 S164 164_S164 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20293967 165_S165 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620119 S165 165_S165 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20293968 166_S166 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620120 S166 166_S166 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20293969 167_S167 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620121 S167 167_S167 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20293970 168_S168 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620122 S168 168_S168 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20293971 169_S169 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620123 S169 169_S169 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20293972 170_S170 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620124 S170 170_S170 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20293973 68_S68 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620125 S68 68_S68 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20293974 171_S171 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620126 S171 171_S171 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20293975 172_S172 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620127 S172 172_S172 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20293976 173_S173 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620128 S173 173_S173 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20293977 174_S174 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620129 S174 174_S174 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20293978 69_S69 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620130 S69 69_S69 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20293979 70_S70 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620131 S70 70_S70 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20293980 71_S71 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620132 S71 71_S71 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20293981 72_S72 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620133 S72 72_S72 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20293982 73_S73 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620134 S73 73_S73 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20293983 74_S74 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620135 S74 74_S74 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20293984 75_S75 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620136 S75 75_S75 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20293985 76_S76 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620137 S76 76_S76 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20293986 59_S59 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620138 S59 59_S59 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20293987 77_S77 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620139 S77 77_S77 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20293988 78_S78 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620140 S78 78_S78 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20293989 79_S79 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620141 S79 79_S79 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20293990 80_S80 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620142 S80 80_S80 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20293991 81_S81 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620143 S81 81_S81 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20293992 113_S113 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620144 S113 113_S113 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20293993 114_S114 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620145 S114 114_S114 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20293994 115_S115 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620146 S115 115_S115 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20293995 116_S116 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620147 S116 116_S116 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20293996 117_S117 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620148 S117 117_S117 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20293997 118_S118 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620149 S118 118_S118 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20293998 63_S63 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620150 S63 63_S63 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20293999 119_S119 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620151 S119 119_S119 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294000 120_S120 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620152 S120 120_S120 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294001 121_S121 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620153 S121 121_S121 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294002 122_S122 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620154 S122 122_S122 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294003 123_S123 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620155 S123 123_S123 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294004 124_S124 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620156 S124 124_S124 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294005 125_S125 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620157 S125 125_S125 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294006 126_S126 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620158 S126 126_S126 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294007 127_S127 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620159 S127 127_S127 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294008 82_S82 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620160 S82 82_S82 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294009 83_S83 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620161 S83 83_S83 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294010 84_S84 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620162 S84 84_S84 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294011 85_S85 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620163 S85 85_S85 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294012 86_S86 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620164 S86 86_S86 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294013 60_S60 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620165 S60 60_S60 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294014 87_S87 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620166 S87 87_S87 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294015 88_S88 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620167 S88 88_S88 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294016 89_S89 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620168 S89 89_S89 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294017 90_S90 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620169 S90 90_S90 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294018 91_S91 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620170 S91 91_S91 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294019 92_S92 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620171 S92 92_S92 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294020 93_S93 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620172 S93 93_S93 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294021 94_S94 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620173 S94 94_S94 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294022 95_S95 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620174 S95 95_S95 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294023 96_S96 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620175 S96 96_S96 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294024 61_S61 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620176 S61 61_S61 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294025 97_S97 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620177 S97 97_S97 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294026 98_S98 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620178 S98 98_S98 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294027 99_S99 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620179 S99 99_S99 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294028 100_S100 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620180 S100 100_S100 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294029 101_S101 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620181 S101 101_S101 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294030 102_S102 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620183 S102 102_S102 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294031 103_S103 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620182 S103 103_S103 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294032 104_S104 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620184 S104 104_S104 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294033 105_S105 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620185 S105 105_S105 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294034 106_S106 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620187 S106 106_S106 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294035 62_S62 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620186 S62 62_S62 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294036 107_S107 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620188 S107 107_S107 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294037 108_S108 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620189 S108 108_S108 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294038 111_S111 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620190 S111 111_S111 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294039 112_S112 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620191 S112 112_S112 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294040 159_S159 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620192 S159 159_S159 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294041 160_S160 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620193 S160 160_S160 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294042 128_S128 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620194 S128 128_S128 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294043 64_S64 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620195 S64 64_S64 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294044 129_S129 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620196 S129 129_S129 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294045 130_S130 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620198 S130 130_S130 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294046 133_S133 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620197 S133 133_S133 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294047 134_S134 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620199 S134 134_S134 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294048 135_S135 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620201 S135 135_S135 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294049 136_S136 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620200 S136 136_S136 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294050 137_S137 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620202 S137 137_S137 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294051 138_S138 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620203 S138 138_S138 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294052 139_S139 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620204 S139 139_S139 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294053 140_S140 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620205 S140 140_S140 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294054 65_S65 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620206 S65 65_S65 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294055 141_S141 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620207 S141 141_S141 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294056 142_S142 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620208 S142 142_S142 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294057 143_S143 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620209 S143 143_S143 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294058 144_S144 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620211 S144 144_S144 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294059 145_S145 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620210 S145 145_S145 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294060 146_S146 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620212 S146 146_S146 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294061 147_S147 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620213 S147 147_S147 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294062 148_S148 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620214 S148 148_S148 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294063 149_S149 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620215 S149 149_S149 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294064 150_S150 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620216 S150 150_S150 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294065 66_S66 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620217 S66 66_S66 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294066 151_S151 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620218 S151 151_S151 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294067 152_S152 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620219 S152 152_S152 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294068 153_S153 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620220 S153 153_S153 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294069 154_S154 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620221 S154 154_S154 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294070 155_S155 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620222 S155 155_S155 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294071 156_S156 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620223 S156 156_S156 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294072 157_S157 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620224 S157 157_S157 AMPLICON METAGENOMIC PCR Illumina MiSeq SRX20294073 158_S158 SRP437110 bp0 The samples were stored in the laboratory at -80C until DNA extraction. For DNA extraction, each membrane was transferred into Lysing Matrix E tube (MP Biomedicals, LLC, Solon, OH) with addition of 2x volume of tris-saline buffer (1 Tris-l; 0,5 EDTA; 5 NaCl; MQ) and homogenized for 5 min at 50 Hz (TissueLyser LT, Qiagen). Then, membranes were enzymatically digested (lysozyme, proteinase K, SDS in total conc. 1%) and DNA was phenol-chloroform extracted (phenol-chloroform 1:1 v/v; chloroform-isoamyl alcohol 24:1 v/v) as reported previously (Gerasimova et al., 2020) [38]. Purification and desalination of DNA was carried out using DNA Clean and Concentrator-5 kit (Zymo Research, USA). Sequencing was performed using Illumina MiSeq platform and MiSeq Reagent Kit V3 2300 bp (Illumina) according to the Illumina workflow (Illumina protocol, Part #15044223, Rev. B). SRS17620225 S158 158_S158 AMPLICON METAGENOMIC PCR Illumina MiSeq