During the recent decade, high-throughput sequencing (HTS) techniques, in particular, DNA metabarcoding, have facilitated increased detection of biodiversity, including harmful algal bloom (HAB) species. In this study, the presence of HAB species and their appearance patterns were investigated by employing molecular and light microscopy-based monitoring in Tokyo Bay, Japan. The potential co-appearance patterns between the HAB species, as well as with other eukaryotes and prokaryotes were investigated using correlation and association rule-based time-series analysis. In total, 40 unique HAB species were detected, including 12 toxin-producing HAB species previously not reported from the area. More than half of the HAB species were present throughout the sampling season (summer to autumn) and no structuring or succession patterns associated with the environmental conditions could be detected. Statistically significant (p < 0.05, rS ranging from −0.88 to 0.90) associations were found amongst the HAB species and other eukaryotic and prokaryotic species, including genera containing growth-limiting bacteria. However, significant correlations between species differed amongst the years, indicating that variability in environmental conditions between the years may have a stronger influence on the microalgal community structure and interspecies interactions than the variability during the sampling season. The association rule-based time-series analysis allowed the detection of a previously reported negative relationship between Synechococcus sp. and Skeletonema sp. in nature. Overall, the results support the applicability of metabarcoding and HTS-based microalgae monitoring, as it facilitates more precise species identification compared to light microscopy, as well as provides input for investigating potential interactions amongst different species/groups through simultaneous detection of multiple species/genera.