The Lianjiang coast in the East China Sea is a typical subtropical marine ecosystem, and shellfish cultivation occupies almost all of the available tidal flats. Many studies have investigated the effects of shellfish cultivation on benthic organisms and sediments, while the impact of shellfish cultivation on plankton ecosystems is still poorly understood. This study investigated the biogeographical patterns of microeukaryotic communities from Lianjiang coastal waters in four seasons using 18S ribosomal RNA gene amplicon sequencing. Microeukaryotes were mainly comprised of Dinoflagellata, Diatomea, Arthropoda, Ciliophora, Chlorophyta, Protalveolata, Cryptophyceae, and Ochrophyta, and presented significant differences in three habitats (the aquaculture area, confluent area, and offshore area) and four seasons. Similarity percentage analysis revealed that Paracalanus parvus, Heterocapsa rotundata, Bestiolina similis, and five additional key taxa contributed to spatio-temporal differences. Seasonal environmental and spatial factors explained 27.47% of microeukaryotic community variation on average, with 11.11% of the variation shared. Environmental variables, particularly depth, pH, and nitrite concentration, were strongly associated with the microeukaryotic community compositions. The neutral community model further demonstrated that stochastic processes were sufficient in shaping substantial variation in microeukaryotic communities across four seasons, which may reveal the remaining unexplained microeukaryotic community variation. We further divided four seasons into the aquaculture stages and non-aquaculture stages, and speculated that aquaculture activities may increase the dispersal limitation of microeukaryotes in coastal waters, especially for the big bodied-microbes like Arthropoda. The results provide a better understanding of the biogeographical patterns, processes, and mechanisms of microeukaryotic communities near shellfish cultivation.