High-throughput sequencing is a powerful tool used for bivalve symbiosis research, but the largest barrier is the contamination of host DNA. In this work, we assessed the host DNA reduction efficiency, microbial community structure, and microbial diversity of four different sample pre-treatment and DNA extraction methods employed in bivalve gill tissue samples. Metagenomic sequencing showed the average proportions of reads belonging to microorganisms retrieved using PowerSoil DNA extraction kit, pre-treatment with differential centrifugation, pre-treatment with filtration, and HostZERO Microbial DNA kit samples were 2.3 ± 0.6%, 2.5 ± 0.2%, 4.7 ± 1.6%, and 42.6 ± 6.8%, respectively. The microbial DNA was effectively enriched with HostZERO Microbial DNA kit. The microbial communities revealed by amplicon sequencing of the 16S rRNA gene showed the taxonomic biases by using four different pre-treatment and DNA extraction methods. The species diversities of DNA samples extracted with the PowerSoil DNA extraction kit were similar, while lower than DNA samples extracted with HostZERO Microbial DNA kit. The results of this study emphasized the bias of these common methods in bivalve symbionts research and will be helpful to choose a fit-for-purpose microbial enrichment strategy in future research on bivalves or other microbe–invertebrate symbioses.