A representative isolate of each morphotype was identified by 16S rRNA gene sequencing analysis (Fredriksson et al., 2013 (link)). For this, bacterial DNA was extracted using a Quick-DNA Fungal/Bacterial Miniprep Kit (Zymo Research, Irvine, CA, USA). The 16S rRNA gene was PCR-amplified in a 25 μl reaction containing 12.5 μl (1×) Radiant 2× PCR Master Mix, 1 μl (5 ng/μl) gDNA as a template, and 1 μl of previously reported primers 27F (5′-AGAGTTTGATCCTGGCTCAG-3′) and 1492R (5′-GGTTACCTTGTTACGACTT-3′) (each at 10 mM). The PCR amplicons were purified with DNA Clean and Concentrator TM-5 (Zymo Research), and sequencing was carried out at LANGEBIO, Mexico. The resulting sequences were compared against the non-redundant database of GenBank (Bethesda, MD, USA; https://www.ncbi.nlm.nih.gov/). Best hits were compared and the highest sequence score (≥99% of the sequence cover for the comparison with the database of NCBI) was used to identify bacterial species. The obtained partial sequences have been deposited in the NCBI GenBank database.